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Abstrac  t 


O  Applied  network  theory  and  marginal  analysis  concepts 
were  utilized  to  design,  computerize,  verify,  and  evaluate, 
three  major  software  modifications  to  the  Network  Repair 
Level  Analysis  (NRLA)  model.  First,  a  preprocessor 
subroutine  using  marginal  analysis  techniques  was  developed 
and  tested  to  reduce  the  computer  processing  requirements  of 
the  program.  Second,  a  new  network  labeling  algorithm  which 
solve  the  max-flow  min-cut  problem  is  presented.  This 
algorithm  performs  100  times  faster  than  the  current 
algorithm  and  73  times  faster  than  the  highly  efficient, 
commercially  available,  primal  networking  code  known  as  GNET . 
Third,  for  the  first  time  a  networking  structure  has  been 
designed  which  allows  for  the  inclusion  of  Centralized 
Intermediate  Repair  Facilities  (CIRF)  in  the  repair  level 
analysis  decision  process. , 

These  products  greatly  expand  the  NRLA  model's  capability 
while  at  the  same  time  improving  its  operational  efficiency. 
Through  their  integration  and  use,  System  Program  Managers 
have  a  comprehensive  analytical  tool  to  effectively  conduct 
repair  level  analysis  and  to  design  more  cost-effective 
logisitical  structures  to  support  the  operation  of  Air  Force 
systems. 


The  NRLA  program  is  hosted  on  the  CREATE  Operating  System 


and  contains  approx i amatel y  5500  lines  of  computer  code.  It 
consists  of  a  main  routine  and  twelve  major  subroutines.  The 


Introduct i on 


Every  program  manager,  who  is  responsible  for  acquiring 
a  new  weapon  system  to  be  used  by  tomorrow/s  United  States 
Air  Force  <USAF>,  faces  a  common  challenge.  It  consists  of 
three,  sometimes  conflicting,  functional  goals;  system 
performance,  system  cost,  and  acquisition  schedule.  His  goal 
is  a  formidable  one.  He  must  attempt  to  maximize  the  perfor¬ 
mance  of  the  system  at  the  least  possible  cost  to  the 
government  while  maintaining  the  scheduled  time-frame  for 
procurement  of  the  system.  A  primary  ingredient  in  the 
successful  accomplishment  of  this  effort  is  an  accurate 
assessment  of  the  "operational  support  requirements  and  limi¬ 
tations  of  the  system"  (13:3)  and  the  timely  consideration  of 
alternative  maintenance  concepts  capable  of  achieving  a 
desired  level  of  system  effectiveness. 

During  the  Vietnam  conflict,  the  importance  of 
developing  a  valid  maintenance  plan  as  an  integral  part  of 
the  system's  engineering  development  became  brutally 
apparent.  Many  times  new  systems  were  deployed  to  meet  the 
changing  threat  only  to  have  their  components  subsequently 
fail  causing  excessive  downtime  for  the  system.  An  example 
of  this  was  the  AN/TRC-87  UHF  radio  set.  The  design  of  this 
radio  set  did  not  anticipate  its  extensive  use  in  a  jungle 
environment,  as  a  result  major  component  failures  occured 


which  required  depot  repair  <1:134).  The  availability, 
dependability,  and  capabi I i ty  of  these  systems  were  ail  sub¬ 
stantially  reduced  because  of  a  lack  of  adequate  attention 
being  given  to  their  reliability,  the  resources  necessary  to 
support  their  repair,  and  specific  identification  of  the 
repair  location  for  their  components.  As  a  direct  result  of 
these  experiences,  the  Department  of  Defense  (DOD)  community 
and  the  Air  Force  have  taken  great  strides  to  ensure  that 
such  decisions  become  an  integral  part  of  the  system's  de¬ 
sign.  The  Integrated  Logistic  Support  <ILS>  Program,  estab¬ 
lished  in  1972  by  AFR  800-8,  requires  that  planning  for  cost- 
effective  logistic  support,  including  analysis  of  repaii — 
level  alternatives,  be  conducted  during  the  acquisition  phase 
of  a  weapon  system  <11:3).  In  conjunction  with  this  guidance, 
•AFSC/AFLCR  800-28  establishes  Air  Force  policies  and 
procedures  with  respect  to  Repair  Level  Analysis  < RLA) ■ 
<14:4).  These  policies  and  procedures  are  designed  to  ensure 
that  alternative  maintenance  concepts  are  considered  during 
conceptual  development  and  at  the  appropriate  times  in  the 
life  cycle  of  the  system.  Additionally,  AFSC/AFLC  Regulation 
800-28  identifies  RLA  as  a  separate  evaluation  factor  during 
the  source  selection  process.  This  is  a  further  indication 
of  the  emphasis  and  recognition  that  is  now  being  given 
repair  level  decision  making  <13:3). 

The  Repair  Level  Anal-sis  process  is  the  procedure  by 
which  economic  comparisons  are  made  between  repair  locations 
to  develop  a  comprehensive  and  cost  effective  maintenance 
plan.  RLA  encompasses  a  wide  variety  of  analytical  techni- 


ques  and  methods.  These  techniques  can  be  used  separately  or 
in  combination  to  provide  the  system  program  offices  with 
economically  based  management  information  as  to  the  "best" 
set  of  repair  level  decisions  for  their  system.  In  an  effort 
to  document  the  strengths  and  weaknesses  of  these  different 
methods  and  to  assist  logistics  analysts  in  determining  the 
appropriate  technique  for  a  particular  application,  the  Air 
Force  has  recently  updated  these  procedures  in  AFLC/AFSC 
Pamphlet  800-4,  Repair  Level  Analysis  Procedures  (revised 
June  1983).  Formerly  known  as  Optimum  Repair  Level  Analysis 
(ORLA)  (12),  RLA  now  encompasses  the  following  methods  (13:5) 


(1)  Network  Repair  Level  Analysis  (NRLA) 

(2)  I  tern  Repair  Level  Analysis  (IRLA) 

(3)  Marginal  Analysis  Repair  Level  Analysis  (MARLA) 

(4)  Equal  Cost  Curves  (ECC) 

(5)  SE/Pipeline  Ratio 


Each  of  these  techniques  are  clearly  explained  in  AFLCP/AFSCP 
800-4.  However,  the  NRLA  conceptual  model,  which  were  intro¬ 
duced  in  1980,  are  recognized  as  the  most  comprehensive  ap¬ 
proach  to  repair  level  analysis  and  has  consequently  received 
the  widest  use  and  attention  (13:5).  Since  NRLA's  imple¬ 
mentation,  a  question  has  arisen  as  to  its  ability  to 
continue  to  provide  reliable  and  consistent  repair  level 
logistics  cost  effectiveness  and  ♦•hs  operational  capability 
of  future  systems.  This  research  has  been  initiated  in  an 
effort  to  answer  this  question.  The  remainder  of  this  chapter 
will  deal  with  the  scope  and  limitations  of  the  NRLA  model 
with  an  emphasis  on  identifying  software  enhancements  to 
correct  identified  shortcomings. 


Background 

Prior  to  1980  repair  level  analysis  tor  Air  Force  sys¬ 
tems  was  conducted  on  an  item  by  item  basis  using  a  method 
called  Optimum  Repair  Level  Analysis  (ORLA)  (12:1).  In  ORLA 
each  Line  Replaceable  Unit  (LRU)  and  Shop  Replaceable  Unit 
(SRU)  in  the  system  would  be  separately  priced  tor  each  ot 
three  repair  options;  <1)  repair  at  intermediate  base,  (2) 
repair  at  depot  and  (3)  discard  or  scrap.  The  repair  option 
associated  with  the  minimum  cost  would  be  selected  tor  that 
particular  LRU  or  SRU.  The' problem  with  this  method  ot 
determining  repair  level  decisions  was  that  the  SE  costs 
associated  with  each  LRU/SRU  repair  were  prorated  and 
estimated  according  to  that  LRU/SRUs  use  ot  the  SE  resource. 
Thus  a  LRU/SRU  which  required  20  percent  ot  the  SE's 
available  repair  time  would  need  to  economically  justify  at 
least  20  percent  of  that  SE  cost.  In  addition  if  all 
LRU/SRUs  requiring  that  SE  were  not  chosen  for  the  same 
repair  level  then  the  percent  utilization  and  prorated  costs 
would  need  to  be  recomputed  based  on  the  new  set  of  LRU/SRUs 
at  that  level.  This  could  then  lead  to  other  LRU/SRUs  not 
being  able  to  support  their  required  SE  costs  which  would 
lead  to  further  proration  of  costs  and  possibly  to  all 
LRU/SRUs  being  repaired  at  the  depot  level  when  in  fact  the 
total  set  of  SE  related  LRU/SRUs  could  easily  justify  base 
level  repair  as  the  most  economical  option  (8:6-7).  A  second 
problem  with  the  ORLA  method  of  determining  repair  level 
decisions  was  that  the  Air  force  had  no  official  computer 


program  to  implement  ORLA.  This  resulted  in  virtual  chaos 
when  evaluating  contractor  estimates,  in  that  contractors  did 
not  use  standardized  software  to  conduct  their  analysis. 

To  solve  these  problems  the  NRLA  program  was  developed 
by  the  Concepts  and  Analysis  division,  XRS/AFALC(AFLC)  giving 
the  Air  force,  for  the  first  time,  the  capability  to  use  a 
systems  approach  in  the  development  and  formulation  of  repair 
level  decisions.  NRLA  was  developed  as  a  FORTRAN  based 
software  package  that  is  transportable  between  computer  sys¬ 
tems  and  is  therefore  useable  by  gove'rnment  contractors  pro¬ 
viding  a  consistent  framework  for  evaluation  of  their 
proposal s. 


NRLA 

The  formulation  of  the  repair  level  analysis  problem  as  a 
network  provides  a  couple  of  advantages  over  previous  RLA 
methodologies.  It  takes  into  account  the  LRU  and  SRU 
indenture  relationships  which  keeps  it  from  making 
inconsistent  decisions.  An  example  of  this  would  be  deciding 
to  scrap  a  LRU  but  repair  its  indentured  SRUs.  The  model 
also  treats  each  piece  of  support  equipment  as  a  common 
resource  which  is  shared  by  a  group  of  LRUs  and/or  SRUs.  The 
repair  level  decisions  are  then  made,  based  on  these  LRU/SRU 
group's  ability  to  economically  support  the  purchase  of  the 
support  equipment  for  a  particular  repair  level.  All  of  the 
repair  decisions  are  determined  simultaneously  for  all  of 
the  failure  modes  of  a  group  of  LRUs  and  their  associated 
SRUs,  thus  the  decisions  that  are  made  are  optimal  for  the 


entire  group  of  items. 

The  actual  formulation  of  the  network  is  shown  in  Figure 
1.  The  nodes  of  the  network  represent  the  LRU  failure  modes, 
the  SRUs,  and  the  SE  resources  required  for  the  LRU/SRU 
repairs.  The  arcs  of  the  network  represent  the  costs 
associated  with  the  different  repair  level  options.  Table  I 
defines  the  specific  arc  costs  for  the  network  in  terms  of 
eleven  types  of  logistical  costs.  The  NRLA  Users  Guide 
available  from  XRS/AFALC<AFLC>  details  the  computation  of 
these  costs  in  terms  of  the  LRU/SRU' s  cost,  MTBF,  and 
standard  maintenance  and  supply  factors.  The  heavy  arcs  are 
dummy  arcs  which  provide  the  LRU/SRU/SE  interdependency  rela¬ 
tionships.  The  capacities  of  the  dummy  arcs  are  set  such 
that  they  will  never  enter  into  the  solution  set. 

The  repair  level  decisions  are  obtained  by  applying  an 
optimization  algorithm  to  the  network  to  solve  the  maximin 
flow  minimum  cut  problem.  This  algorithm  provides  a  network 
solution  with  a  cut  set  that  describes  the  unique  optimum  set 
of  minimum  cost  repair  level  decisions  for  the  entire  system. 
It  should  be  pointed  out  that  the  NRLA  model  in  no  way 
attempts  to  compute  a  total  life  cycle  cost  for  a  particular 
set  of  repair  decisions  but,  'includes  only  those  costs  which 
directly  impact  the  repair  level  decision*.  Examples  of 
costs  which  are  not  included  are:  repair  in  place  costs,  and 
the  costs  associated  with  removing  the  failed  LRU/SRU  from 


LOGISTIC  FACTORS 


1.  Support  Equipment 

(la)  Acqu i si  si t i on 

(lb)  Ops  &  Mai nt 

( lc)  Fac i 1 i t i es 

2.  Tech.  Data  Acquis. 

3.  Maint  Training 

4.  Repair  Labor 

5.  Item  Entry 

6.  Supply  Admin. 

7.  Repair  Material 

8.  Packing  &  Shipping 

9.  Base  Spares  Quantity 

10.  Depot  Spares  Quant  it 

11.  Replacement  Spares 


Virtually  any  model  makes  some  simplifying  assumptions 
which  enables  it  to  calculate  a  solution  or  make  some  deter¬ 


mination/conclusion  about  the  system  being  studied.  The  NRLA 
model  is  no  exception  to  this  rule.  The  following 
assumptions  have  been  identified  as  relevant  to  the  solutions 
derived  by  NRLA  (8:3-6). 

<1)  The  user  specifies  the  number  of  bases  and  the  number 
of  end  items  per  base  (assumed  to  be  equal  for 
all  bases).  All  depot  repair  of  a  particular  LRU/SRU 
is  accomplished  at  the  same  depot  location. 

(2)  Base  level  maintenance  system  data  are  equal  for  all 
bases  and  all  types  of  repair  tasks.  The  correspond¬ 
ing  depot  data  factors  are  constant  also. 

(3)  Supply  system  data  factors  are  constant  for  all  LRUs 
and  SRUs  analyzed,  ie.  shipping  times  from  depot  to 
any  CONUS  location  are  equal,  and  the  same  is  true 
for  shipments  from  depot  to  any  overseas  location. 

(4)  Only  a  single  set  of  technical  data  is  purchased 
from  the  contractor. 

(5)  Preventative  and  scheduled  maintenance  actions  are 
not  addressed  by  the  model . 

(6)  The  model  explicitly  evaluates  each  LRU  failure  mode 
for  a  repair  level  decision,  however  SRU  failure 
modes  are  assumed  to  be  similar  enough  to  allow  for 
considering  them  all  equal. 

(7)  Maintenance  man-hours  for  repair  and  SE  utilization 
are  assumed  to  be  equal . 

(8)  The  depot  stock  level  of  SRUs  is  designed  to  satisfy 
base  level  demands.  The  stock  level  supports  base 
level  SRU  remove  and  replace  maintenance  actions 
but  not  similar  depot  actions. 


LifflLU&jff.H.* 

The  NRLA  program  as  it  exists  today  has  proven  to  be  a 
valuable  tool  in  such  areas  as  determining  optimal  repair 


levels  in  support  o-f  the  provisioning  process,  identification 
and  justification  of  LRU/SRU  support  equipment  requirements, 
and  assessment  of  a  system's  sensitivity  to  cost  and  reliabi- 
1 i ty  growth.  However,  limitations  in  the  programs  efficiency 
and  capability  have  resulted  in  the  need  to  enhance  and  in 
some  cases  expand  the  software  of  the  program  to  accommodate 
larger  and  more  complex  systems  and  to  be  able  to  evaluate 
new  maintenance  concepts. 

Approx imatel y  50-60  contractors  and  system  program  of- 
fices  are  currently  using  the  NRLA  model  as  implemented  in 
1980.  User  comments  indicate  that  approx  imatel  y  60'/.  of  the 
analyses  that  they  perform  using  NRLA  are  conducted  on  a 
piecemeal  basis  with  subsystems  being  analyzed  individually. 
This  is  due  to  the  fact  that  it  is  easier  and  more 
convenient  to  analyze  and  work  with  sub-systems.  However, 
the  main  reason  which  drives  users  to  analyze  smaller 
sections  of  the  system  is  computer  resources.  As  the  system 
being  analyzed  gets  larger  and  more  complex  it  requires  more 
computer  storage  space  to  load  the  program  and  more  computer 
time  to  run  it.  Both  storage  requirements  and  run-time 
requirements  are  used  to  prioritize  computer  jobs  with  the 
result  that  as  the  job  gets  bigger  the  turn  around  time  on 
the  analysis  gets  longer,  an  undesireable  result.  Also  once 
the  analysis  of  the  subsystems  is  complete  these  results  are 
then  manually  cross-ref erenced  to  determine  the  final  systems 
set  of  repair  level  decisions. 

A  second  limitation  in  the  NRLA  model  results  from  the 


I 
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•volution  of  •  new  maintenance  concept,  Centralized 
|  Intermediate  Repair  Facilities  (CIRF).  This  concept  is  a 

result  of  the  higher  LRU/SRU/SE  costs  associated  with  new 
systems.  In  many  instances  the  SE  costs  are  too  high  to 
|  allow  for  base  level  repair  and  at  the  same  time  the  pipe¬ 

line  costs  associated  with  depot  repair  are  also  much  higher 
than  is  desired.  CIRF  provides  for  a  single  central  repair 
I  facility  servicing  several  bases  thereby  eliminating  some  of 

the  costs  associated  with  each  of  the  above  repair  options. 
The  NRLA  model  currently  does  not  allow  for  consideration  of 
a  CIRF  maintenance  concept  in  its  optimum  repair  level 
anal ysi s. 

Thesis  Objectives 

The  main  objectives  of  the  thesis  effort  deal  with 
developing  enhancements  to  the  NRLA  model  which  will  provide 
solutions  to  the  shortcomings  as  identified  by  the  users;  ie 

<1>  The  excessive  computer  time  necessary  to  execute  the 
program  for  large  systems 

(2)  Inability  of  the  program  to  consider  the  option  of  a 
Centralized  Intermediate  Repair  Facility  (CIRF) 

The  objectives  fall  into  two  distinct  areas;  (1)  enhancements 
which  will  allow  the  NRLA  model  to  analyze  larger  and  more 
complex  systems  and  (2)  enhancements  which  will  allow  the 
NRLA  model  to  evaluate  additional  maintenance  concepts,  pri¬ 
marily  CIRF.  The  first  objective  area  can  be  further  broken 


down  into  two 


(1)  enhancements  which  increase  the 
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eff  i  c  i  ency  o-f  the  algorithm  which  solves  the  max-flow  min-cut 
problem  and  (2)  enhancements  which  reduce  the  size  of  the 
problem  input  to  the  NRLA  model.  A  short  description  of  each 
of  these  objective  areas  and  the  proposed  solution 
methodology  follows.  A  comprehensive  treatment  of  these 
areas  is  contained  in  Chapter  III. 

The  computational  time  requirements  for  the  solution  of 
the  repair  level  problem  by  the  NRLA  model  are  driven  by  two 
major  factors,  one  internal  to  the  program  and  one  external. 
The  internal  factor  deals  with  the  computational  efficiency 
of  the  labelling  algorithm  subroutine  used  to  solve  the  max- 
flow  min-cut  problem.  Labeling  in  the  context  of  network 
flows  is  the  means  by  which  the  algorithm  identifies  which 
nodes  in  the  network  it  has  visited.  The  labels  prevent  the 
algorithm  from  revisiting  a  node  thus  ensuring  that  the  flow 
augmentation  path  will  be  identified  in  a  finite  number  of 
steps.  The  external  factor  affecting  the  computation  time  is 
the  size  of  the  system  to  be  analyzed  in  terms  of  LRUs,  SRUs, 
and  SE.  Each  of  these  areas  is  addressed  in  the  thesis 
effort . 

The  current  labeling  algorithm  utilized  in  the  NRLA 
model  is  one  developed  by  Ford  and  Fulkerson  in  the  early 
1960s  (17).  The  approach  taken  is  to  perform  a  breadth  first 
search  with  labeling  as  you  proceed  to  the  network  sink. 
Once  the  sink  is  reached  the  flow  is  augmented  on  the  path 
and  the  labels  erased  and  the  process  started  over.  In  a 
large  system  the  majority  of  the  processing  time  is  used 


performing  the  labeling  process  which  can  be  very  ineffi¬ 
cient.  The  approach  taken  to  reduce  this  processing  time  is 
to  use  a  depth  first  search,  as  conceptually  described  by 
Horowitz  and  Sahni  <19:268),  for  flow  augmentation  until  the 
maximum  flow  is  found.  This  is  then  followed  by  a  breadth 
first  search  to  provide  the  correct  labels  for  the  identifi¬ 
cation  of  the  unique  minimum  cost  set  of  repair  level 
dec i si ons. 

In  addition  to  the  depth  first  approach  for  increasing 
the  algorithm  efficiency  a  second  method  of  solving  the 
problem  using  a  large  scale  commercially  available  algorithm 
called  GNET  will  be  explored.  This  method  is  expected  to 
provide  a  viable  alternative  for  the  solution  of  systems  of 
extreme  size,  such  as  the  Bi  bomber  and  the  MX  missile  The 
GNET  algorithm  is  a  primal  algorithm  which  has  been  used  to 
solve  extremely  large  transportation  and  transshipment  prob¬ 
lems  (7).  The  NRLA  program  will  be  modified  to  include  GNET 
as  a  subroutine  to  solve  the  max-flow  problem,  this  will  once 
again  be  followed  by  the  Ford  Fulkerson  labeling  process  to 
identify  the  unique  min  cost  solution  needed  by  NRLA. 

Creating  a  mechanism  to  control  the  second  factor,  that 
of  system  size,  was  initially  much  more  difficult  to  concep¬ 
tualize.  Normally,  the  configuration  of  a  proposed  weapon 
system  does  not  dramatically  change  once  the  preliminary 
design  review  <PDR)  between  the  SPO  and  the  contractor  is 
conducted.  Therefore  the  types  and  number  of  items  (LRUs, 
SRUs,etc.)  which  make-up  the  total  system  are  given  as  fixed 


inputs  to  the  program.  However,  a  possible  solution  to  this 
dilemma  is  the  use  of  another  RLA  technique  called  Marginal 
Repair  Level  Analysis  < MARLA > .  Simply  stated,  the  approach  is 
to  reduce  the  number  of  items  which  have  to  be  analyzed  by 
NRLA.  This  will  be  accomplished  by  incorporating  a  marginal 
repair  level  analysis  subroutine  into  the  existing  program. 
This  software  package,  acting  as  a  NRLA  preprocessor,  would 
calculate  ”  marginal  values  for  reparables,  and  then  deter¬ 
mine,  based  on  those  marginal  values,  whether  an  item  should 
be  scrapped,  depot  repaired,  intermediate  repaired  a(13:67> 
or  CIRF  repaired. 

The  final  aspect  of  NRLA  which  will  be  investigated  is 
its  ability  to  analyze  a  CIRF  designed  maintenance  plan.  The 
structure  of  the  NRLA  model  is  based  on  the  traditional 
three-level  maintenance  concept  (base,  intermediate,  depot> 
as  detailed  in  AFR  60-5.  Recently,  an  alternative  maintenance 
approach  known  as  Centralized  Intermediate  Repair  Facilties 
(CIRF)  has  emerged  as  an  efficient  way  to  effectively  support 
the  repair  of  a  weapons  system.  Normally  each  base  has  its 
own  dedicated  intermediate  repair  shop,  which  diagnose, 
repair,  and  replace  system  LRUs  and  SRUs.  On  the  other  hand, 
the  CIRF  concept  utilizes  a  centrally  located  intermediate 
repair  facility  to  service  multiple  bases  (Figure  2) .  In 
order  to  ensure  that  a  comprehensive  economic  analysis  is 
performed  when  determining  the  optimal  set  of  repair  level 
decisions  an  expansion  and  modification  of  the  program  is 
necessary  to  evaluate  the  cost  related  impacts  of  using  a 
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Figure  2.  Centralized  Intermediate  Repair  Facility  Concept 


CIRF  maintenance  approach.  This  will  be  the  second  of  the  two 
primary  objectives  of  this  thesis  effort,  the  development  of 
a  network  structure  capable  of  analyzing  a  CIRF  repair 
option.  If  possible  this  capability  will  be  incorporated 
into  the  existing  NRLA  software  program. 


Or oan i zat i on 

Chapter  II  encompasses  a  survey  of  the  current  litera¬ 
ture  available  relating  to  networking  algorithms  with  special 
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1 1  .  Literature  Review 


I n  troduc  t i on 

The  real-world  nature  of  the  thesis  subject  directed  the 
research  effort  into  two  principle  areas:  first  from  a  prac¬ 
tical  viewpoint,  a  study  and  understanding  of  current  USAF 
policies  and  procedures  regarding  Repair  Level  Analysis  (RLA) 
was  necessary  and  second,  from  a  theoretical  perspective  the 
investigation  of  alternative  state-of-the-art  optimization 
techniques  which  could  be  used  to  solve  the  multi-item  multi¬ 
echelon  repair  level  decision  problem  was  required.  It  was 
essential  that  each  of  these  subjects  be  addressed  to  ensure 
that  products  generated  from  this  effort  be  not  only  accepta¬ 
ble  in  terms  of  their  technical  content  and  structure,  but 
also  be  compatible  with  existing  current  USAF  logistics 
policy  and  doctrine. 

Repair  Level  Analysis  -  A  Macro  View 

In  an  effort  to  assist  Air  Force  acquisition  managers  in 
the  successful  design  and  implementation  of  efficient  mainte¬ 
nance  structures,  the  Department  of  the  Air  Force  has  estab¬ 
lished  specific  guidelines  and  procedures  to  be  followed  by 
both  the  government  and  contractors  in  the  design,  analysis, 
and  operation  of  the  RLA  program  (14:1).  As  stated  in 
AFSC/AFLC  Regulation  800-28,  Repair  Level  Analysis  (RLA) 


» 


the  RLA  program  has  two  major  objectives: 


<1>  Oes i on-or i en ted  RLA.  This  is  the  preliminary 
analysis  that  begins  in  the  conceptual  phase  of  the 
program  and  continues  through  the  critical  design 
review.  Its  goal  is  to  evolve  a  design  that 
considers  the  economics  of  support  alternatives  and 
produces  an  economical  life  cycle  cost  profile. 


<2>  Prov i si  on i no-or i en ted  RLA. This  analysis 
begins  during  the  full-scale  development  (FSD) 
phase  and  continues  into  the  production  program. 
The  objectives  of  this  phase  of  RLA  are  to  assign 
the  maintenance  portion  of  the  source,  maintenance, 
recoverabi 1 i ty  ( SMR)  codes  and  complete  maintenance 
planning  during  the  provisioning  process . * ( 1 4 : 1 ) 


Both  of  these  portions  of  the  RLA  program  are  conducted  using 
mathematical  models  and  techniques  to  determine  the  most 
appropriate  economic  level  of  repair  for  system  components. 
However,  as  can  be  seen  from  the  causal  diagram  in  Figure  3, 
there  are  a  multitude  of  factors  which  affect  the  operation 
of  a  logistics  structure.  This  is  not  only  reflected  in 
terms  of  its  costs,  but  also  in  terms  of  its  requirements 
for;  specialized  equipment,  manpower  with  certain  maintenance 
skills,  and  dedicated  facilities  at  the  repair  site. 
Unfortunately,  individual  identification  and  assessment  of 
these  factors  on  a  system's  operational  effectiveness  is  only 
an  initial  step  in  the  process  of  developing  a  total  systems 
perspective  of  the  repair  level  problem. 

Simulation  models  can  be  effective  tools  to  evaluate 
various  intergrated  logistics  management  approaches  because 
of  their  ability  to  treat  complex  interactions,  time-depen- 
dent  behavior,  and  system  feedback  mechanisms;  however,  they 
do  not  provide  'optimal'  solutions  but  rather  'acceptable' 
solutions.  Conversely,  purely  analytical  approaches,  while 
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proving  accurate  and  appropriate  -for  establishing  stockage 
requirements  during  the  initial  provisioning  process,  are 
subject  to  criticism  because  o-f  their  limiting  assumptions. 
The  most  damaging  of  these  is  that  of  applying  steady-state 
approaches  to  dynamic  systems.  Practically  speaking,  once 
the  system  begins  to  operate  structural  shifts  may  occur 
requiring  changes  to  the  repair  structure  (22:391-392). 

Aware  of  the  inherent  limitations  experienced  when  ap¬ 
plying  analytical  optimization  techniques  to  dynamic  systems 
such  as  the  multi-item  multiechelon  repair  system,  the  Air 
Force  has  initiated  several  efforts  to  assess,  on  a  compara¬ 
tive  basis,  the  life  cycle  costs  (LCC)  associated  with 
selection  of  a  repair  level  for  a  particular  item  within  a 
system.  Initial  efforts  were  centered  on  single  item  repair 
analysis;  however,  as  has  been  explained  earlier  in  Chapter 
I,  the  inability  to  adequately  address  assignment  of  the 
fixed  costs  related  to  such  items  as  support  equipment 
severely  damaged  the  credi tabi 1 i ty  and  validity  of  this 
techn i que . 

During  the  early  1970s,  the  Air  Force  began  to  investi¬ 
gate  the  feasibility  of  using  network  theory  as  a  basis  for 
solving  the  repair  level  problem.  Much  of  the  motivation  for 
this  effort  was  generated  by  an  article  by  J.M.U.  Rhys' 
entitled,  "A  Selection  Problem  of  Shared  Fixed  Costs  and 
Network  Flows'  (20:3).  Although  Rhys  does  not  specifically 
address  the  fact  that  the  repair  level  decision  problem  could 
be  structured  as  a  special  case  of  the  shared  fixed  cost 


Figure  4.  Shared-Fixed  Cost  Network 

problem,  he  is  generally  credited  with  making  the  major 
transition  -from  theory  to  application.  The  shared  fixed  cost 
problem  can  be  described  as  the  problem  of  selecting  a  set  of 
activities  from  some  larger  set  of  activities,  which  require 
the  allocation  of  fixed  overhead  costs  that  cannot  be 
specifically  related  to  one  particular  activity  (21:201). 

Using  the  duality  principle  of  linear  programming,  Rhys 
hypothesized  and  proved  that  the  shared  fixed  cost  problem 
could  be  solved  by  creating  a  directed  network  as  shown  in 
Figure  4.  By  determining  the  maximum  flow  on  this  network 
the  minimum  cost  decisions  are  simultaneously  identifying  the 
minimum  cut  set.  The  network  has  a  single  source,  S  and  a 
single  sink,  Z.  Sets  of  arcs  ESP]  and  CQZJ  are  defined  with 
each  arc  in  SP  corresponding  to  an  activity  (benefit),  P,  and 
each  arc  in  QZ  to  a  facility  (cost),  Q.  Each  arc  in  SP  is 
allocated  a  capacity  p,  where  p  is  the  net  gain  associated 
with  option  P,  and  each  arc  in  QZ  is  allocated  a  capacity  c, 


with  option  P,  and  each  arc  in  QZ  is  allocated  a  capacity  c, 
where  c  is  the  net  loss  associated  with  option  Q.  (21:201). 

The  MITRE  Corporation  of  Bedford,  Mass,  has  successfully 
employed  this  technique  on  a  variety  of  Air  Force  programs  to 
identify  optimum  levels  of  repair  for  weapons  systems  items 
(20).  Although  MITRE'S  efforts  were  commendable,  as  late  as 
1978,  a  standardized  software  package  implementing  these 
procedures  on  an  Air  Force  wide  level  was  still  missing. 
Being  the  office  of  primary  responsibility  (OPR)  in  the  USAF 
for  the  development  of  analitical  techniques  such  as  this, 
the  Concepts  and  Analysis  Division,  Air  Force  Acquisition 
Logistics  Center  (XRS/AFALC)  initiated  efforts  to  design, 
validate,  and  distribute  just  such  a  product.  By  June  of 
1980,  a  FORTRAN  IV  computer  program  known  as  the  Network 
Repair  Level  Analysis  (NRLA)  program  had  been  produced.  In 
conjunction  with  this  program,  two  documents,  The  NRLA  Users 
Gu i de  and  The  NRLA  Programmers  Guide,  were  published  to 
assist  in  the  understanding  of  the  structure,  limitations, 
and  capabilities  of  the  NRLA  program.  The  Users  Guide  pro¬ 
vides  a  general  description  of  the  model  design,  the  programs 
operation  in  terms  of  input  and  execution  requirements,  an 
explanation  of  the  LRU,  SRU,  and  SE  cost  computations,  as 
well  as,  a  presentation  of  the  network's  formulation.  The 
programmers  guide  is  designed  as  an  aid  to  understand  the  * 
programs  structure,  logic,  input  and  output  operations  and 
the  organization  of  data  so  that  modification  and/or  correc¬ 
tions  can  be  made  *  (2:1).  One  of  the  most  useful  portions  of 
the  programmers  guide  is  Appendix  A:  the  variables 
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dictionary,  which  is  an  alphabetical  listing  o-f  the  FORTRAN 
variable  names  and  arrays  used  throughout  the  NR LA  program. 

After  reviewing  the  existing  literature  relating  to  the 
development  and  application  of  networking  theory  to  repair 
level  analysis  in  the  Air  Force,  state-of-the-art  solution 
algorithms  were  researched.  These  different  techniques  were 
investigated  to  determine  the  potential  for  and  the 
appropriateness  of  their  application  within  a  specialized 
structure  such  as  the  RLA  network.  Initially  the  Max-fto 
Min-cut  problem,  which  is  the  theoretical  basis  for  the  NRLA 
formulation,  was  studied.  The  various  solution  procedures 
which  appeared  to  have  potential  are  discussed  in  the 
f ol 1 owi ng  pages. 


Max-flow  Min-cut  Problem 

The  max-flow  min-cut  theorem  which  was  initially 
established  in  the  late  50's  and  early  60's  by  Ford  and 


Fulkerson  is  the  basis  for  much  of  the  success  that  has  been 


experienced  in  formulating  an  allocation  problem  as  a 
directed  network.  Given  a  network  with  a  source  node  desig¬ 
nated  by  's'  and  a  sink  node  designated  by  't',  the  theorem 
is  as  follows:  aFor  any  network  the  maximal  flow  value  from 
s  to  t  is  equal  to  the  minimal  cut  capacity  of  all  cuts 
separating  s  and  t  *  (17:11).  The  following  example  should 
aid  in  understanding  of  this  theorem.  Consider  the  directed 
network  shown  in  Figure  5.  In  this  example  the  maximum  flow 
on  the  network  is  4  units,  with  1  unit  of  flow  on  path  s-1- 


Figure  5.  Example  Max-flow  Min-cut  Network 


s-2-l-t.  The  set  of  arcs  (s,l),  (2,1),  and  <2,t>  has 
capacity  4  and  is  the  minimum  cut  for  the  network.  This  set 
of  arcs  comprises  the  critical  flow  path  in  the  network;  to 
increase  flow  on  the  network  the  capacity  of  one  of  these 
arcs  must  be  increased.  Thus  the  capacity  of  the  minimum  cut 
is  equal  to  the  maximum  flow. 

Mathematical  Formulation 

The  relationship  between  the  max-flow  problem  and  the 
min-cut  problem  is  founded  in  the  duality  relationships  of 
the  two  problems.  The  mathematical  formulation  of  these 
problems  will  be  presented  here,  for  an  in  depth  treatment  of 
the  duality  relationships  refer  to  one  of  the  following 
works;  Ford  and  Fulkerson  (17:26-30),  Bazaraa  and  Jarvis 
(3:473-477),  or  Jenson  and  Barnes  (2:147-153). 


Probl em .  Consider  a  network  with  n  nodes  end  m 


A.  ",  VA  %  »  >  ’  . 


Max-f 1 ow 

arcs,  with  each  ere  (i,j>  haying  e  lower  bound  of  0  end  en 
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where  the  sums  end  inequalities  ere  taken  over  ell  of  the 
arcs  in  the  network. 

Before  proceeding  to  the  min-cut  problem,  the  concept  of 
a  cut-set  needs  to  be  explicitly  defined.  Let  N  be  the  set  of 
nodes  contained  in  a  network.  Let  Xi  be  any  subset  of  N  con¬ 
taining  node  1  (the  source)  but  not  node  n  (the  sink). 
Similarly  let  X  equal  N  -  X,  ,  the  subset  of  N  containing 
node  n  but  not  node  1.  Then  the  set  of  arcs  connecting  these 
two  subsets  is  called  the  cut-set.  For  other  than  trivially 
small  networks,  to  determine  the  number  of  possible  cut-sets 
becomes  a  problem  of  a  combinatorial  nature.  The  min-cut 
problem  seeks  to  identify  the  cut  set  with  the  minimum 
capac i ty . 

Min-cut  Problem.  The  dual  to  the  max-flow  problem  is 
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the  min-cut  problem.  Two  additional  variable  sets  are  added 

in  the  dual,  *  which  corresponds  to  the  conservation  equa- 

i 

tions  and  &  which  can  be  thought  of  as  "identifying  var- 

ij 

iables".  They  are  identifying  variables  because  in  the 

optimal  solution  if  &  >0  then  arc  <i,j>  becomes  a  member 

ij 

of  the  bottleneck  set  of  arcs  defining  the  cut-set.  It  is 
not  necessary  for  the  cut-set  to  be  a  unique  set  of  arcs. 
The  min-cut  problem  can  be  defined  as  follows: 

Minimize  Z\  X  "  c  & 
i*1  jal  ij  ij 

subject  to  It  -  *  =  1 

n  i 

<  -  It  +  S  >0 

i  j  ij 

&  >  0 

ij 

Labeling  Algorithms 

The  max-flow  problem  is  usually  solved  using  some  type  of 
labeling  algorithm.  There  are  a  multitude  of  these 
algorithms  available  to  solve  this  problem;  however,  many  of 
these  are  the  result  of  simply  changing  the  decision  rules 
for  selecting  the  flow  augmenting  paths.  The  classical 
approach  normally  used  to  solve  this  problem,  was  developed 
by  Ford  and  Fulkerson  and  uses  a  "Breadth  First  Search" 
algorithm  as  its  basis.  This  approach  has  been  fully 
developed  and  applied  to  the  max-flow  problem.  Additionally, 
several  refinements  to  this  algorithm  have  been  developed  and 
will  be  discussed.  Another  approach  to  this  problem  bases  its 
solution  procedure  on  a  "Depth  First  Search"  algorithm  as  its 
basis.  The  conceptual  basis  for  the  depth  first  approach  is 
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presented  in  Horowitz  and  Sahn i  (19:268-269).  However,  the 
literature  review  did  not  reveal  an  application  o-f  this 
search  algorithm  to  solving  the  max-flow  problem;  therefore 
an  algorithm  to  accomplish  this  was  developed  as  part  of  the 
thesis  effort.  A  brief  description  of  these  algorithms  will 
be  presented  with  an  in-depth  description  reserved  for 
Chapter  V. 

Before  proceeding  to  the  labeling  algorithms,  a  brief 
review  of  common  terminology  is  in  order.  In  both  algorithms 
labels  are  assigned  to  the  nodes  of  the  network.  Initially, 
all  nodes  except  the  source  node  are  unlabeled.  Uhen  a  label 
is  assigned  to  a  node  the  label  contains  two  pieces  of  infor¬ 
mation:  the  node  the  flow  came  from,  and  the  amount  of  flow 
potentially  available.  Consider  the  label  <l+,5>,  this  label 
indicates  that  5  units  of  flow  is  available  from  node  1.  The 
label  <2-, 3)  indicates  that  3  units  of  flow  which  had 
previously  been  flowed  to  node  2  has  the  potential  to  be 
returned  to  another  node  for  redirection.  It  should  be  noted 
that  the  labels  only  indicate  the  potential  amount  of  flow 
which  may  pass  between  nodes.  The  actual  amount  of  flow  and 
the  path  it  will  take  are  not  determined  until  a  path  to  the 
sink  node  is  found. 

The  node's  label  identifies  the  node  as  being  in  one  of 
three  states.  State  1  is  unlabeled,  state  2  is  labeled  but 
unscanned  and  state  3  is  labeled  and  scanned.  The  process  of 
scanning  occurs  when  a  node  has  all  of  the  arcs  originating 
or  terminating  at  it  checked  for  available  additional  flow. 


admissible  arcs. 


Acini  ssible  arcs  occur  in  two  forms:  forward  arcs  and  reverse 
arcs.  An  admissible  forward  arc  originates  at  the  node  being 
scanned*  has  a  destination  node  which  is  unlabeled,  and  has 
current  flow  less  than  capacity.  An  admissible  reverse  arc 
terminates  at  the  node  being  scanned,  has  a  source  node  which 
is  unlabeled,  and  has  current  flow  greater  than  the  lower 
bound  on  the  arcs  capacity. 

Breadth  First  Algorithm.  Initially,  all  nodes  are 
unlabeled  (state  1)  except  the  source  node  which  is  labeled 
and  unscanned  (state  2) .  The  source  is  then  scanned  and  all 
admissible  arcs  are  identified.  Their  associated  nodes  are 
labeled  and  their  state  changed  to  state  2.  After  each  node 
is  labeled,  a  check  is  made  to  determine  if  the  sink  node  has 
been  labeled.  If  it  has,  the  algorithm  goes  to  a  flow 
augmenting  procedure;  if  not,  it  continues.  When  the  source 
has  been  scanned,  a  labeled  but  unscanned  node  is  selected 
(usually  in  numerically  ascending  order)  and  is  scanned. 
This  continues  until  the  sink  node  is  labeled  or  there  are  no 
available  nodes  in  state  2.  If  there  are  no  nodes  in  state 
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2,  then  the  maximal  flow  has  been  identified  and  the 
algorithm  stops.  UJhen  the  sink  node  is  labeled,  the  flow 
augmenting  algorithm  comes  into  play.  It  identifies  the 
amount  of  flow  which  could  be  flowed  to  the  sink  from  its 
label.  It  then  uses  the  labels  to  trace  the  path  back 
through  the  network  to  the  source.  Each  arc  on  this  path  has 
its  flow  adjusted  by  the  amount  which  actually  reached  the 
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sink  node . 


Once  this  -flow  augmentation  is  complete  all 
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labels  are  erased  and  the  algorithm  restarted. 

Since  arc  capacities  are  restricted  to  being  only 
integer  values,  termination  of  this  algorithm  will  occur  in  a 
•finite  number  of  iterations.  Ford  and  Fulkerson  (17:121) 
provide  an  example  that  shows  that  the  algorithm  may  not 
terminate  or  terminates  with  a  sub-optimal  flew  when 
irrational  arc  capacities  are  allowed.  However,  a  problem 
can  still  arise  when  using  integral  capacities.  This 
involves  the  situation  where  the  total  number  of  flow 
augmentations  is  equal  to  the  maximum  flow  on  the  network. 
Edmonds  and  Karp  (16:250)  present  an  example  where  the  flow 
augmentation  at  each  iteration  is  only  one  unit,  this  leads 
to  an  upper  bound  being  created  for  the  number  of 
augmentations  equal  to  the  max-flow  of  the  network. 

Two  refinements  of  the  algorithm  are  presented  by  Edmonds 
and  Karp  (16:251-253)  which  can  reduce  the  upper  bound  on 
augmentations  required.  The  first  refinement  involves  deci¬ 
sion  rules  for  selecting  the  shortest  paths  for  flow 
augmentations.  They  show  that  if  each  flow  augmentation  is 
done  along  a  path  with  the  fewest  number  of  arcs  then  a  max- 
flow  will  be  identified  in  at  most  l/4(n  -n)  augmentations. 

This  result  is  derived  based  on  the  length  of  the  paths  to  a 
bottleneck  arc  from  the  source  and  the  sink,  and  the  max 
number  of  occurances  of  bottleneck  arcs  in  a  network  of  n 
nodes.  A  further  refinement  to  this  process,  includes  using  a 
rule  where  the  weights  of  the  paths  are  calculated  based  on 


their  lengths  and  these  weightings  are  then  used  to  aid  in 
the  selection  of  a  path.  Another  rule  involves  selecting  the 
path  with  the  -fewest  number  of  reverse  arcs. 


A  second  refinement  presented  uses  a  decision  rule  which 

selects  the  path  with  the  maximum  possible  number  of  flow 

augmentations  (16:253-255).  This  refinement  provides  that 

the  maximum  number  of  augmentations  necessary  to  attain  a 

maximal  flow  will  be  at  most  1+Log  f*(t,s),  where  f*(t,s> 

M/OH) 

denotes  the  value  of  maximal  flow.  Both  of  these 
refinements  are  intuitively  appealing  in  that  the  selection 
of  a  shortest  path  or  the  max-flow  per  augmentation  seem  to 
be  good  rules  of  thumb  for  obtaining  max-flow  in  the  fewest 
number  of  iterations.  Additional  considerations  for  speeding 
up  the  process  could  involve:  not  erasing  the  existing  set  of 
labels  after  each  iteration,  or  simply  saving  the  location  of 
the  last  node  with  excess  potential  on  the  flow  augmenting 
path.  This  could  save  a  great  deal  of  time  which  is  expended 

doing  redundant  labeling  and  should  significantly  improve 

computational  times. 

Depth  First  Algorithm.  As  in  the  'Breadth  First 
Algorithm”  the  initial  states  of  the  nodes  in  the  'Depth 
First  Algorithm*  are  source  node  in  state  2  and  all  other 
nodes  in  state  1.  The  search  for  a  flow  augmenting  path 
initiates  at  the  source  by  identifying  an  admissible  forward 

arc.  The  node  associated  with  this  arc  is  labeled  and  its 

state  is  changed  to  state  2.  A  search  is  then  made  from  this 
node  for  an  admissible  arc  and  the  process  is  repeated  until 
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the  sink  is  labeled  or  no  admissible  arcs  are  available  to 


augmenting  program  is  called  and  flow  is  augmented  along  the 
path.  As  the  flow  is  augmented  on  each  arc  a  check  is  made 
to  determine  if  there  is  any  potential  flow  remaining  at  any 
of  the  nodes  on  the  path.  If  there  is  no  flow  potential  at  a 
node,  that  node's  label  is  erased  and  it  is  returned  to  state 
1.  When  this  augmentation  is  complete,  the  algorithm  then 
starts  with  the  last  labeled  node  on  the  path  and  proceeds  to 
find  another  path  to  the  sink.  If  the  sink  cannot  be  reached 
from  the  last  node  on  the  path;  and  there  are  no  admissible 
arcs  at  this  node,  then  the  algorithm  backs  up  to  the 
previous  node  on  the  path  and  checks  for  admissible  arcs. 
This  process  continues  until  the  sink  is  labeled  or  until  the 
source  node  is  revisited,  at  which  point  a  new  admissible  arc 
is  chosen  from  the  source  and  the  algorithm  starts  again. 
When  there  are  no  further  admissible  arcs  at  the  source  and 
the  last  path  has  been  retraced  back  to  the  source,  the 
maximal  flow  has  been  found  and  the  algorithm  halts. 

To  ensure  that  the  maximal  flow  is  obtained  and  to 
identify  the  unique  set  of  labels  associated  with  this  flow, 
it  is  necessary  to  make  three  passes  through  the  network  with 
this  algorithm.  During  the  first  pass  only  forward  arcs  are 
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erased  at  the  end  of  the  first  and  second  passes  so  that  all 
possible  paths  can  be  explored  during  the  subsequent  passes. 
The  second  pass  is  used  to  identify  all  additional  paths 
through  the  network  on  which  flow  can  be  augmented.  These 
paths  may  contain  reverse  arcs.  The  third  pass  is  used  to 
identify  the  unique  set  of  labels  associated  with  the  minimum 
cut  set. 

Minimum  Cut  Aloorithm.  As  explained  earlier,  the 
optimum  set  of  repair  level  decisions  can  be  determined  by 
finding  the  minimum  cut  set  of  the  RLA  network.  Because  of 
this  relationship,  an  algorithm  which  could  directly  solve 
the  min-cut  problem  would  be  acceptable  for  use  with  RLA 
networks.  A  search  of  the  literature  indicated  that  the  most 
promising  algorithm  of  this  type  was  developed  by  Dessouky 
and  Phillips  (15)  and  is  called  the  "Cut  Search  Algorithm”. 
In  this  algorithm  the  min-cut  set  is  located  directly  by  a 
cut  seeking  procedure. 

The  min-cut  set,  which  is  designated  by  K*<s,t>,  is 
identified  by  a  two  stage  process  which  divides  the  set  of 
network  nodes,  designated  by  N,  into  three  groups  T,  W,  and 
S.  The  first  stage  of  the  process  starts  by  assigning  the 
sink  node  to  set  T.  Set  T  is  then  expanded  toward  the  source 
by  computing  the  value  of  its  cut  set  defined  by  K(T,T>  where 
N  =  T  +  T.  At  each  iteration  a  calculation  is  performed  for 
every  node  which  is  a  member  of  T  and  has  an  arc  connecting 
it  to  T  to  determine  if  it  can  be  added  to  T.  This  stage 
terminates  when  T  contains  only  the  source  node  or  if  no 
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nodes  in  T  can  be  added  to  T. 


If  T  contains  only  the  source 


node  then  K<T,T>  =  and  the  minimum  cut  set  has  been 
identified.  If  T  cannot  be  expanded  to  include  all  nodes 
except  the  source  then  the  second  stage  is  initiated.  In 
this  stage  the  sets  U  and  S  initially  contain  only  the  source 
node.  Set  W  is  expanded  toward  se t  T  until  U  =  T  at  which 
point  the  min-cut  set  is  contained  in  S.  S  is  contained  in 
Ul;  however,  it  is  expanded  only  by  adding  nodes  of  W,  which 
will  reduce  the  value  of  the  cut-set  as  defined  by  (S,S). 
Therefore  at  termination  K*<Sft>  =  K(S,S>. 

Computational  results  comparing  this  algorithm  with  the 
Out-of-Kilter  Algorithm  (OKA)  developed  by  Fulkerson  and 
Danzig  are  presented  in  Dessouky  and  Phillips"  article 
(15:403).  These  results  indicate  that  the  Cut  Search 
algorithm  is  more  efficient  than  the  OKA.  As  an  example;  for 
a  network  consisting  of  16  nodes  and  240  arcs  the  OKA  took 
5.33  seconds  of  cpu  time  while  the  cut  search  algorithm  took 
1.58  seconds  (15:403). 

Although  the  Cut  Search  Algorithm  appears  to  be  very 
efficient,  investigation  of  its  applicability  to  the  RLA 
problem  was  not  addressed  during  this  research  effort.  This 
algorithm  was  reserved  for  investigation  in  the  event  that 
the  depth  first  labeling  algorithm  proved  less  efficient  than 
anticipated  and  to  allow  for  the  investigation  of  the  appli¬ 
cability  of  generalized  network  primal  algorithms  described 
in  the  next  section. 
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Generalized  Network  Primal  Algorithms.  In  recent  years 
significant  advances  have  been  made  in  the  development  and 
implementation  of  a  class  of  algorithms  which  can  solve 
large  scale  capacitated  transportation  and  transshipment 
problems.  These  algorithms,  often  referred  to  as  primal 
network  codes,  can  also  be  applied  to  the  solution  of  the 
max-flow  problem.  Because  of  this  characteristic  and  the 
potential  for  applying  NRLA  to  large  size  RLA  systems  such  as 
the  Bl-Bomber  and  MX-Missile  programs,  investigation  of  one 
of  these  primal  network  codes  seemed  appropriate.  The  code 
chosen  for  use  in  this  context  was  a  commercially  available 
code  called  GNET.  This  code  was  selected  for  two  reasons; 
(1)  it  was  currently  available  on  the  computer  systems 
accessable  to  AFIT  and  (2)  it  was  the  subject  of  an  excellent 
article  in  Management  Science  by  Bradley,  Brown  and  Graves 
(7)  which  covered  in  detail  its  design,  implementation  and 
use.  This  type  of  in-depth  coverage  of  a  commercially 
marketed  algorithm  is  usually  very  difficult  to  find  or  is 
non-ex i stent . 

Due  to  the  complexity  of  the  GNET  code,  and  the  fact  that 
it  was  applied  as  a  subroutine  in  NRLA  to  solve  the  max-flow 
problem  only  a  brief  description  of  the  GNET  code  will  be 
given.  For  an  in-depth  explaination  of  this  particular 
application  of  primal  networking  one  should  refer  to  the 
above  referenced  article.  Basically,  the  GNET  code  solves 
the  general  linear  programming  <LP>  problems  (transportation, 
transshipment,  and  max-flow)  by  specializing  these  LP  prob- 
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network  is  the  bounded  variable  revised  simplex  method.  The 
uniqueness  of  this  approach  however,  lies  in  the  fact  that 
the  network  which  is  developed  by  the  algorithm  takes  advan¬ 
tage  of  the  sparsity  of  the  arc  node  matrix  and  uses  upper 
tr i angu 1 ar i zed  basis  representations  to  quickly  solve  the 
network . 

Results  published  by  Glover  and  Klingman  indicate  that 
when  this  code  is  used  for  finding  solutions  to  transporta¬ 
tion  and  transshipment  problems,  that  it  is  30%  to  40%  faster 
than  the  OKA  code  in  the  solution  of  transportation  and 
transshipment  problems.  Bradley,  Brown  and  Graves  indicate 
that  it  is  currently  believed  that  primal  implementations  are 
faster  and  require  less  storage  than  OKA  or  other  algorithms 
used  to  solve  these  type  problems  (7:3). 
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III.  Me thodol oc 


introduct i on 


Initially,  a  study  was  made  of  the  logic  and  rationale 
for  each  part  of  the  NRLA  program.  The  NRLA  program, 
originally  developed  on  the  AFLC  Honeywell  635  computer  using 
the  CREATE  time  sharing  system  is  "composed  of  a  main 
routine,  a  block  data  subroutine,  plus  twelve  additional 
subroutines"  (2:1).  The  logic  diagram  in  Figure  6  displays 
the  interrelationships  of  the  various  subroutines  to  the  main 
program.  Appendix  B  provides  a  brief  explai nation  of  their 
respective  functions.  After  several  discussions  with  both 
the  developers  and  users  of  the  model  (6,9),  it  was  deter¬ 
mined  that  wherever  possible  a  structured  programming 
approach  using  FORTRAN  subroutines  would  be  used  to 
incorporate  new  or  enhanced  features  into  the  model.  By  using 
this  type  of  approach,  there  would  be  minimal  impact  on 
current  users  of  the  model.  Additionally,  it  was  recognized 
that  a  strategy  of  this  type  would  probably  facilitate  the 
verification  and  validation  phases  of  the  design  process. 
With  this  general  concept  in  mind,  the  identification  of  the 
necessary  tasks  related  to  accomplishing  the  thesis  objec¬ 
tives  was  initiated.  In  conjunction  with  this  effort,  a  flow 
diagram  was  developed  to  order  and  prioritize  these  tasks. 
This  was  done  to  better  visualize  how  they  the  total  process. 
This  conceptual  process  resulted  in  the  flow  diagram  shown  in 
Figure  7.  To  successfully  achieve  the  thesis  objectives  of: 
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reducing  the  NRLA  program's  execution  time,  and  integrating 
the  CIRF  repair  level  option,  -four  new  subroutines  were 
designed  and  evaluated  -for  possible  integration  into  the 
existing  NRLA  structure.  These  are  called  MARLA,  MSETNT , 
NMXFLO,  and  CRFSET  and  are  indicated  by  the  dashed  lines  in 
Figure  8. 

The  MARLA  subroutine  employs  the  concept  o-f  marginal 
value  analysis  in  an  attempt  to  determine  i  -f  optimal  deci¬ 
sions  can  be  made  -for  individual  items.  As  can  be  seen  -from 
Figure  8,  the  MARLA  subroutine  functions  as  a  preprocessor  or 
filter  to  the  MSETNT  subroutine.  In  effect,  when  MARLA  makes 
an  optimal  decision  for  an  LRU  or  SRU  failure  mode,  it  has 
reduced  the  total  number  of  items  which  need  to  be  further 
analyzed.  It  was  believed  that  gains  in  time  efficiency 
would  be  achievable  from  this  enhancement;  the  network  could 
be  built  faster  in  MSETNT,  and  the  problem  to  be  solved  in 
MAXFLO  would  be  of  smaller  magnitude.  This  assumes  that  the 
execution  time  requirements  for  the  MARLA,  MAXFLO,  and  MSETNT 
subroutines  would  be  less  than  what  is  currently  required  for 
the  SETNET  and  MAXFLO  subroutines  in  NRLA. 

The  second  area  which  exhibited  a  potential  opportunity 
for  reducing  NRLA's  execution  time  was  in  the  MAXFLO  subrou¬ 
tine.  MAXFLO's  purpose  is  to  determine  the  minimum  cost  set 
of  optimal  repair  level  decisions.  The  solution  technique 


currently  employed  is  a  direct  application  of  the  two  stage 
labeling  procedure  developed  by  Ford  &  Fulkerson  (17:22>. 
Implementation  of  a  new  labeling  procedure  appeared  likely  to 
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yield  significant  time  savings  in  this  area  o-f  the  program. 
The  alternative  labeling  technique  which  was  developed  is 
constructed  in  NMXFLO  and,  as  will  be  explained  later,  can  be 
used  as  an  alternative  to,  or  in  conjunction  with,  the  exist¬ 
ing  MAXFLO  subroutine. 

The  final  objective  of  incorporating  a  CIRF  maintenance 
concept  as  a  fourth  repair  alternative  was  proposed  to  be 
accomplished  by  using  a  third  subroutine  called  CRFSET.  Sub¬ 
routine  CRFSET  was  used  with  the  existing  SETNET  subroutine 
to  determine  if  the  CIRF  repair  option  could  be  the  optimal 
decision  for  a  particular  item.  It  compares  the  repair  level 
recommendation  generated  by  SETNET,  (depot,  base,  or  scrap) 
with  the  CIRF  option  by  building  a  new  unique  network  (See 
Figure  8). 

Although  integration  of  each  of  these  subroutines  would 
necessitate  modifications  to  other  portions  of  the  NRLA 
program,  variable  and  array  integrity  was  to  be  maintained  to 
the  maximum  extent  possible.  Additionally,  the  individual 
subroutines  were  developed  with  the  objective  of  being 
compatible  with  existing  program  data  structures,  as  well  as, 
with  linkages  between  data  elements  (i.e.  the  use  of 
pointers).  Separately,  these  subroutines  accomplish  specific 
subobjectives  which  contribute  to  the  successful  achievement 
of  the  previously  outlined  thesis  objectives.  The  following 
three  chapters  will  discuss  in  detail  the  concept  develop¬ 
ment,  computerization,  verification,  validation,  and 
analysis  phases  of  each  of  these  three  major  enhancements. 


IV .  Marginal  Analysis  Enhancement 

Concept  Development.  The  application  of  marginal  cost 
analysis  techniques  to  solve  the  repair  level  decision  prob¬ 
lem  is  conceptually  based  upon  the  idea  that  the  analyst  can 
use  existing  cost  in-formation  that  is  readily  available  to 
quickly  identify  the  optimal  level  o-f  repair  for  a  specific 
item.  AFSC/AFLC  Pamphlet  800-4,  Repair  Level  Analysis  Proce¬ 
dures.  identifies  marginal  analysis  as  a  simple  and  efficient 
way  to  synthesize  the  repair  level  problem  into  a  more 
workable  format  (13:32-33).  An  appealing  feature  of  using 
marginal  analysis  for  repair  level  decision  making  is  its 
ability  to  make  decisions  without  prorating  the  cost  of  the 
support  equipment  which  would  be  necessary  to  repair  an  item. 
This  precludes  the  possibility  of  making  inconsistent  or 
suboptimal  decisions.  It  is  for  these  reasons  that  develop¬ 
ment  of  a  marginal  analysis  approach  was  selected  as  a 
possible  solution  to  solving  the  problem  associated  with 
NRLA's  computer  run  time. 

To  repair  or  replace  an  item  or  component  of  an  Air  Force 
system  requires  the  expenditure  of  materiel  and  personnel 
resources,  as  well  as  the  utilization  of  specialized  and 
sophisticated  diagnostic  test  equipment.  The  costs  which  are 
generated  from  this  process  can  be  broadly  categorized  into 
two  types  of  costs:  (1)  Pipeline  costs  <P>  -  which  are 
defined  as  those  costs  incurred  directly  to  repair  the  item, 
and  (2)  Support  Equipment  costs  (SE)  -  which  are  defined  as 
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Table  II 


Repair  Level  Costs 


I  .  Pipeline  Costs 

A .  Initial  Spares 

B.  Packing  and  Shipping 

C.  Replenishment  Spares 

D.  Repair  Labor 

E.  Item  Management 

F.  Technical  Management 

G.  Training  of  Maintenance  Personnel 


II.  Support  Equipment  Costs 

A.  Acqu i si t i on 

B.  Operation  and  Maintenance 

C .  Fac i 1 i t i es 


III.  Scrap  Costs 

A.  Initial  Spares 

B.  Replenishment  Spares 


those  costs  related  to  SE  use  that  cannot  be  assigned  to  a 
specific  item.  Table  II  gives  an  itemized  listing  of  each  of 


these  categories  of  cost.  Once  these  costs  have  been  deter¬ 
mined  for  each  possible  level  of  repair,  a  marginal  analysis 
can  be  made  of  the  SE  resources  at  each  repair  location.  In 
the  context  of  comparing  the  four  repair  options  of  Depot, 
Base,  CIRF,  or  Scrap,  if  the  scrap  cost  of  the  item  is  known 
to  be  less  than  the  pipeline  costs  of  the  item  at  the  other 
locations,  then  the  scrap  decision  will  always  be  the  optimal 
decision  for  that  item.  In  this  situation,  the  scrap  deci¬ 
sion  is  said  to  dominate  the  other  three  alternatives 


(13:32) . 


Equations  (1)  through  (3)  represent  the  mathemati- 


cal  interpretation  of  the  marginal  analysis  comparisons. 
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pipeline  costs  to  scrap  the  item, 
pipeline  costs  to  depot  repair  the  item, 
pipeline  costs  to  Cl RF  repair  the  item, 
pipeline  costs  to  base  repair  the  item, 
cost  o-f  the  depot  support  equipment, 
cost  o-f  the  CIRF  support  equipment, 
cost  o-f  the  base  support  equipment. 


In  order  to  determine  the  marginal  value  o-f  the  SE 
resources  for  each  level,  the  pipeline  costs  associated  with 
repair  of  the  item  are  subtracted  from  the  pipeline  costs  to 
scrap  the  items 


MV  =  P  -  P 
SD  S  D 

^sc  =  Ps  ~  Pc 

^SB  PS  ”  PB 


<4> 

<5> 

(<6) 


From  the  above  computations,  it  can  be  seen  that  in  this 
initial  example  that  the  marginal  value  of  the  SE  at  every 
level  was  negative.  This  would  mean  that  there  would  not  be 
any  cost  savings  to  be  realized  from  buying  the  necessary  SE 
to  repair  the  item.  In  fact,  before  the  purchase  of  SE 
resources  could  be  economically  justified,  either  the 
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pipeline  scrap  cost  would  have  to  increase  or  the  pipeline 
costs  of  one  of  the  repair  options  would  have  to  decrease 
until  one  of  the  SE  marginal  values  became  zero. 
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Unfortunately,  many  situations  do  not  lena  themselves  to 
such  a  straight  forward  type  of  analysis.  A  much  more  com¬ 
prehensive  set  of  decision  criteria  must  be  developed  when 
dealing  with  SE  that  possess  a  positive  marginal  value.  Based 
on  the  scenario  of  having  three  possible  repair  locations,  as 
well  as,  the  option  to  discard  the  item,  a  comparison  must 
now  be  made  not  only  between  the  pipeline  costs  of  the  item; 
but  also  between  the  total  costs  to  repair  the  item  for 
each  of  the  alternative  levels.  Table  III  displays  in  a 
systematic  fashion  the  necessary  cost  comparisons  which  must 
be  made  to  ascertain  whether  or  not  a  repair  option  can  be 
eliminated.  If  any  of  the  statements  are  true  for  a  given 
repair  level,  that  particular  repair  option  is  eliminated. 
To  interpet  the  meaning  of  these  relationships  with  respect 
to  the  marginal  value  of  the  SE,  a  joint  marginal  value  is 


calculated  as 
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The  joint  marginal  value  is  always  taken  as  the  minimum 
of  the  three  values  because  it  identifies  the  maximum  amount 
of  funds  that  should  be  allocated  for  the  SE  at  that  repair 


location.  As  soon  as  the  cost  of  the  SE  exceeds  this  amount, 
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Marginal  Analysis  Decision  Criteria 
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the  depot  repair  option, 
the  CIRF  repair  option, 
the  base  repair  option, 
total  repair  costs  at  depot, 

total  repair  costs  at  CIRF. 

total  repair  costs  at  base. 


it  would  no  longer  be  cost  effective  to  purchase  the  support 
equipment  to  repair  the  item  (13:71).  A  deductive  type  of 
process  is  used  to  isolate  the  optimal  decision.  Only  if 
three  o-f  the  -four  choices  can  be  eliminated  does  the  marginal 
analysis  approach  identi-fy  the  optimal  decision.  However,  if 
marginal  analysis  were  to  be  used  as  a  separate  analysis 
program  the  utility  of  its  approach  shhould  not  be  based 
solely  on  its  ability  to  identify  optimal  decisions.  But 
rather  a  more  accurate  measure  would  be  the  total  number  of 
repair  options  it  could  exclude  from  the  decision  making 
process . 

Computer i zat i on/1 n tear at i on .  To  implement  the  concept  of 
marginal  repair  level  analysis  a  FORTRAN  IV  program  called 
MARLA  was  developed.  Appendix  D  contains  a  computerized 
listing  of  the  MARLA  subroutine.  This  subroutine  is  designed 
to  be  called  by  the  MAIN  subroutine  subsequent  to  the 
computation  of  the  LRU/SRU  pipeline  repair  costs  which  are 
performed  in  the  LRUCMP,  FMCMP,  and  SRUCMP  subroutines,  but 
pr  ror  to  the  computation  of  the  SE  requirements  which  are 
made  in  subroutine  SECMP.  The  reason  that  the  MARLA 
subroutine  should  be  called  in  this  manner  is  two  fold: 
first,  to  accurately  perform  the  marginal  analysis  MARLA 
needs  the  pipeline  costs  for  each  of  the  items,  and  second, 
any  SE  that  is  purchased  as  a  result  of  a  MARLA  optimal 
decision  should  be  considered  when  determining  the  additional 
SE  requirements  for  the  remaining  items. 


Figure  9.  Example  of  Proposed  Arc  Elimination  Process 


The  most  difficult  interface  that  had  to  be  designed 
for  integration  of  the  MARLA  subroutine  was  the  modification 
of  subroutine  SETNET  to  accept  the  MARLA  decisions.  When 
optimal  decisions  were  made  by  MARLA,  it  eliminated  the 
requirement  to  construct  the  arcs  which  would  be  necessary  in 
MAXFLO  to  analyze  the  item.  The  example  in  Figure  ? 
indicates  the  arcs  that  could  be  erased  from  a  network  if  an 
optimal  decision  is  made;  in  this  case  if  an  optimal  decision 
had  been  made  for  LRU  #1,  eight  out  of  fifteen  of  the  network 
arcs  (#1  , #4, #6, #7, #9, #10  ,#1!  are  no  longer  needed. 
Using  this  idea,  modifications  were  made  to  subroutine  SETNET 


to  give  it  the  capability  to  reconfigure  a  new  network  based 


strictly  on  just  those  items  still  requiring  a  repair  level 
decision.  This  new  subroutine  was  called  MSETNT .  By 
implementing  this  change,  greater  efficiencies  were  expected 
to  be  realized  during  the  programs  operation.  The  code  which 
was  developed  to  accomplish  this  task  is  contained  in 
Appendix  E. 

Because  item  unique  information  such  as  unit  cost  and 
SE  utilization  hour  requirements  are  stored  in  separate  ar¬ 
rays  for  LRU  and  SRU  failure  modes,  MARLA  was  divided  into 
two  major  sections  of  code;  one  to  analyze  LRU  items  and 
another  to  handle  SRU  items.  Proper  control  of  the  NRLA 
exclusion  arrays  throughout  both  of  these  portions  of  the 
program  was  determined  to  be  extremely  important  when  using 
MARLA  as  a  preprocesser  to  NRLA.  There  were  three  primary 
reasons  for  this.  First,  if  a  repair  level  had  been  excluded 
by  the  user,  it  is  not  necessary  to  perform  a  marginal  analy¬ 
sis  for  that  particular  combination  of  item  and  repair  level. 
Second,  when  the  actual  network  is  constructed  in  subroutine 
SETNET,  any  repair  levels  that  have  been  excluded  have  their 
respective  arc  costs  set  to  a  very  large  cost  called  JUMBO. 
This  prevents  them  from  ever  entering  the  cut  set  of  optimal 
decisions.  However,  the  MARLA  decision  criteria  were  not 
sufficient  to  exclude  a  repair  option  from  further  considera¬ 
tion  by  NRLA.  Consequently,  the  FORTRAN  code  which  was 
written  to  accomplish  the  marginal  analysis  was  reequired  to 
embody  this  situation.  The  third  effect  MARLA's  integration 
had  on  the  exclusion  arrays  occured  in  subroutine  OUTPUT. 


Here  the  exclusion  arrays  are  checked  again  to  see  it  the 


user  has  eliminated  any  o-f  the  repair  options.  For  those 
items  where  repair  levels  have  been  omitted,  the  program  will 
print  the  word  'XCLD' .  Two  programming  changes  were  required 
in  subroutine  OUTPUT  to  ensure  that  this  operation  was 
performed  for  only  user  exclusions  and  not  for  exclusions 
generated  by  MARLA. 

Another  important  aspect  of  subroutine  MARLA  is  the  iden¬ 
tification  and  acquisition  of  the  correct  quantities  of  SE 
resources  necessary  to  repair  the  item.  Although  when 
performing  the  marginal  analysis  the  potential  SE  costs  are 
calculated  for  each  item,  the  actual  procurement  of  these 
resources  should  not  be  made  unless  an  optimal  decision  has 
been  made  for  that  item.  For  this  reason  the  logic  of  the 
MARLA  subroutine  does  not  adjust  the  SE  arrays  for  cost  and 
hours  until  the  optimal  decision  arrays,  OPDECL  and  OPDECS, 
are  screened  to  determine  if  in  fact  a  valid  requirement  does 
exists  for  that  particular  piece  of  support  equipment.  MARLA 
is  structured  to  buy  the  minimum  number  of  types  of  SE  needed 
for  each  repair  level  based  on  the  total  system  requirement. 
This  was  accomplished  by  using  the  logic  described  in  Figure 
10.  The  pivotal  factor  in  this  procedure  was  the  accurate 
determination  of  the  SE  hours  currently  available  for  each 
type  of  SE  at  each  repair  level.  This  precluded  the  purchase 
of  more  SE  than  was  actually  needed.  The  array,  HRAVSE,  was 
created  and  integrated  into  the  common  MARLA  block  statement 
to  store  the  available  hours  for  each  piece  of  SE  and  to 


«*  «* 


allow  -for  these  hours  to  be  transferee!  to  the  subroutine, 


SECMP.  This  ensured  that  any  SE  hours  still  available  could 
be  applied  against  the  utilization  hours  required  for  the 
remaining  items  to  be  analyzed. 

L>er i f i cat i on/Val i dat i on .  To  verify  the  operational 
accuracy  of  the  MARLA  enhancement,  two  undertakings  were 
necessary.  Not  only  was  it  essential  to  confirm  the  logic  of 
the  MARLA  subroutine,  but  it  was  also  critical  to  separately 
authenticate  the  major  changes  that  had  to  be  made  to  the 
SETNET  subroutine.  It  was  decided  that  these  tests  should  be 
conducted  independently  to  provide  the  maximum  assurance  that 
both  portions  of  the  program  were  performing  correctly. 

Initial  testing  of  the  MARLA  subroutine  revealed 
significant  program  logic  flaws.  By  modifying  the  subroutine 
so  that  it  could  be  run  as  a  separate  program  errors  were 
quickly  identified  and  corrected.  This  would  not  have  been 
possible  if  MARLA  haa  v«Jf'  immediately  integrated  into  the 
NRLA  program.  The  sample  data  set  which  was  used  for  this 
testing  included  a  total  of  18  LRU/SRU  failure  inodes  with 
three  different  types  of  SE  at  each  of  the  repair  locations. 
Prior  to  running  the  program,  the  correct  repair  decisions 
were  manually  calculated  using  the  decision  criteria  in  Table 
III  on  page  46.  The  results  from  the  MARLA  program  were  then 
compared  to  these  recommendations.  This  iterative  process 
continued  until  the  MARLA  decisions  matched  identically  with 
the  true  decisions.  The  results  obtained  from  this  test 


indicated 


that  the  MARLA  subroutine 


was  f unc  t i onal 1 y 


correctly  both  computationally  and  logically. 


A  different  design  was  developed  to  verify  the 
accuracy  of  the  SETNET  modifications.  A  predefined  network 
was  utilized  to  determine  if  MSETNT  would  eliminate  the  the 
appropriate  arcs  for  a  given  LRU/SRU  optimal  decision.  By 
presetting  the  OPDECL  and  OPDECS  arrays  with  a  particular 
optimal  decision,  the  arcs  to  be  eliminated  from  the  network 
could  easily  be  identified.  A  trace  of  the  network, 
consisting  of  each  arc's  source  node  and  destination  node, 
was  output  to  determine  if  the  correct  arcs  had  been  removed. 
Initially,  the  program  eliminated  only  a  portion  of  the  arcs 
it  should  have.  However,  by  continuing  to  use  the  trace  from 
the  network,  all  of  the  logic  and  syntax  coding  errors  were 
satisfactorily  resolved. 

During  this  process,  an  important  point  was 
discovered.  If  the  optimal  decision  for  an  LRU  is  'base", 
the  arcs  representing  the  transportation  costs  for  replace¬ 
ment  SRUs  <446,449,4410  in  Figure  8,  pg.48)  could  also  be 
eliminated.  However,  since  the  costs  associated  with  these 
arcs  are  related  to  the  repair  of  the  SRU,  it  was  realized 
that  they  should  be  added  to  the  SRU  depot  and  scrap  repair 
arcs  <443  and  418,  Figure  8)  so  that  the  SRU  costs  would  not  be 
underest imated. 

Initial  validation  of  both  the  MARLA  and  MSETNT  sub¬ 
routines  was  accomplished  by  consulting  with  the  original 
developers  of  the  NRLA  program  <6>.  However,  complete  valida¬ 
tion  of  the  MARLA  enhancement  will  only  be  determined  after 


NRLA  users  throughout  the  Air  Force  and  defense  contracting 
community  have  had  the  opportunity  to  exercise  this  option 
and  provide  -feedback  as  to  its  performance  and  accuracy. 

Anal  ys  i  s/Resu  Its.  The  -feasibility  of  using  the  MARLA 
subroutine  as  a  preprocesser  to  the  NRLA  program  was  analyzed 
in  terms  of  its  effect  upon  total  processing  time  and  addi¬ 
tional  storage  requ i remen ts .  A  data  set  which  had  previously 
been  analyzed  by  the  original  NRLA  program  was  used  with  the 
enhanced  program.  In  this  way,  a  standard  of  comparison 
could  be  established  to  measure  the  effect  on  program  time 
and  storage  requirements.  Additionally,  the  repair  level 
recommendations  produced  could  be  compared  for  consistency. 

When  used  as  external  preprocessor,  the  MARLA  program 
produced  results  identical  to  the  NRLA  program.  However, 
when  integrated  as  a  subroutine  in  the  main  program,  the 
repair  recommendations  generated  did  not  agree.  It  was 
observed  that  although  the  pipeline  logistical  costs  for  each 
repair  level  matched  exactly  for  both  programs,  the  SE  cost 
estimates  differed  significantly.  It  was  determined  by  -cross 
checking  the  outputs  from  both  programs  that  this  was 
problably  the  reason  for  the  disparity  in  the  model's 
recommendat ions. 

Two  other  significant  results  were  observed  from  com¬ 
paring  the  output  of  the  two  programs.  These  related  to  the 
processing  time  and  the  number  of  arcs  constructed  in  each  of 
the  networks.  These  program  parameters  are  directly  propor¬ 
tional  to  each  other;  the  greater  the  number  of  arcs  in  the 


network,  the  longer  it  take  the  program  to  execute.  For  the 
sample  data  set,  which  consisted  of  23  LRU  -failure  modes,  10 
SRUs,  and  12  pieces  o-f  SE ,  the  original  NRLA  network 
consisted  o-f  361  arcs  and  required  .0149  hours  to  execute. 
For  the  enhanced  program  the  output  indicated  that  the 
network  contained  only  12  arcs  and  processed  in  .0034  hours. 
Although  these  results  showed  a  substantial  improvement  over 
the  original  network,  a  high  level  o-f  con-fidence  could  not  be 
placed  in  their  validity  because  of  the  programs  inability  to 
produce  decisions  consistent  with  the  original  NRLA  program. 

Storage  requirements  for  the  new  program  were  also 
analyzed.  To  incorporate  the  MARLA  subroutine  nine  new 
arrays  were  required.  Seven  of  these  were  needed  to  perform 
the  marginal  analysis  while  the  two  remaining  ones,  HRAVSE 
and  TSECST,  were  used  to  collect  the  SE  cost  and  available 
hours,  respectively.  The  number  of  arrays  was  a  function  of 
the  number  of  repair  locations.  The  size  of  the  arrays  used 
in  the  marginal  analysis  was  determined  by  the  total  number 
of  failure  modes  and  SRUs  to  be  analyzed.  Appendix  A 
contains  a  listing  of  each  of  these  arrays  and  their 
function  in  the  program. 

Based  on  the  limited  test  results  which  were  achieved 
with  MARLA  during  this  study,  further  efforts  are  necessary 
in  order  to  fully  assess  its  capability  as  a  preprocessor  to 
the  NRLA  networking  process.  A  better  evaluation  could  then 
be  made  of  its  utility  in  terms  of  reducing  the  overall 
program  processing  time  versus  the  additional  storage 


requirements  to  achieve  this. 


55 


.  Max— flow  Min-cut  Labeling  Enhancement 

NRLA  Network  Formulation 

Network  Structure.  Understanding  the  operation  of  the 
NRLA  program  requires  an  in-depth  knowledge  o-f  its  network 
structure  and  how  this  structure  models  the  cost  -factors 
associated  with  the  repair  level  decisions.  The  ten  decision 
cost  -factors  and  the  eleven  logistics  cost  -factors  which 
comprise  them  are  shown  in  Table  1  on  page  8.  As  can  be 
seen  in  this  table,  the  decision  cost  -factors  are  simply  the 
summation  o-f  the  logistic  costs  associated  with  a  LRU's  or 
SRU's  repair  level  decision.  This  consolidation  o-f  costs 
allows  their  use  as  capacities  -for  the  arcs  o-f  the  repair 
level  network.  The  network  -formulation  was  examined  briefly 
in  Chapter  I;  however,  the  details  of  how  the  network  is 
built  and  how  the  various  cut-sets  (repair  level  decisions) 
are  determined  were  not  covered  at  that  time.  The  following 
section  covers  these  topics  in  detail. 

The  NRLA  program  builds  the  network  structure  by  first 
assigning  numbers  to  the  nodes  in  a  specific  order  and  then 
laying  in  the  arcs  associated  with  these  nodes.  As  shown  in 
Figure  1  on  page  7,  the  nodes  are  assigned  numbers  in  the 
following  order  source  node,  depot  support  equipment,  depot 
LRU  failure  modes,  depot  SRUs,  base  LRU  failure  modes,  base 
SRUs,  base  support  equipment,  and  finally  the  sink  node. 
This  careful  ordering  of  the  nodes  allows  the  program  to 
develop  forward  and  backward  pointers  to  the  vector  arrays 
containing  the  LRU,  SRU,  and  SE  relationships. 
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Next  the  arcs  are  added  to  the  network.  This  is  also 
done  in  an  ordered  manner  starting  with  the  source  node  and 
working  through  the  nodes  one  at  a  time.  For  simplicity  the 
LRU  -failure  mode  nodes  will  be  re-fered  to  as  LRU  nodes  during 
this  discussion.  The  arcs  -from  the  source  to  the  depot  SE 
(depot  SE  cost)  are  added  -first,  -followed  by  the  arcs  to  the 
depot  LRUs  (depot  LRU  repair  costs)  and  the  SRUs  (depot  SRU 
repair  costs).  The  next  arcs  added  are  the  depot  SE  to  LRU 
and  SRU  JUMBO  (large  capacity)  arcs.  These  arcs  shown  as  the 
thicker  arcs  in  the  -figure  represent  the  relationships 
between  the  individual  depot  SE  and  the  set  o-f  LRU/SRUs  which 
it  repairs.  Because  these  arcs  have  JUMBO  capacity  they 
cannot  be  capacitated  (-flow  equals  capacity);  therefore,  they 
allow  the  total  set  o-f  LRU/SRUs  to  share  the  burden  o-f  the  SE 
cost . 

The  arcs  originating  at  the  depot  LRUs  are  added  next. 
First  the  depot  LRU  to  depot  SRU  (DEC1)  arc  is  added.  Then 
the  depot  LRU  to  base  LRU  arc  (scrap  cost).  The  DEC1  arcs 
represent  a  portion  o-f  the  additional  costs  associated  with 
the  decisions  to  base  repair  the  LRU  and  depot  repair  the 
SRU.  The  rest  o-f  these  costs  are  contained  on  the  DEC2  arcs. 
The  costs  are  split  to  allow  -for  representation  o-f  the  two 
way  trip  associated  with  SRU  depot  repair  verses  the  one  way 
trip  if  the  SRU  is  scrapped  when  the  LRU  is  base  repaired. 

The  SRU  scrap  cost  arcs  are  added  next.  They  originate 
at  the  depot  SRU  and  terminate  at  the  base  SRU.  Moving  on  to 
the  base  LRU  nodes  the  following  arcs  are  added.  First  an 


arc  -from  the  base  LRU  to  its  associated  base  SRU  (DEC2  costs) 
and  then  an  arc  -from  the  base  LRU  to  the  sink  node  (base  LRU 
repair  costs).  Following  these  the  arcs  emanating  -from  the 
base  SRUs  are  added.  The  -first  arc  is  a  JUMBO  arc  -from  the 
base  SRU  to  its  associated  base  LRU.  This  arc  ensures  that 
the  decisions  to  base  repair  the  SRU  and  scrap  or  depot 
repair  the  LRU  cannot  be  made.  Next  an  arc  -from  the  base  SRU 
to  the  sink  node  is  added  (SRU  base  repair  costs). 

The  last  nodes  dealt  with  are  the  base  SE  nodes.  They 
have  several  arcs  associated  with  them.  The  -first  arcs  added 
are  the  JUMBO  arcs  -from  the  set  o-f  base  LRUs  to  the  SE  which 
are  used  in  their  repair.  Next  the  JUMBO  arcs  -from  the  set 
of  base  SRUs  to  the  SE  are  added,  and  finally  the  base  SE 
cost  arc  is  added  from  the  SE  to  the  sink.  Each  base  SE  is 
considered  separately  and  all  the  above  arcs  added  prior  to 
going  to  the  next  SE.  The  JUMBO  arcs  perform  the  same 
function  for  the  base  SE  as  they  did  for  the  depot  SE,  that 
of  inter-relating  the  LRU/SRU  sets  to  support  the  SE  costs. 

Minimum  Cuts.  Once  the  network  is  built,  the  set  of 
repair  level  decisions  are  determined  by  applying  a  max-flow 
min-cut  algorithm  to  identify  the  min-cut  set.  Figure  11 
displays  the  seven  potential  RLA  network  cuts.  Table  IU 
identifies  the  decisions  associated  with  these  cuts  and  the 
costs  incurred  for  each  decision.  Recall  that  a  cut  divides 
the  network  into  two  sets  of  nodes:  one  containing  the  source 
node  and  one  containing  the  sink.  The  arcs  of  the  cut-set  are 
the  chokepoints  of  the  network  flow.  Cut-set  1  in  Figure  11 


TABLE  IV 

RLA  Decision  Costs 


RLA  DECISION  COSTS 


Decision 


Included  Costs 


LRU 

SRU 

SE 

DEP 

BASE 

DEC  1 

DEC  2 

Depot 

Depot 

Yes 

No 

No 

No 

Base 

Base 

No 

Yes 

No 

No 

Scrap 

Scrap 

No 

No 

No 

No 

Depot 

Scrap 

Yes 

No 

No 

No 

Base 

Scrap 

No 

Yes 

No 

Yes 

Base 

Depot 

Yes 

Yes 

Yes 

Yes 

Scrap 

Depot 

Yes 

No 

Yes 

No 

can  be  seen  to  contain  the  -following  arcs:  (1,2).  (1,3),  and 
(1,4).  These  represent  the  costs  associated  with  the  depot 
SE,  the  depot  LRU  repair  costs,  and  the  depot  SRU  repair 
costs.  This  indicates  that  when  the  maximum  flow  was  deter¬ 
mined  in  the  network,  these  arcs  were  capacitated;  therefore, 
the  optimal  repair  decisions  are  depot  repair  of  both  the  LRU 
and  SRU. 

The  following  notes  on  the  indicated  cuts  are  necessary 
to  ensure  that  the  cut-sets  are  not  misinterpreted.  JUMBO 
arcs  cannot  be  capacitated;  therefore,  all  cuts  crossing 
JUMBO  arcs  do  not  include  these  arcs.  Cut-set  4  does  not 
include  the  DEC1  arc  because  it  originates  at  an  unlabeled 
node  and  terminates  at  a  labeled  one;  to  be  part  of  the  cut- 


set  the  reverse  must  be  true 


Cut  7  represents  the  anomolous 


situation  where  the  LRU  is  scrapped  but  the  SRU  is  depot 
repaired.  The  cost  of  the  SRU  is  included  in  the  LRU  cost, 
therefore  a  LRU  scrap  decision  should  include  scrapping  the 
SRU.  The  structure  of  the  network  cannot  exclude  the 
possibility  of  this  cut  occuring.  If  it  does  occur,  a  valid 
decision  can  be  reached  by  running  the  model  twice:  the  first 
time  excluding  the  LRU  scrap  option,  and  the  second  time 
forcing  the  LRU  and  SRU  scrap  decisions.  The  optimal 
repair  level  decisions  are  associated  with  the  minimum  of 
these  costs. 

Now  that  the  RLA  network  is  built  and  the  potential  cut¬ 
sets  explained,  it  is  necessary  to  examine  the  efficiency  of 
the  algorithm  used  to  solve  the  network  max-flow  min-cut 
problem.  The  next  section  covers  the  labeling  algorithms 
applied  to  this  problem. 

Max-flow  Min-cut  Algorithmic  Improvements 

Conceptual  Development.  A  primary  concern  expressed  by 
the  users  of  the  NRLA  model  is  the  amount  of  computer 
resources  used  if  the  sensitivity  analysis  option  is  selected 
(6,9).  The  NRLA  program  allows  the  user  to  perform 
sensitivity  analysis  on  a  selected  range  of  costs  and/or 
failure  rates,  based  on  mean  time  between  failures  (MTBF) , 
for  each  LRU,  LRU  failure  mode,  and  SRU.  This  allows  the 
user  to  accomplish  a  "what  if"  type  analysis  to  determine  the 
best  repair  level  options  in  the  event  that  one  of  these 
variables  changes  during  the  systems  development.  To 
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program 


accomplish  this  sensitivity  analysis,  the  NRLA 
modifies  each  item's  cost  or  MTBF  and  then  re-sol ves  the 
network  based  on  this  modification.  It  then  resets  the  net¬ 
work  and  performs  the  next  modification.  To  get  the  true 
impact  of  this,  consider  a  sample  system  which  consists  of  23 
LRU  failure  modes  and  10  SRUs  with  both  cost  and  MTBF  sensi¬ 
tivity  analysis  selected.  In  this  example  the  sensitivity 
analysis  program  would  modify  and  re-solve  the  network  over 
300  times.  Each  time  the  network  is  modified,  subroutine 
MAXFLO  is  called  to  re-solve  it.  Therefore,  the  computational 
efficiency  of  the  algorithm  used  in  MAXFLO  to  solve  the  max- 
flow  min-cut  problem  plays  a  major  role  in  determining  the 
total  computer  run  times  for  the  model. 

Based  on  this  information,  the  next  step  in  the  research 
effort  involved  analyzing  the  current  labeling  algorithm  for 
potential  time  saving  modifications.  To  assist  in  this 
analysis  a  detailed  description  of  the  Breadth  First  Labeling 
algorithm  is  presented  next  <1 7> . 

Breadth  First  Labeling  Algorithm 

Notation  :  f..  =  flow  from  node  i  to  node  j 

|J 

c..  =  capacity  from  node  i  to  node  j 

w 

Concept  :  Assign  labels  to  the  nodes  of  the  network  such  that 
a  flow  path  may  be  determined  on  the  network.  Augment 
flow  on  these  paths  until  the  maximum  flow  has  been 
attained.  The  labels  have  the  form  (i+,e)  or  <i-,e>, 
where  i  is  a  node  of  the  network  and  e  is  a  positive 


integer  or  •  according  to  the  rules  specified  in  proce¬ 
dure  A.  The  state  o-f  a  node  is  modi-fied  by  the 
algorithm  according  to  whether  it  is  unlabeled  (state 
1),  labeled  and  unscanned  (state  2),  or  labeled  and 
scanned  (state  3). 

Procedure  A  :  Start  by  assigning  the  source  ,  s  ,  a  label  o-f 

(-,«>,  this  indicates  that  it  is  receiving  infinite  flow 

from  an  external  source.  All  other  nodes  are  unlabeled 

(state  1)  and  the  state  of  the  source  node  is  state  2. 

In  general,  select  any  node  i  which  is  in  state  2. 

Suppose  it  is  labeled  <k±  ,  e(i>>.  Scan  node  i  for 

admissible  arcs.  Recall  that  admissible  forward  arcs 

require  the  destination  node  j  to  be  unlabeled  and 

<  c  ,  while  admissible  reverse  arcs  require  the 
'•>  ij 

source  node  j  to  be  unlabeled  and  f  >0.  For  all 

admissible  forward  arcs  assign  the  appropriate  node  j 

the  label  (i+,e(j>>  where  e(j)  =  min  (  e(i>,  c.  -f  ). 

'J  ij 

For  all  admissible  reverse  arcs  assign  the  appropriate 

node  j  the  label  (i-,e(j)>,  where  e(j>  =  Min  (e(i),f..>. 

J  I 

As  each  node  is  labeled  a  check  is  made  to  determine  if 

the  sink  has  been  labeled.  If  it  has  proceed  to 

Procedure  B,  if  not  continue  the  labeling  process.  If 
node  i  is  completely  scanned  and  the  sink  is  unlabeled 
change  node  i's  state  to  state  3  (labeled  and  scanned). 
Repeat  the  general  step  until  either  the  sink  node  t  is 
labeled,  or  until  no  additional  labels  can  be  assigned 
and  the  sink  is  unlabeled.  If  the  sink  is  labeled  go  to 


Procedure  B,  i  f  not,  stop  -  the  maximum  -flow  has  been 
-found. 

Procedure  B  :  The  sink  node  t  has  been  labeled  <j+,e<t>>. 
Trace  the  labeled  path  back  to  the  source  node  s  while 
augmenting  the  -flow  on  the  arcs  o-f  this  path.  In 
general,  i f  node  j  is  labeled  <i+,e(x))y  replace  f 

with  e<t>  +  -f  ,  and  i -f  j  is  labeled  <i-,e<x>>,  replace 

ij 

f  wi  th  -f  e<t).  Proceed  to  node  i  and  continue 

ji  ji 

this  -flow  augmentation  until  the  source  is  reached, 

erase  all  the  current  labels  and  return  to  Procedure  A. 

Analysis  o-f  the  above  labeling  algorithm  indicates  that 
it  uses  the  labels  it  produces  very  inefficiently.  The 
following  simple  example  will  more  clearly  illustrate  this 
problem.  Consider  the  network  shown  in  Figure  12.  The 
source  node  is  connected  to  100  nodes  which  are  in  turn 

connected  to  the  sink.  At  each  iteration  of  the  above 
labeling  algorithm  all  100  arcs  from  the  source  to  the  inter¬ 
mediate  nodes  will  be  checked  for  admissibility  and  the  nodes 
labeled  if  appropriate.  After  the  source  is  completely 
scanned  the  algorithm  proceed  to  the  next  labeled  but 
unscanned  node  and  then  scans  it.  In  this  example  one  of  two 
things  can  happen  when  an  intermediate  node  is  scanned; 

either  the  sink  is  labeled  or  the  node  is  marked  as  scanned 
and  a  new,  labeled  but  unscanned,  node  is  sought.  When  the 
sink  is  labeled,  flow  is  augmented  on  the  path  from  the 
source  to  the  sink,  all  the  labels  are  erased,  and  the 

process  is  restarted.  If  the  algorithm  has  to  scan  n  inter- 


mediate  nodes  before  it  -finds  a  path  to  the  sink,  then  there 
can  be  as  many  as  100-n  nodes  which  are  labeled  and  then 


erased  without  being  used.  Now  suppose  that  all  arcs  of  the 
network  initially  have  flow  less  than  capacity.  Given  this 
initial  condition,  for  this  example  the  algorithm  would 
produce  99  excess  labels  on  the  first  iteration,  98  on  the 
second  and  so  forth  ,  for  a  total  of  4950  excess  labels. 
Clearly  this  is  a  less  than  optimal  use  of  the  labels 
produced . 

One  way  to  eliminate  this  problem  involves  using  a  Depth 
first  search  algorithm  in  the  labeling  process.  A  review  of 
the  literature  indicated  that  the  application  of  depth  first 
search  concepts  to  the  labeling  process  has  not  been 
attempted  before.  Conceptually  what  a  depth  first  search 
does,  is  move  through  the  network  from  the  source  to  the 
first  level  and  then  directly  on  to  the  second  level  and  so 
forth  until  the  sink  is  reached.  A  level  in  this  context 
consists  of  the  set  of  nodes  whose  shortest  path  from  the 
source  is  the  same  length,  where  length  is  given  in  total 
number  of  arcs  on  the  path.  In  the  above  example  this  would 
be  equivalent  to  labeling  node  2  and  then  immediatly  labeling 
the  sink  from  node  2,  thus  eliminating  the  labeling  of  nodes 
3  through  101.  RLA  networks  are  well  suited  to  this  type 
search,  in  that  only  two  levels  separate  the  source  and  sink 
in  most  cases.  The  above  stated  advantages  of  a  depth  first 
search  concept  led  to  its  use  in  developing  a  new  labeling 
algorithm  to  solve  the  max-flow  problem.  A  detailed 
description  of  this  algorithm  follows. 


Depth  First  Labeling  Algorithm 


Notation  :  -f  =  -flow  -from  node  i  to  node  j 

c  =  capacity  from  node  i  to  node  j 
ij 

Concept  :  Assign  labels  to  the  nodes  of  the  network  such  that 
a  flow  path  may  be  determined  on  the  network.  Augment 
flow  on  these  paths  until  the  maximum  flow  has  been 
attained.  Labels  are  assigned  to  the  nodes  on  the  path 
from  the  source  to  the  sink  only.  Nodes  are  scanned 
only  until  an  admissible  arc  is  found,  at  which  point 
that  node  is  labeled  and  scanned  for  an  admissible  arc. 
The  state  of  a  node  is  modified  according  to  whether  it 
is  unlabeled  state  1  or  labeled  state  2.  Nodes  are  not 
labeled  as  scanned  or  unscanned. 

Procedure  A  :  Initially  all  nodes  are  unlabeled  (state  1>. 
Proceed  by  labeling  the  source  node  s,  (-,•>,  and 
setting  it  to  state  2. 

First  Pass  (forward  Paths  Only) 

Starting  at  the  source  check  for  an  admissible  forward 

arc  and  label  its  associated  node  i  (s+,c  -f  >,  assign 

*'  si 

node  i  to  state  2  and  scan  it  for  an  admissible  forward 
arc.  In  general  proceed  in  this  fashion  until  the  sink 
is  labeled  or  the  last  node  on  the  path  has  no 
admissible  forward  arcs.  If  the  sink  is  labeled  go  to 
Procedure  B.  If  there  are  no  admissible  forward  arcs 
back  up  to  the  previous  node  on  the  path  and  scan  it  for 
an  admissible  forward  arc,  if  there  is  one  proceed  with 
the  general  rule,  if  not  continue  backing  up.  When  no 


admissible  forward  arcs  are  available  at  the  source 
erase  all  labels  and  go  to  the  second  pass. 

Second  Pass  (Forward  and  Reverse  Arcs  allowed) 

Once  again  label  the  source  (-,«►)  and  set  it  to  state 
2.  Scan  the  source  for  an  admissible  arc  and  label  its 
associated  node.  Scan  this  node  for  an  admissible  arc 
and  label  its  associated  node  according  to  whether  it  is 
a  forward  or  reverse  admissible  arc.  In  general  continue 
in  this  fashion  until  the  sink  is  reached  or  no 
admissible  arcs  are  available  at  the  last  node  on  the 
path.  If  the  sink  is  labeled  go  to  Procedure  B.  If  no 
achnissible  arc  is  available  back  up  to  the  previous  node 
on  the  path  and  proceed.  Uhen  no  admissible  arcs  are 
available  at  the  source  erase  all  labels  and  proceed  to 


the  Third  Pass.  At  this  point  the  maximum  flow  has  been 
i dent i f i ed. 

Third  Pass  (labeling  pass) 

Repeat  the  procedure  used  for  the  second  pass  until  no 
admissible  arcs  are  available  at  the  source.  Uhen  this 
occurs  stop  -  the  max-flow  with  correct  labels  to 
identify  the  min-cut  has  been  found. 

Procedure  B  :  The  sink,  node  t,  has  been  labeled  (j+,e(t)>. 
Trace  the  labeled  path  back  to  the  source,  node  s,  while 
augmenting  the  flow  on  the  arcs  of  this  path.  In 
general,  if  node  j  is  labeled  (i+,e(x>>,  replace  f 
with  e(t>  +  f  ,  and  if  j  is  labeled  (i-,e(x)>,  replace 

u 

f..  with  f,  -  e(t).  As  the  flow  is  augmented  on  this 


path,  each  node  is  checked  to  see  i -f  it  has  potential 


-flow  remaining.  If  e(j>  >  e(t)  and  no  previous  node  on 
the  path  has  had  potential  remaining,  then  this  node  is 
stored  as  the  new  starting  point  when  the  flow 
augmentation  is  complete.  If  e(j>  <  e(t)  then  change 
the  nodes  state  to  1  and  erase  its  label.  Nodes  with 
potential  remaining,  stay  in  state  2  and  retain  their 
labels.  When  the  flow  augmentation  is  complete  return 
to  the  appropriate  pass  in  procedure  A  and  continue  the 
pass  by  scanning  either  the  stored  node  or  the  source. 

Several  factors  contribute  to  the  overall  efficiency  of 
the  above  labeling  algorithm.  The  algorithm  pushes  the  flow 
to  the  sink  on  the  first  available  path,  thereby  eliminating 
any  unused  labels.  It  retains  the  position  of  the  deepest 
node  on  the  with  flow  potential  remaining.  This  allows  for  a 
head  start  when  starting  the  search  for  the  next  flow 
augmentation  path.  It  also  stores  an  array  of  the  arcs  on 
the  path  from  the  source  to  the  sink.  This  eliminates  the 
need  to  identify  these  arcs  by  individually  scanning  each 
node  on  the  path.  The  breadth  first  search  must  identify  the 
path  in  this  way  because  the  entire  path  is  not  identified 
until  the  sink  is  labeled. 

Computer i zat i on/Inteorat i on .  The  coding  of  the  depth 
first  labeling  algorithm  was  accomplished  using  FORTRAN  IV. 
This  version  of  FORTRAN  was  used  to  ensure  compatibility  with 
the  existing  program  and  to  allow  for  transpor tab i 1 i ty 


between  computer  systems.  The  fully  developed  subroutine  is 
called  NMXFLO  and  is  listed  in  Appendix  F.  The  only  part 
missing  from  this  listing  are  the  Common  statements  from  the 
NRLA  program.  NMXFLO  was  designed  to  replace  MAXFLO,  and  as 
such  can  be  used  in  place  of  or  in  addition  to  MAXFLO  in  the 
existing  NRLA  model.  The  addition  of  one  array  (NAPATH)  and 
two  variables  were  the  only  additions  to  NRLAs  storage 

requ i rements. 

Ver i f i cat i on/Ual i dat i on ■  The  verification  of  this 
subroutine  was  accomplished  by  developing  a  network  data  set 
with  known  arc  capacities  and  flows.  This  data  set  was  then 
input  to  NMXFLO  and  the  resulting  arc  flows  compared  with  the 
expected  results.  In  this  manner  initial  logic  flaws  in  the 
algorithm  were  identified  and  corrected.  When  the  algorithm 
was  fully  debugged  and  operating  as  anticipated  it  was 

applied  to  a  validation  phase  of  testing  by  incorporating  it 
in  the  NRLA  model  and  using  it  to  solve  the  max-flow  min-cut 
problem.  The  results  of  this  testing  were  compared  with 
solutions  using  MAXFLO  (the  original  subroutine)  for 
validity.  In  all  instances  NMXFLO  produced  identical 
resul ts. 

Anal ys i s/Resu Its.  The  final  analysis  performed  on  the 
NMXFLO  subroutine  dealt  with  determining  its  efficiency  at 
solving  the  max-flow  min-cut  problem.  To  accomplish  this, 
the  NRLA  model  was  run  using  data  sets  supplied  by  the  B-l 

SPO  and  AFALD/XRS.  The  execution  times  for  the  NMXFLO  and 

MAXFLO  subroutines  were  calculated  using  the  CREATE  computer 
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Processing  Time  in 

seconds 

Col  d) 

(new) 

Data 

Set 

* 

LRU 

M  M  M  M 

SRU  SE  NODES  ARCS 

Breadth 

First 

Depth 

First 

GNET 

without  sensitivity  analysis 

1 

23 

10  12  100  361 

0.432 

0.038 

3.041 

2 

52 

0  106  212  470 

8.454 

0.082 

6.054 

with 

sensitivity  analysis 

1 

23 

10  12  100  361 

37.764 

7.977 

** 

2 

52 

0  106  212  470 

495.939* 

15.419 

** 

* 

run  t 

me  exhausted 

** 

GNET  was  not  compared  because  of  the  way  i t  was 
implemented  in  the  model. 

system  processing  time  program.  This  program  clocks  the 
actual  processing  time  in  hours.  The  time  hacks  were  taken 
at  the  beginning  o f  the  subroutines  and  just  before  the 
return  to  the  main  program.  The  resulting  times  were  then 
multiplied  by  3600  to  convert  the  times  to  computer 
processing  seconds.  The  results  of  these  tests  are  shown  in 
Table  V.  As  can  be  seen,  the  depth  first  search  algorithm  is 
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s  i  gn  i  f  i  can  1 1  y  -faster  than  the  the  breadth  -first  search 
algorithm.  The  largest  system  tested  showed  the  depth  -first 
algorithm  with  100  fold  increase  in  efficiency  over  the 
breadth  first  algorithm.  The  results  obtained  when  GNET  was 
used  to  solve  the  max-flow  problem  in  NRLA  are  also  contained 

in  Table  V.  In  this  case,  GNET  did  not  overcome  the  ! 

efficiency  of  the  breadth  first  algorithm  until  the  system 
size  increased  to  the  maximum  size  tested,  and  it  never 

approached  the  efficiency  exhibited  by  the  depth  first 

algorithm.  GNET  was  not  tested  for  the  sensitivity  option 
due  to  the  manner  in  which  is  was  implemented  in  the  NRLA 
program.  The  modular  form  in  which  GNET  was  implemented  did 
not  allow  it  to  retain  information  on  the  network  structure. 

This  would  result  in  it  taking  the  total  network  solution 
time  at  each  step  of  the  sensitivity  analysis.  The  other 
algorithms  took  advantage  of  previous  information  and, 
therefore,  any  comparison  between  them  in  this  area  would  be 
inaccurate.  Based  on  these  results  it  is  strongly 
recommended  that  the  depth  first  search  algorithm  be 
incorporated  as  the  primary  MAXFLO  subroutine  in  the  current 


NRLA  model 


*JI  ,  Centralized  Intermediate  Repair  Facility  Enhancement 


r9z 


Concept  Development 

To  address  the  feasibility  of  including  a  CIRF 
maintenance  option  in  the  NRLA  program  required  an  in-depth 
study  of  the  existing  network  formulation.  The  rigid 
structure  of  the  existing  network  was  not  adaptive  to  the 
insertion  of  a  fourth  repair  option.  However,  it  was 
hypothesized  that  by  constructing  a  second  network,  which 
would  represent  the  optimal  decision  from  the  first  network 
and  the  CIRF  alternative,  that  a  comparative  economic 
analysis  could  be  made. 

To  successfully  implement  this  concept  into  the  program, 
five  major  tasks  had  to  be  accomplished.  First,  a  definitive 
description  of  the  meaning  of  CIRF  repair  was  required. 
Without  this,  it  would  be  very  difficult  to  validate  the 
computerization  of  the  CIRF  repair  option.  The  second  major 
requirement  was  the  identification  of  the  necessary  input 
data  elements  to  support  a  CIRF  developed  network.  The  third 
objective  related  to  the  development  of  CIRF  cost  equations 
to  calculate  pipeline  as  well  as  support  equipment  costs. 
Fourthly,  the  various  network  structures  that  could  be 
generated  from  the  decisions  of  the  first  network  had  to  be 
constructed,  coded,  and  analyzed.  Finally,  the  impact  from 
the  integration  of  each  of  these  modifications  had  to  be 
assessed. 

The  CIRF  repair  option  is  based  on  a  generalized 
i n terpe tat i on  of  how  the  CIRF  maintenance  concept  is 
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currently  employed  by  the  majority  o-f  the  system  program 


offices.  Extension  of  the  existing  NRLA  assumptions  were 
made  with  one  exception.  Presently,  the  NRLA  model  assumes 
an  uniform  distribution  for  the  basing  of  the  total  number  of 
weapons  systems  being  analyzed.  However,  in  the  CIRF 
scenario,  there  are  additional  deployment  factors  to  be 
considered.  Although  the  number  of  systems  at  each  base  is 
the  same,  the  overriding  question  is  the  number  of  CIRF 
locations  and  the  number  of  bases  that  are  being  serviced  by 
that  repair  location.  As  an  illustration,  in  the  common 
situation  where  there  are  two  centralized  servicing  facili¬ 
ties,  one  in  the  CONUS  and  one  overseas,  the  number  of  bases 
linked  to  each  CIRF  can  create  an  imbalanced  workload.  In  an 
attempt  to  compensate  for  this  imbalance,  it  was  decided  that 
an  existing  systems  variable  called  "OS'',  the  fraction  of  the 
force  overseas,  would  be  applied  to  SE  utilization  require¬ 
ments.  The  application  of  this  variable  limits  the  number  of 
CIRF  locations  to  two;  one  overseas  and  one  in  the  CONUS.  It 
is  applied  by  taking  the  standard  CIRF  factors  and 
multiplying  them  by  the  overseas  fraction.  These  numbers  are 
then  used  to  calculate  the  SE  requirements  for  the  overseas 
CIRF  location.  The  same  logic  is  applied  to  the  CONUS  CIRF 
SE  requirements,  only  now  the  multiplication  factor  is 
(l-OS).  This  results  in  a  more  accurate  assessment  of  the 
true  SE  requirements  for  the  CIRF  alternative. 

Input  Data  Requirements 

To  generate  a  CIRF  repair  option  a  list  of  the 


necessary  CIRF  variables  was  compiled.  Once  this  had  been 
done,  the  appropriate  READ  and  WRITE  statements  were  located 
in  the  program  and  modified  to  accept  these  additional 
variables.  A  list  of  these  new  variables  and  their  meaning 
is  contained  in  Appendix  A.  One  particular  situation  that 
occurred  during  this  process  related  to  the  changes  that  had 
to  be  made  to  accomodate  the  CIRF  support  equipment. 
Currently,  the  NRLA  program  uses  the  first  digit  of  the  value 
stored  in  the  SECODE  array  to  identify  both  the  type  of  SE 
and  the  level  at  which  it  is  used.  The  digits  and  their 
meaning  are  as  follows: 


Common  SE  (Depot) 
Peculiar  SE  (Depot) 
Add"!  SE  (Depot) 

SE  Software  (Depot) 
Common  SE  (Base) 
Peculiar  SE  (Base) 
Add' l  SE  (Base) 

SE  Software  (Base) 


Since  the  SECODE  array  is  formatted  for  a  four  digit  integer 
value  this  allowed  a  user  to  input  up  to  999  different  types 
of  SE  in  any  of  these  classifications.  To  differentiate 
between  the  CIRF  SE  and  the  depot  or  base  SE,  and  to  also 


maintain  the  capability  to  have  a  maximum  of  999  pieces  of 
SE,  the  SECODE  values  were  expanded  to  five  digit  numbers. 


The  first  two  digits  of  the  SECODE  value  were  used  for  the 
CIRF  SE  as  follows: 


Common  SE  (CIRF) 
Pecu 1 i ar  SE  (CIRF) 
Add' 1  SE  (CIRF) 

SE  Software  (CIRF) 
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Thus  a  piece  of  support  equipment  that  was  previously 
identified  with  the  number  2013  will  now  be  identified  as 


V  % 
%',V 


t.. 


W 


V- 


20013.  If  this  piece  of  equipment  needed  to  be  used  at  the 
CIRF  location  its  identification  number  would  be  92013. 

Cost  Requirements 

Each  of  the  NRLA  cost  subroutines  < FMCMP ,  LRUCMP,  and 
SRUCMP)  were  expanded  to  incorporate  a  capability  to  compute 
CIRF  related  pipeline  and  SE  costs.  The  CIRF  equations  were 
developed  to  be  compatible  with  the  existing  depot,  base,  and 
scrap  cost  equations.  Using  the  CIRF  input  variables,  the 
CIRF  pipeline  costs  were  calculated  based  on  the  eleven 
logistical  costs  in  Table  I  on  page  8.  Two  additional  cost 
arrays,  CIRFLC  and  CIRFSC,  were  needed  to  store  these  costs 
for  LRUs  and  SRUs,  respectively. 

To  compute  the  CIRF  SE  costs,  the  SE  requirements  for 
each  of  the  CIRF  locations  were  separately  determined.  These 
requirements  were  then  added  together  to  obtain  the  total 
needs  of  the  system.  This  value  was  then  multiplied  times 
the  unit  cost  of  the  SE  < SE  array  CADB)  to  obtain  the  total 
CIRF  acquisition  costs.  The  standard  cost  factor  for  opera¬ 
tions  and  maintenance  contained  in  SE  array  CODB  was  added  to 
the  acquisition  cost  along  with  any  costs  that  were  known  to 
be  incurred  for  special  facilities  < SE  array  FDB) .  To  obtain 
a  per  base  cost  the  summation  of  these  costs  was  then  divided 
by  the  total  number  of  bases  in  the  system.  These  changes 
are  contained  in  subroutine  SECMP  in  Appendix  H.  Having 


calculated  the  costs  for  a  CIRF  repair  option,  it  was  now 


necessary  to  design  the  correct  RLA  networks  for 
option. 


CIRF  Network  Conceptual  Development 

Initially,  the  incorporation  of  CIRF  into  the  NRLA  model 
was  visualized  as  requiring  only  a  simple  modification  of  the 
network  structure  to  include  the  CIRF  option.  However,  this 
proved  to  be  false.  A  major  restriction  on  tnis 
modification  was  that  it  had  to  be  accomplished  without 
loosing  the  ability  to  identify  the  optimal  repair  level 
decisions  through  the  min-cut  set.  Upon  investigation,  it 
became  apparent  that  the  inclusion  of  CIRF  in  this  manner 
could  not  be  accomplished.  There  is  no  network  structure 
that  retains  the  ability  to  identify  a  unique  repair  level 
decision  from  the  four  repair  options.  This  is  because  of 
the  unique  use  of  a  maximal  flow  algorithm  to  solve  the 
minimal  cost  problem. 

The  next  option  considered,  involved  exploring  the 
feasibility  of  making  the  decisions  in  a  two-pass  scenario. 
The  initial  pass  would  compare  the  standard  depot-scrap-base 
configuration  and  the  second  pass  would  then  compare  these 
optimal  decisions  against  the  CIRF  option.  This  method  of 
incorporating  CIRF  proved  to  be  feasible  and  was  developed 
for  inclusion  in  NRLA. 

An  initial  problem  to  be  overcome  with  this  method  was 
how  to  handle  the  seven  potential  SRU/LRU  decision  pairs 
which  can  result  from  the  first  pass.  In  each  instance  the 
network  generating  subroutine  must  identify  the  SRU/LRU 


decision  pair  and  build  an  appropriate  network  structure  so 
that  the  resulting  comparison  with  CIRF  will  provide  valid 


repair  level  decisions.  The  potential  for  seven  differently 
constructed  subnetworks  brings  up  the  problem  of  how  to 
identify  the  optimal  decisions  in  these  different  structures. 
The  solution  to  these  problems  led  to  the  development  of  the 
set  of  seven  subnetworks  shown  in  Figures  13(a)  thru  13(g). 
As  can  be  seen  in  these  figures,  the  basic  RLA  network 
structure  is  retained  for  each  subnetwork.  One  major 
difference  is  that  both  the  depot  and  base  SE  are  now 
included  on  the  right-hand  side  of  the  network  structures. 

CIRF  Network  Construction.  Constructing  the  networks 
displayed  in  Figure  13  required  a  major  modification  of  the 
SETNET  subroutine.  The  manner  in  which  the  nodes  were 
assigned  changed  very  little.  The  source  is  assigned  first, 
followed  by  the  CIRF  SE,  then  the  CIRF  LRUs  and  SRUs.  These 
nodes  replace  the  depot  equivalents  in  the  original  RLA 
networks.  The  next  nodes  assigned  are  the  LRU  and  SRU  nodes 
associated  with  the  optimal  decisions  from  the  first  NRLA 
pass.  These  nodes  take  the  place  of  the  base  LRU  and  SRU 
nodes  in  the  original  RLA  networks.  Next  the  depot  SE  nodes 
are  assigned,  followed  by  the  base  SE  nodes,  and  then  the 
sink  node.  This  structure  requires  that  the  total  number  of 
nodes  be  increased  by  the  number  of  CIRF  SE. 

The  manner  in  which  the  arcs  are  connected  to  these  nodes 
is  the  key  ingredient  in  obtaining  the  optimal  repair  level 
decisions.  The  arcs  originating  at  the  source  and  CIRF  SE 


are  the  same  as  those  in  the  original  RLA  network.  The  -first 
arcs  requiring  special  attention  are  those  that  originate  at 
the  Cl RF  LRUs. 

The  arcs  originating  -from  the  CIRF  LRUs  are  assigned  in 
the  -following  order.  First  the  arc  connecting  the  LRU  to  the 
associated  CIRF  SRU,  and  then  the  arc  connecting  the  LRU  to 
the  NRLA  decision  LRU.  The  first  arcs  capacity  will  be  DEC  3 
if  the  LRU  decision  from  the  first  pass  is  base  repair,  and 
JUMBO  otherwise.  The  DEC  3  costs  are  those  pipeline  costs 
associated  with  base  repair  of  the  LRU  and  CIRF  repair  of  the 
SRU.  The  JUMBO  capacity  is  used  to  exclude  undesirable 
decisions  such  as  scrapping  or  depot  repairing  the  LRU  and 
CIRF  repairing  the  SRU.  The  second  arc  is  associated  with 
the  scrap  decision  for  the  LRU  and  as  such  will  have  a 
capacity  of  JUMBO  if  the  first  pass  decision  was  not  scrap. 
If  the  decision  was  scrap,  then  the  the  capacity  of  this  arc 
is  the  cost  associated  with  the  scrap  decision. 

The  CIRF  SRU  arcs  are  assigned  next.  The  first  arc 
connects  the  SRU  to  the  CIRF  LRU  node.  The  capacity  of  this 
arc  is  DEC  2  and  it  corresponds  to  the  pipeline  costs  of 
replacing  the  SRU  when  the  LRU  is  CIRF  repaired  and  the  SRU 
is  scrapped.  The  next  arc  originating  from  the  CIRF  SRU  is 
connected  to  the  NRLA  SRU  decision  node.  As  with  the  LRU 
above,  this  arc's  capacity  will  be  JUMBO  if  the  SRU  decision 
was  not  scrap.  If  the  decision  was  scrap,  then  a  check  is 
made  to  determine  which  decision  was  made  for  the  LRU.  If 
the  LRU  decision  was  base  repair,  then  the  DEC  2  costs  are 


added  to  the  SRU  scrap  cost  and  this  is  used  as  the  arc 


capacity.  I f  the  LRU  decision  was  not  base,  then  the  SRU 
scrap  cost  alone  is  used  as  the  arcs  capacity. 

The  next  arcs  entered  are  those  associated  with  the  NRLA 
LRU  decisions.  The  -first  arc  connects  the  LRU  node  to  the 
sink.  I-f  the  NRLA  decision  was  scrap,  then  the  capacity  of 
this  arc  is  JUMBO.  If  it  was  base  or  depot,  then  the  corre¬ 
sponding  LRU  repair  cost  is  used  as  its  capacity.  Following 
this,  the  JUMBO  capacity  arcs  connecting  the  LRU  to  the  SE 
are  assigned.  The  only  arcs  assigned,  connect  the  LRU  to  the 
SE  associated  with  the  NRLA  LRU  repair  decisions.  As  an 
example,  if  the  LRU  decision  was  base  repair,  then  the  arcs 
connecting  the  LRU  node  to  the  base  SE  would  be  assigned. 

The  NRLA  SRU  arcs  are  assigned  next.  The  first  arc 
assigned  connects  the  NRLA  SRU  to  the  NRLA  LRU.  This  arc  is 
assigned  in  only  two  instances;  (1)  if  the  NRLA  decisions 
were  base  repair  of  both  the  SRU  and  LRU  and  (2)  if  the  NRLA 
decisions  were  depot  repair  of  both.  Uhen  the  decisions  are 
base  repair,  the  arc  capacity  is  JUMBO.  This  ensures  that 
the  SRU  is  not  base  repaired  when  the  LRU  is  CIRF  repaired. 
When  the  decision  is  depot  repair  of  both,  then  the  arc 
capacity  is  the  sum  of  the  DEC  i  and  DEC  2  costs.  This  adds 
the  appropriate  DEC  costs  for  CIRF  LRU  repair  and  depot  SRU 
repair.  Following  this,  the  arcs  to  the  sink  are  added.  The 
capacity  of  these  arcs  are  JUMBO  if  the  NRLA  decision  was 
scrap  the  SRU.  If  both  the  LRU  and  SRU  are  repaired  at  the 
same  facility  the  corresponding  SRU  repair  cost  is  used  as 
the  capacity  for  this  arc.  In  the  instances  where  the  LRU  is 
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decisions  are  scrap  or  base  repair  and  the  SRU  decision  is 
depot  repair,  then  the  arc  capacity  is  the  sum  of  the  SRU 
depot  repair  cost  plus  the  DEC  1  and  DEC  2  costs.  The  last 
arcs  added  that  originate  from  the  NRLA  SRUs  are  the  JUMBO 
capacity  arcs  to  the  appropriate  SE  resources. 

The  first  arcs  added  to  the  network  connect  the  base  and 
depot  SE  resources  to  the  sink.  The  capacity  of  these  arcs 
are  the  SE  costs. 

Although  the  above  set  of  network  structures  is  fairly 
complicated,  it  does  provide  a  unique  set  of  repair  level 
decisions  for  the  DEPOT-SCRAP-BASE-CIRF  repair  options.  A 
vital  element  of  this  structure  is  the  standardization  of  the 
CIRF  repair  nodes  and  arcs  on  the  left  side  of  the  network. 
This  allows  the  output  subroutine  to  identify  changes  in  the 
repair  level  decisions  by  just  checking  the  state  of  the  CIRF 
LRU  and  SRU  nodes.  If  these  nodes  are  in  state  1 
(unlabeled),  then  the  LRU  or  SRU  optimal  repair  level  deci¬ 
sion  is  CIRF.  If  the  nodes  are  in  state  2  (labeled)  then  the 
original  NRLA  decision  is  optimal. 

CIRF  Network  Cut  Sets.  The  differences  in  the  structures 
of  the  subnetworks  lead  to  different  cut-sets  for  each  struc¬ 
ture.  Figures  14(a)  thru  14(g)  depict  the  feasible  cut-sets 
possible.  Table  V  presents  the  costs  associated  with  these 
cuts.  The  same  restrictions  apply  to  these  cut  sets  that  ap¬ 
plied  to  the  original  NRLA  cut  sets.  The  JUMBO  capacity  arcs 
cannot  be  part  of  the  cut  set,  and  the  cut  arcs  must  origi¬ 
nate  at  a  labeled  node  and  terminate  at  an  unlabeled  one. 


2b  2a 


<c>  Scrap/Scrap 


6c  6a 


<f>  Base/Depot 


7b  7a 


<g)  Scrap/Depot 


TABLE  L'l 


CIRF  Network  Decision  Costs 


Decision 

I  n  c  1 

LRU 

5RU 

DEP 

BE 

BASE 

CIRF 

Depot 

Depot 

X 

CIRF 

Depot 

X 

X 

CIRF 

CIRF 

x 

Base 

Base 

X 

Base 

CIRF 

X 

X 

CIRF 

CIRF 

X 

Scr  ap 

CIRF 

CIRF 


Depot  Scrap 
CIRF  Scrap 
CIRF  CIRF 


Base 


icrap 


Base  CIRF 
CIRF  Scrap 
CIRF  CIRF 


Base  Depot 
Base  CIRF 
CIRF  Depot 
CIRF  CIRF 


Scrap  Depot  X 
CIRF  Depot  X 
CIRF  CIRF 


The  inclusion  of  a  CIRF  capability  in  the  NRLA  program 


required  extensive  changes  to  the  program.  Subroutines 
OUTPUT  and  SECMP  were  modified  to  accept  and  work  with  the 
new  SE  codes.  These  subroutines  also  identify  when  CIRF  is 
included  in  the  analysis  so  that  the  appropriate  calculations 
can  be  performed.  The  extensive  reports  prepared  by  NRLA 
were  a  significant  hinderance  in  achieving  complete  program 
integration.  The  code  listing  for  OUTPUT  and  SECMP  are 
contained  in  Appendixes  G  arid  H.  Subroutine  CRFSET  creates 
the  new  network  structures  and  a  listing  of  it  is  contained 


in  Appendix 

C.  Once  again 

al  1 

changes  and  code 

were 

accomp 1 i shed 

using  FORTRAN  IV 

for  compatibility. 

I n  add i t i on 

to  the  above 

major  changes, 

numerous  smaller 

changes 

were 

made  in  subroutines  LRUCMP, 

FMCMP, 

and  SRUCMP 

The 

MAIN 

program  also 

required  changes  to 

i ncorporate 

CIRF. 

Th  i  s 

portion  of  the  program  also  has  a  significant  report  writing 
capability  which  required  modification. 

Ver i f i cat i on/Val idation 

The  verification  of  subroutine  CRFSET,  which  builds  the 
CIRF  networks,  was  accomplished  in  two  stages.  The  first 

stage  involved  the  development  of  a  set  of  network  data  which 
required  the  subroutine  to  build  ail  possible  network 
structures.  The  second  stage  of  the  verification  required 
running  the  CRFSET  subroutine  in  conjunction  with  SETNFT  *•  ' 
NMXFLO  to  ensure  that  the  input/output  interphases  wor>e: 
expected  and  to  further  identify  any  potential  prob'e"  , 
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The  data  set  developed  -for  this  test  included  all  possible 
cut  set  combinations  so  that  all  potential  logic  errors  could 
be  identified. 

The  data  set  developed  for  the  first  stage  was  designed 
to  build  all  seven  possible  network  structures.  The  arcs  in 
these  seven  structures  were  assigned  capacities  equal  to 
their  arc  number.  The  arc  numbers  were  based  on  the  expected 
arc  assignments  from  the  developed  program  logic.  This  set 
was  then  read  into  a  test  program  and  CRFSET  was  called  and 
implemented.  The  resulting  network  structures  were  then 
compared  with  the  expected  results  for  accuracy  of  execution. 
When  it  was  verified  that  the  program  was  building  the 
networks  correctly  the  second  stage  of  testing  began. 

The  second  stage  required  three  data  sets  to  test  the 
entire  series  of  cut  set  combinations.  During  this  test  the 
data  sets  were  read  into  the  test  program,  and  processed  by 
SETNET  to  build  the  original  RLA  network.  This  network  was 
then  solved  by  NiXFLO  and  the  results  printed  out  and  checked 
to  ensure  that  the  input  data  set  was  producing  the  expected 
set  of  repair  level  decisions  for  input  to  CRFSET. 
Following  this,  CRFSET  was  called  to  build  the  CIRF  network 
structures.  These  were  then  input  to  NMXFLO  for  solution. 
The  results  of  the  network  max-flow  solution  was  printed  out 
and  compared  with  the  projected  results.  The  results  of 
these  tests  conclusively  proved  that  subroutine  CRFSET  was 
processing  the  input  data  and  producing  the  correct  network 
structure  in  all  cases.  The  final  task  dealing  with  the  CIRF 
capability,  was  to  integrate  it  into  the  NRLA  program.  This 


task  proved  to  be  much  more  complicated  than  initially 
ant i c i pated. 


There  were  two  primary  -factors  which  increased  the 
complexity  o-f  the  integration  project.  The  first  factor  was 
the  extensive  work  required  to  update  the  report  capabilities 
of  the  program  to  include  CIRF.  This  included  changing  much 
of  the  logic  of  the  code  so  that  the  new  SE  codes  could  be 
accepted.  The  second  factor  which  is  directly  tied  to  this 
effort  is  the  user  friendliness  of  the  CREATE  diagnostic 
system.  On  several  occasions  during  the  debugging  of  the 
program  it  was  noted  that  the  diagnostic  outputs  from  the 
CREATE  system  shed  little  light  on  the  nature  of  the  problem 
encountered.  Logic  error  reports  were  usually  obscure  or 
non-existent.  In  one  instance  an  entire  subroutine  was 
missing  from  the  program  being  tested,  yet  no  recognizable 
diagnostic  error  showed  up  to  indicate  this  deficiency. 

Validation  of  the  CIRF  repair  level  option  requires  final 
integration  of  CIRF  into  the  NRLA  model.  In  addition,  the 
cost  equations  developed  for  CIRF  repair  must  be  validated  by 
the  appropriate  cost  analysts  in  AFALC. 

Anal ysi s/Resul ts 

The  final  analysis  of  the  CIRF  capability  requires  total 
integration  into  the  model  and  validation  of  the  cost 
equations  developed.  The  capability  to  include  CIRF  in  the 
RLA  networks  is  completely  verified  and  will  provide  a 
significant  analytical  tool  for  logistics  analysts  when  it 
has  been  integrated  into  the  NRLA  model. 
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VI I .  Cone  1  usi  ons/Recommendat  i  ons 


General  Comments. 

The  NRLA  program  is  a  unique  application  of  networking 
to  solve  the  repair  level  analysis  problem.  It  is  extensively 
used  by  the  Air  Force  logistics  community  to  aid  in 
determining  the  optimal  repair  level  decisions  for  the 
components  of  new  systems  being  procured  by  the  Air  Force. 
It  is  anticipated  that  models  of  this  type  will  continue  to 
play  an  increasingly  important  role  in  this  process  and  as 
such,  improvements  to  these  models  provide  a  significant 
benefit  to  the  Air  Force. 

During  the  course  of  this  thesis  effort  several  areas  of 
the  NRLA  program  were  identified  as  being  candidates  for 
further  study  and  modification.  In  general,  it  is  felt  that 
the  current  implementation  of  the  program  should  be  rewritten 
using  a  more  current  version  of  FORTRAN.  This  will  provide 
significant  benefits  in  the  versatility  of  the  programming 
techniques  available  for  further  enhancements  to  the  model. 
The  sensitivity  analysis  section  of  the  model  requires  study 
to  determine  if  it  can  be  implemented  with  the  CIRF 
capability.  The  data  structures  of  the  model  provide  an  area 
for  future  enhancements  using  new  and  more  efficient 
techniques  such  as  those  found  in  Jensen  and  Barnes  (4). 
Modification  of  the  data  structures  will  require  a  major 
revision  of  the  model  and,  therefore,  should  provide  an 
excellent  opportunity  to  upgrade  the  model  to  the  newer 


version  of  FORTRAN  as  recommended  above. 


Recommendat i ons 


MARLA .  Although  the  MARLA  enhancement  has  been  verified 
and  integrated  into  the  model,  conflicting  repair  decisions 
were  experienced.  It  appears  that  these  inconsistencies  can 
be  resolved  with  further  effort.  Based  on  the  possible 
efficiencies  that  could  be  realized  from  the  use  of  MARLA  as 
a  preprocessor  to  the  NRLA  model,  it  is  recommended  that  this 
effort  be  undertaken.  This  recommendation  is  compatible  with 
the  Air  Force  guidance  as  outlined  in  AFLC/AFSC  Pamphlet  800- 
4  (13).  Following  this,  MARLA  should  be  full/  tested  by  the 
users  to  determine  the  extent  of  the  benefits  available  from 
its  use.  For  very  large  systems  MARLA  may  be  the  best  tool 
available  to  reduce  the  size  of  the  data  input  to  NRLA.  In 
any  case  MARLA  can  be  utilized  as  a  separate  program  to 
perform  marginal  analysis  by  developing  i t  as  a  program  for 
use  on  micro  computer  systems. 

Label ino  Improvements.  The  Depth  First  Search  Labeling 
algorithm  developed  for  use  in  the  NMXFLO  subroutine  is  an 
original  variation  of  currently  accepted  networking 
algorithms.  Based  on  the  results  achieved  during  this  study 
a  high  level  of  confidence  was  developed  in  its  ability  to 
outperform  the  existing  procedure  as  well  as  other 
alternative  methods.  From  this  perspective  it  is  recommended 
that  this  enhancement  be  immediately  implemented  for  use  in 
the  NRLA  model  and  distributed  to  current  users  of  the  model 
as  a  modular  replacement  for  the  MAXFLO  subroutine.  The  use 
of  this  subroutine  should  significantly  reduce  the  computer 


resources  expended  during  NRLA  analysis,  thus  saving  the  Air 
Force  time  and  money. 

Cl RF  The  ability  to  analyze  a  CIRF  repair  option  is  a 
sorely  needed  capability  in  the  logistics  planning  area.  By 
using  a  new  series  of  network  structures  an  answer  has  been 
formulated  to  this  problem.  Verification  of  this  effort  has 
been  completed  and  full  integration  into  the  NRLA  program  is 
recommended  at  this  time.  A  special  consideration  related  to 
the  CIRF  effort  is  the  validation  of  the  cost  equations  by 
qualified  cost  analysts  in  AFALC. 


Cone  1  us  i  or 


In  summary,  for  the  integrated  logistics  support  process 
to  be  successful,  requires  that  repair  level  analysis  be 
skillfully  performed.  Without  question,  the  problem  solving 
technique  employed  in  the  NRLA  model  produces  economically 


based  optimal  repair  level  recommendations.  However,  it 
should  be  remembered  that  the  validity  of  these  computer 


generated  decisions  is  directly  dependent  upon  the  accuracy 
of  the  data  input  and  the  cost  estimating  relationships  used. 

By  accomplishing  the  stated  thesis  objectives, 
substantial  progress  has  been  achieved  in  increasing  both  the 


efficiency  and  the  capabilities  of  the  NRLA  model.  Because 
of  these  improvements,  it  is  now  appropriate  to  investigate 


methods  to  incorporate  the  other  non-economic  factors  which 
influence  the  repair  level  decision.  A  possible  approach 
could  be  the  development  of  an  objective  function  consisting 
of  the  elements  <  economic  and  non-economic)  which  are 


identified  as  critical  to  the  systems  performance.  This 
fuctional  relationship  could  then  be  used  to  establish  the 
appropriate  "cost"  for  each  arc's  capacity  in  the  network. 
In  this  way  such  factors  as  a  system's  mobility  and 
deployment  requirements,  availability  of  skilled  maintenance 
personnel,  or  operational  readiness  requirements  could  be 
considered  as  part  of  the  analytical  optimization  process 
used  by  NRLA.  Proper  integration  of  these  effects  would 
better  reflect  the  total  dynamics  of  the  multi-item,  multi¬ 
indenture,  multi-level  system  which  has  been  studied  during 
this  effort.  If  these  recommendations  are  implemented  a  more 
responsive  and  versatile  tool  will  be  available  to  Air  Force 
decision  makers  in  the  planning  and  management  of  the 
logistical  support  structures  necessary  for  tomorrow's 
weapons  systems. 


CAA 


CIRFLC 

CIRFSC 

CLR 

CMMH 

CMiHS 

CRCTC 

CRCTO 

CRCTPL 

CRCTSL 

HRAVSE 


Appendix  As  Glossary  of  Variables 


-  available  work  time  per  month  -for  a  centralized 
intermediate  level  maintenance  man  (man-hours/ 
month,  Maintenance  System  Data  Record). 

-  an  array  of  values,  one  for  each  LRU  failure  mode, 
each  of  which  is  the  sum  of  certain  life  cycle 
logistics  costs  associated  with  intermediate  level 
repair  of  the  LRU. 

-  an  array  of  values,  one  for  each  SRU,  each  of  which 
is  the  sum  of  certain  life  cycle  logistics  costs 
associated  with  intermediate  level  repair  of  the 
SRU. 

-  hourly  labor  rate  for  centralized  intermediate 
level  maintenance  man  (*/hour,  Maintenance  System 
Data  Record) . 

-  the  number  of  maintenance  man-hours  required  for 
repair  of  an  LRU  if  the  repair  is  done  at 
Cl RF  level  (man-hours/repair,  LRU  Failure  Mode  Data 
Record) . 

-  the  number  of  maintenance  man-hours  required  for 
repair  of  an  SRU  if  the  repair  is  done  at 
the  CIRF  level  (man-hours/repair,  LRU  Failure  Mode 
Data  Record.) 

-  the  elapsed  time  from  removal  of  a  failed  LRU  at 
CONUS  CIRF  until  the  item  could  become  a  service¬ 
able  spare  in  depot  stock,  it  includes  the  time 
required  for  base  to  depot  transpor tat i on  and  the 
depot  shop  flow  time  required  for  repair  (months, 
LRU  Data  Record) . 

-  the  elapsed  time  from  removal  of  a  failed  LRU  at  an 
overseas  CIRF  until  the  item  could  become  a 
serviceable  spare  in  depot  stock  (months,  LRU  Data 
Record) . 

-  the  expected  number  of  unserviceable  LRU  assets  in 
the  CIRF  repair  pipeline  (No. LRU, Computed) . 

-  the  number  of  spare  LRUs  to  be  purchased  to  satisfy 
LRU  demands  expected  to  occur  during  the  CIRF 
repair  cycle  time  (No.  LRUs, Computed) . 

-  an  array  for  each  item  of  support  equipment  which 
stores  the  available  equipment  hours. 


NAPATH  -  an  array  which  is  used  to  store  the  arcs  on  the 
path  from  the  source  to  the  sink. 

NTCLFM  -  an  array  of  values,  one  for  each  LRU  failure  mode, 

each  of  wh i ch  is  the  minimum  number  of  CIRF  level 

maintenance  personnel  to  be  trained  for  the  LRU 
repair  task  (number  of  people,  LRU  Failure  Mode 
Data  Record) . 

NTCS  -  an  array  of  values,  one  for  each  SRU,  each  of  which 

is  the  minimum  nimber  of  CIRF  level  maintenance 

personnel  to  be  trained  for  the  SRU  repair  task 
(number  of  people,  SRU  Data  Record) 

SCRCTC  -  SRU  CIRF  repair  cycle  time  for  CONUS  bases;  the 
elapsed  time  from  removal  of  a  failed  SRU  (from  the 
LRU)  at  a  CONUS  base  until  the  item  could  become  a 
serviceable  spare  in  depot  stock  (months,  SRU  Data 
Records) . 

SCRCTL  -  the  number  of  spare  SRUs  to  be  purchased  to  satisfy 
SRU  demands  expected  to  occur  during  the  SRU  CIRF 
repair  cycle  time  (No. SRU, Computed) . 

SCRCTO  -  SRU  Cl Rf  repair  cycle  time  for  overseas  bases;  the 
elapsed  time  from  removal  of  a  failed  SRU  from  the 
LRU  at  an  overseas  base  until  the  item  could  become 
a  serviceable  spare  in  depot  stock  (months, SRU  Data 
Record) . 

SCRCTP  -  the  expected  number  of  unserviceable  SRU  assets  in 
the  CIRF  repair  pipeline  (No. SRU, Computed) . 

SEARY  -  a  support  equ i pment/LRU/SRU  cross  reference  array, 
which  stores  the  pointers  of  the  SE  for  a 
particular  failure  mode  or  SRU. 

TBCST  -  a  cost  array  which  collects  the  total  cost  of 

repair  at  the  base  level  for  an  LRU  or  SRU. 

TCCST  -  a  cost  array  which  collects  the  total  cost  of 

repair  at  the  CIRF  for  an  LRU  or  SRU. 

TDCST  -  a  cost  array  which  collects  the  total  cost  of 

repair  at  the  depot  level  for  an  LRU  or  SRU. 

TFD  -  training  factor  for  depot;  the  expected  number  of 

times  that  formal  maintenance  training  will  be 
required  for  depot  personnel  (dimensionless, 
Computed) . 

-  total  life  cycle  repair  demands  for  an  LRU  at  each 
base  (No. repair  demands, /L. C. , Computed) . 


TLCD 


TL  CDF  -  total  1  i  -f  e  cycle  repair  demands  -for  a  particular 
failure  mode  of  an  LRU  at  each  base  (No. repair 
demands, /L. C. ,  Computed). 

TRB  -  annual  turnover  rate  for  intermediate  level  maint¬ 

enance  personnel  (fraction  of  personnel 
replaced/year,  Maintenance  System  Data  Record). 

TRC  -  the  expected  training  cost,  instruction  and  mater¬ 

ials,  for  LRU  repairs  (  ♦/Man/week,  LRU  Failure 
Mode  Data  Record) . 

TSECB  -  a  cost  array  used  to  collect  the  SE  costs  at  the 
base  level  for  a  particular  LRU  or  SRU. 

TSECC  -  a  cost  array  used  to  collect  the  SE  costs  at  the 
CIRF  for  a  particular  LRU  or  SRU. 

TSECD  -  a  cost  array  used  to  collect  the  SE  costs  at  the 
depot  level  for  a  particular  LRU  or  SRU. 

TSECST  -  the  total  cost  associated  with  the  acquisition  and 
operation  of  a  particular  piece  of  support  equ ip- 
men  t . 


m 


FMCMP  - 


OUTPUT  - 


RESET  - 


SECMP  - 


Appendix  6:  NRLA  Subroutine  Descriptions 


DECIDE  - 


LRUCMP  - 


MAXFLO  - 


MTBFSA  - 


this  subroutine  is  called  -from  UCLSA  and  MTBFSA. 
Its  purpose  is  to  specifically  determine  the 
decision  changes  identified  by  the  sensitivity 
analysis  routines  and  document  the  changes  on  an 
output  file. 

this  subroutine  is  called  by  MAIN,  MTBFSA,  and  UCLSA 
to  compute  repair  level  option  costs  for  an  LRU 
failure  mode.  These  costs  are  saved  in  arrays 
DEPOLC,  SCRPLC,  CIRFLC,  and  BASELC.  The  expected 
number  of  maintenance  man-hours  required  monthly  for 
the  failure  mode  is  computed  for  depot,  CIRF,  and 
base  and  saved  in  FSEUHD ,  FSEUHC,  and  FSEUHB 
respect i vel y . 

this  subroutine  is  called  by  MAIN,  MTBFSA  and  UCLSA. 
Its  function  is  to  compute  inventory  stock  levels 
and  life  cycle  SCRAP  option  costs  for  an  LRU. 

this  subroutine  is  called  by  MAIN,  MTBFSA,  and 
UCLSA.  Its  purpose  is  to  determine  the  maximum  flow 
through  the  network  and  the  minimum  cut  set.  This 
cut  set  identifies  the  minimum  cost  set  of  repair 
level  decisions. 

this  subroutine  is  called  from  MAIN  to  determine  the 
effects  of  changes  in  the  MTBF  for  an  LRU.  The 
effects  could  be  changes  in  failure  mode  repair 
level  decisions,  SRU  repair  level  decisions,  and/or 
changes  to  SE  decisions. 

this  subroutine  is  called  from  MAIN  to  print  the 
optimal  solution  results.  This  includes  support 
equipment  decisions  plus  LRU  and  SRU  repair  level 
decsi si ons. 

this  subroutine  is  called  by  MAIN,  MTBFSA,  and 
UCLSA.  Its  purpose  is  to  exploit  the  known 
structure  of  the  RLA  network  so  that  advance  flow 
can  be  placed  on  the  network  arcs  to  reduce  the 
solution  time  for  the  MAXFLO  subroutine. 

this  subroutine  is  called  by  MAIN  and  MTBFSA.  It 
determines  the  quantity  of  each  SE  resource 
potentially  required  at  depot,  CIRF,  and  base  level 
plus  the  life  cycle  cost  for  the  resources. 
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SETNET  - 


subroutine  is  called  by  MAIN  after 


costs 


and  dependenc''  relationships  are  known.  The 
function  of  SETNET  is  to  use  this  information  to 
construct  a  matrix  representation  of  the  RLA 
network,  determine  and  save  pointers  for  use  by  the 
solution  algorithm  and  save  pointers  for  use  by -the 
sensitivity  analysis  subroutines. 


SORT 


SRUCMP  - 


this  subroutine  is  called  from  SETNET.  As  ists  name 
implies,  its  function  is  to  sort  values  supplied  by 
SETNET.  Specifications  to  SORT  are  contained  in  its 
arguments  list. 

this  subroutine  is  called  by  MAIN,  MTBFSA,  and 
UCLSA  to  compute  SRU  related  life  cycle  costs. 
These  costs  are  saved  in  the  SRU  related  arrays 
DEPOSC,  SCRPSC ,  CIRFSC,  BASESC,  BRLSS,  BRLDS,  and 
BRLCS. 


Appendix  C:  CRFSET  Subroutine 


10  SUBROUTINE  CRFSET(MS0L2 ) 

20C 

30C 

40  ICIRF=2 

50C 

60C  CODE  FROM  STATEMENT  l‘  2000  THROUGH  I)  2210  UTILIZES  THE  ABOVE 

70C  COMPUTED  ITEM  $,  SE  COSTS  TO  PRODUCE  AN  RLA  COST  NETWnRK . 

80C 

90C  AS  A  PRELIMINARY  STEP  DETERMINE  NODE  NUMBERS  FOR  KEY  NODES  IN 

100C  THE  NETWORK.  THE  VARIABLES  ARE  CODED  WITH  THE  PREFIXES  'LC' 

110C  FOR  'LAST  CIRF'  AND  'LN'  FOR  'LAST  NRLA  INPUT' 

120C 

130  2000  LCSE-NCSE+1 
140  LCLRU-LCSE+LFMS 

150  LCSRU-LCLRU+SFMS 

160  LNLRU-LCSRU+LFMS 

170  LNSRU-LNLRU+SFMS 

180  LNDSE-LNSPU+NDSE 

190  LNBSE=LNDSE+NBSE 

200  LNODE-LN3SE+1 

210  I F ( LNODE .LE . MAXNOD )  GO  TO  2005 

220  WRTTE(6,2002) 

230  2002  FORMAT ( ' 0  NODE  VECTORS  ARE  TOO  SMALL  —  STOP') 

240  STOP 

250C 

260C 

270C  SET  !■  OF  ARCS  TO  ZERO  &  2  OF  NODES  TO  1 

280  2005  NARCS-0 

290  NNODF.S-1 

300  NBSEPO-NDSF.+NBSE 

310  NBSFP1 =NDSE+NBSE+l 

320C 

330  IF(NCSE.LT.l)  GO  TO  2015 

340C  CREATE  ARCS  FROM  SOURCE  TO  CIRF  SF  —  CIP.F  SE  COST  ARCS 

350  DO  2010  I2010-MBSEP1 .NUMSER 

360  NAPCS-NARCS+1 

370  NNODES-NNODES+l 

380  SRCE (NARCS ) *1 

390  DEST (NARCS )-NNODES 

400  CAP (NARCS )-SEC0ST( 12010) 

410  2010  CONTINUE 
420C 

430C  CREATE  ARCS  FROM  THE  SOURCE  TO  THE  CIRF  LRU  FAILURE  MODE 

440C  NODES  —  ITEM  RELATED  COSTS  FOR  LRU  REPAIR  AT  CIRF 

450  2015  DO  2020  12020-1, LFMS 

460  NAPCS»NARCS+1 

4  70  "INODES -NNODES+1 

480  SRCE (NARCS )-l 

490  DEST (NARCS )-NN0DES 

500  CAPCiARcsi-rTPFT  riT?n?n> 


510  2020  CONTINUE 
520C 

330C  CREATE  ARCS  FROM  THE  SOURCE  TO  THE  CIPF  SRU  NODES  —  ITEM 

540C  RELATED  COSTS  FOR  SRU  REPAIR  AT  CIRF 

550  IF(SFMS.LT.l)  GO  TO  2040 

560  DO  2030  12030=1  ,SF^S 

570  NARCS *NARCS+1 

5S0  NNODES=NNODF.S+1 

590  SRC E (NARCS )=l 

600  DEST (MARCS )=MMODES 

6LO  CAP (MARCS )=CIRFSC(I2030) 

620  2030  CONTINUE 
630C 

640C  CREATE  ARCS  FROM  THE  CIRF  SE  MODES  TO  THE  CIRF  LRU/SRU  MODES 

650  2040  IF(NCSE.LT.l)  CO  TO  2070 
660C  FOR  EACH  CIRF  SE  — 

670  DO  2060  I2060-NBSEP1 .NUMBER 

680C  SEPF  CONTAINS  A  POINTER  TO  THE  SF  CROSS  REFERENCE  (SEXREF) 

690  IPTR=SEPF(I2060) 

700  IF(IPTR.EQ.O)  CO  TO  2060 

7 IOC  SEXREF  CONTAINS  THE  *  OF  AN  LRU  FAILURE  MODE  (ITEM  >  0)  OR 

720C  THE  NUMBER  OF  AN  SRU  (ITEM  <  0) 

730  2050  ITEM-SEXREF(IPTR) 

740  IF(ITEM.LT.O)  ITEM-LFMS-ITEM 

750  NARCS =NARCS+1 

760  SRC E (NARCS )=I 206 O-NBSEPO+l 

770C  ITEM  IS  A  POINTER  TO  AN  LRU  OR  SRU  AND  WHEN  ADDEn  TO  LCSF. 

780C  GIVES  THE  APPROPRIATE  NODE  NUMBER 

790  DEST (NARCS )=LCSE+ITEM 

800  CAP (NARCS ) =JUH30 

8  IOC 

820C  NXTITM  IS  A  POINTER  TO  THE  NEXT  ENTRY  IN  SEXREF  WHICH  IS  AN 

830C  LRU  OR  SRU  REQUIRING  THE  CURRENT  SE,  I.E.,  SF  V  12060. 

840  NPTR-NXTITM(IPTR) 

850  IF(NPTR.EO.O)  CO  TO  2060 

860  IPTR-NPTR 

870  GO  TO  2050 

880  2060  CONTINUE 

R90C 

900C  CREATE  THE  ARCS  EMANATING  FROM  THE  CIRF  LRU  FAILURE  MODE  NODES 

910  2070  DO  2080  12080=1 ,LFMS 
920C  CHECK  FOR  AN  SRU 

930  IPTR=SRUPTR( T2080) 

940  IF(IPTR.EO.O)  CO  TO  2075 

950C  CREATE  AN  ARC  TO  THE  CIRF  SPU  NODE  —  COSTS  UNIOUE  TO  BASE 

960C  REPAIR  OF  LRU  6.  CIRF  REPAIR  OF  SRU 

970  NARCS “NARCS+1 

980  SRCF (NARCS )=LCSE+I 2080 

990  DEST ( NARCS ) “LCLRU+IPTR 

1000  CAP (NARCS ) =3RLCS ( IPTR) 

101 o  IF  (OPDFCL ( 1 20^0) .NE.LOCAT ( 3 ) )  CAP (NARCS ) = JUMBO 

1020C  CREATE  AN  ARC  TO  THE  URL \  LRU  NODES  —  COST  OF  SCRAPPING  THE  L 


1030  20  75  NARCS-NARCS+l 
1040  SRCE (NARCS ) -LCS E+I2080 

1050  DEST (MARCS ) -LCS RU+I 2080 

1060  CAP (MARCS )=SCRDLC (12080) 

1070  IF  (OPDECL(I20R0) .ME.LOCAT(2) )  CAP (MARCS)* JUMBO 

1030  2080  CONTINUE 

1090C 

1100C  CREATE  ARCS  FROM  THE  CIRF  SRU  NODES  TO  TOE  NRLA  SRU  NODES  - 

11  IOC  COST  OF  SCRAPPING  THE  SRU 

1120  IF(SFMS.LT.l)  GO  TO  2100 

1130  DO  2090  12090-1, SFMS 

1140C  CREATE  ARCS  FROM  THE  CIRF  SRU  MODES  TO  THE  CIRF  LRU 

1150C  NODES,  TO  ELIMINATE  INFEASIBLE  LRU /SRU  REPAIR  LEVEL 

1160C  MATCHES  SUCH  AS  LRU-D,  SRU-C;  LRl’-S ,  SRU-C;  LRU-C,  SRU-B 

1170  IF  (OPDECS ( 12090 ) . EO.LOCAT(2 ) . AND .OPDECL(LRUPTR( 12090 ) ) . 

1180  &EQ.LOCAT( l ))  GOTO  2082 

1190  IF  (OPDECS (12090) .EQ.LOCAT(2) . AND .0?DECL(LRUPTR( 12090 ) ) .EO. 

1200  &LOCAT(2) )  GOTO  2082 

1210  COTO  2035 

1220  2082  NARCS-NARCS+l 

1230  SRCE (NARCS ) -LCLRU+I2090 

1240  DEST (MARCS ) -LCSE+LRUPTR (12090) 

1250  CAP (NARCS )*BRLSS (12090) 

1260C  CREATE  ARCS  FROM  CIRF  SRU  TO  NRLA  SRU  —  SCRAP  COSTS 

1270  2085  NARCS-NARCS+l 

12S0  SRCE (NARCS )-LCLRU+I2090 

1290  DEST (MARCS ) -LNLRU+I2090 

1300  CAP ( MARCS )*SCRPSC( 12090) 

1310  IF  (OPDECS ( 12090) .NE.LOCAT(2) )  CAP(NARCS) -JUMBO 

1320C  ADD  DEC2  COSTS  TO  SCRAP  CSOT  FOP.  BASE  SCRAP  DECISION 

1330  IF  (OPUECL(LRUPTR(I2090) ) . EO.LOCAT(3) . ANO.OPDECS(I2O90) .EO. 

1340  &LOCAT ( 2 ) )  CAP(NARCS)-CAP(NARCS)+BRLSS(I2090) 

1350  2090  CONTINUE 
1360C 

1370C  CREATE  ARCS  EMANATING  FROM  THE  NRLA  LRU  NODES 

1380  2100  DO  2120  12120-1, LFMS 

1390C  CREATE  AN  ARC  FROM  THE  NRLA  LRU  MODE  TO  THE  LAST  NOD  (SINK) 

1400C  ITEM  RELATED  COST  OF  NRLA  REPAIR  OF  THE  LRU 

1410  2110  NARCS-NARCS+l 

1420  S  P.CE  (NARCS  ) -LCS  RU+I  2 1 20 

1430  DEST (NARCS )«LNODE 

1440  CAP (NARCS) -JUMBO 

1450  IF  (OPDECL(I2120).EO.LOCAT(3))  CAP(NARCS)-BASELC(I2120) 

1460  IF  (OPDECL(I2120) .EQ.LOCAT ( 1 ) )  CAP(NARCS)=DEPOLC(I2120) 

1470  2120  CONTINUE 
1480C 

1490C  CREATE  ARCS  EMANATING  FROM  THE  MRLA  SRU  NODES 

1500  IF  (SFMS.LT.l)  GOTO  2140 

1510  DO  2130  12130-1, SFMS 

1520C  CREATE  AN  ARC  FROM  THE  NRLA  SRU  TO  MRLA  LRU  IF  THE 

1 5 3 )C  MRLA  DECISION  IS  EITHER  LRU-D,  S°U-D  OR  LRU-P  ,  SRU-B 

1540  IF  (OPDECL (L°UPTR( 12130)). FO. LOCAT( 1 ) . AND. OPDECS ( 1 2 1 30 ) . EO. 


1550  &LOCAT ( 1 ) )  GOTO  2132 

1550  IF  (0PDECL(LRUPTR(I2130))  .F.O.LOCAT(3)  .AND  .OPDF.CS  (12130)  .EO. 

1570  &LOCAT (3 ) )  GOTO  2132 

1580  GOTO  2135 

1590  2132  MARCS-NARCS+1 

1600  SRCE(MARCS)=LNLRU+I2130 

1610  DEST  (MARCS  )  -LCS  RU+I.RUPTR  (12130) 

1620  CAP (NARCS ) -JUMBO 

1630  IF  (OPOECS(I2130) .F0.LOCAT(3) )  GOTO  2135 

1640  CAP (11  ARCS )  -BRLSS  (121 30 )+RRLDS ( 12 1 30 ) 

1650C 

1660C  CREATE  AM  ARC  FROM  THE  N REA  SR!!  MODE  TO  THE  LAST  NODE  (SINK) 

1670C  ITEM  RELATED  COST  OF  NRLA  REPAIR  FOR  THE  SRU 

1630  2135  NARCS-NARCS+l 

1690  SRCE (MARCS )-LNLPU+I 21 30 

1700  DEST (MARCS )-LMODE 

1710  CAP (NARCS ) -JUMBO 

1720  IF(OPDECS  (12130)  .EO.LOCAT(3) )  CAP (NARCS )-RASESC (12 130) 

1730  IF  (OPDECS(I2130).EO.LOCAT(1))  CAP(NARCS)-DEPOSC(I2130) 

1740  IF  (OPDECL (LRUPTR  ( 1 2 1 30 ) ) . ME . LOCAT ( 1 )  . AMD . OPDF.CS  ( 1 2 1  30 ) . F.Q . 

1750  &LOCAT ( 1 ) )  CAP (NARCS ) -CAP ( MARCS )+BRLS  S ( 1 2 1 30 )+BRLDS ( 1 2 1 30 ) 

1760  2130  CONTINUE 
1770C 

1730C  CREATE  ARCS  TO  AND  FROM  THE  URL A  SE  NODES 

1790C  ADD  ARCS  FOR  THE  NRLA  DEPOT  SE  DECISIONS 

1800  2140  IF  (MDSE.LT.l)  GOTO  2160 
18 IOC  FOR  EACH  DEPOT  SE  — 

1820  DO  2150  12150-1 .NOSE 

1*30  IPTR-SEPF(I21 50) 

1340  IF  (IPTR.EQ.O)  GOTO  2149 

1850  2145  ITEM-S  EXREF ( IPTR ) 

1860  IF  (ITEM. LT.O. AND. OPDECS(TABS  (ITEM)  )  .NF..LOCAT(  1 )  )  GOTO  2148 

1870  IF  ( ITEM. GT.Q. AND. OPDECL (ITEM) . ME. LOCAT ( 1 ) )  COTO  2148 

1830  IF  (ITEM. LT.O)  ITF.m-LFMS-ITCM 

1390  NARCS-NARCS+l 

1900  SRCE (NARCS )-LCSRU+ITEM 

1910  DEST  (MARCS  )-LMSP.U+I 2 150 

1920  CAP (MARCS) -JUMBO 

1930  2 14*  NPTR-NXTITM(IPTR) 

1940  IF  (NPTR.EO.O)  GOTO  2149 

1950  IPTR-NPTR 

i960  GOTO  2145 

1970C  ARCS  FROM  DE°OT  SE  TO  SINK 

1980  2149  MARCS -MARCS+l 

1990  SRCE(MARCS)«LMSRU+I2150 

2000  DEST (MARCS )=LMODE 

2010  CAP (NARCS ) -SECOST (12150) 

2020  2150  CONTINUE 

2030C  ADD  ARCS  FOR  NPL\  BASF.  DECISIONS 

2040  2160  IF  (NBSE.LT. 1 )  GOTO  21*0 
2050C  FOR  EACH  BASF.  SE 

2060  MDSFP1 -NDSE+1 


2 n 70  DO  21  70  121  70-vDSEPl  .NESEPO 

2080  IPTR-SEPF(I2170) 

2090  IF  (IPTR.EQ.O)  GOTO  2169 

2100  2165  ITEM-SEXREF(TPTR) 

2110  IF  ( ITEM. LT.O .AND .OPDECS  ( IA8S  ( ITEM)  )  .NF.  .LOCAT (  3 ) )  GOTO  2168 

2120  IF  (ITEM. CT.9. AND. OPDF.CL  ( ITEM)  .NE. LOCAT( 3 ) )  ^OTO  216S 

2130  IF  ( TTFM.LT.O)  ITEM-LFMS-ITEM 

2140  NARCS-NARCS+1 

2150  s  rc  e  (narcs  )*lcsru+itf.m 

2 1  60  PEST  ('.’ARCS )  -LNSRU+I2  1  70 

2170  CAP  (  "ARCS  )  -JUMBO 

2180  2168  NPTR-NXTlTM(irTR) 

2190  IF  (UPTR.EO.O)  COTO  2160 

2200  IPTP-NPTR 

2210  COTO  2165 

2220C  ARCS  FROM  BASE  SE  TO  SINK 

2230  2169  NARCS-NARCS+l 

2240  SRCF(VARCS)-LMS:”.’+I2170 

2250  PEST (NARCS )»LNOPE 

2260  CAP (NARCS )-SEC0ST (12170) 

2270  2170  CONTINUE 
2280  2180  NNODES-LNODF. 

2290C 

2300C  EXPLICITLY  SORT  THE  SRCF.-DEST-CAP  TRIPLES  INTO  ASCENDING  ORDER 

23IOC  RY  SRCE  (FLOW  IS  USED  BY  SORT  TO  STORF.  ORDERING  POINTERS) 

2320  CALL  SORT(SRCE, DEST, CAP, FLOW, NARCS .MAXARC.l ) 

2330C  SAVE  POINTERS  FOR  FORWARD  SCAN  IN  FWDSP 

2340  DO  2190  12190-1 .NARCS 

2350C  UTILIZE  THE  VALUES  IN  SRCE  IN  REVERSE  ORDER.  (IF  -’ODE  *1*  IS  A 

2360C  SOURCE  FOR  r'l'  ARCS  THEN  FWPSP(l)  WILL  SUCCESSIVELY  GET  THE 

2370C  VALUES:  N,  N-l,  N-2 ,  .  .  .  ,2,  1.  SIMILARLY,  FNPSP(K)  LULL 

2380C  HAVE  THE  VALUE  'J'  WHERE  J  CORRESPONDS  TO  THE  FIRST  OCCURENCE 

2390C  OF  NODE  *  K  IN  SRCE.) 

2400  J2190«NApCS+l-I2190 

2410  K2190-SRCE( J2190) 

2420  FWPSP(K2190)-J2190 

2430  2190  CONTINUE 

2440  FWDSP(LNODE)«NARCS+i 

2450  LNODM 1 -LNODE- 1 

2460  DO  2195  12’  .LNODUl 

2470  J2195-LNO m 

2480  IF(FL’PSPf.j:  '  0)  FWDSP(J2195)-FNDSP(J2105+l ) 

2490  2195  CONTINUE 
2500C 

2510C  USE  SORT  7  .  .  r  ITLY  SEQUENCE  THE  DEST  ENTRIES  INTO  ASCENDING 

2520C  ORDER.  THE  ORDER  WILL  BE  SPECIFIED  BY  THE  ENTRIES  IN  BKPTR . 

2530C  THUS,  THE  ENTRIES  IN  BACKSP  ARE  POINTERS  TO  SRCE-DFST-CAP 

2540C  TRIPLES  BASED  ON  TUE  DF.ST  VALUES. 

2550  CALL  SORT (DEST ,0 ,0 , BKPTR .NARCS .MAXARC ,0 ) 

2560C  SAVE  POINTERS  FOR  BACKWARD  SCAN  IN  BACKSP 

2570  BACKSP ( 1  )*0 

2580  DO  2200  12200-1 .NARCS 


2590  J22O0=BKPTR(I22OO) 

2603  K2200=DEST(J2200) 

2610C  RACKS P (K)  IS  THE  POSITION  IN  BKPTR  WHICH  HOLDS  THE  LAST  POINTER 

2620C  TO  A  SRCE-nEST-CAP  TRIPLE  HAVING  NODE  *  K  AS  THE  DRST.  THE  FIRST 

2630C  POINTER  TO  A  TRIPLE  HAVING  NODE  11  K  AS  THE  DEST  IS  BKPTR (K-l  )  +  l . 

2640  BACKS? (K2 200 >*12200 

2650  2200  CONTINUE 
2660  DO  2222  12222*2 ,LNODE 

2670  IF(EACKSP(I2222) .EQ.O)  RACKSP(I2222) 

2680  2222  CONTINUE 
2690C 

2700C  IF(MSGL2.LT.l)  GO  TO  2205 
2710C  PRINT  2202, 'SRCF.  ' , (SRCE( I) ,1=1 .NARCS ) 

2720C  2202  FORMAT < '  ',A6,('  ',2013)) 

2730C  PRINT  2202, 'DEST  ' , (DEST(I) ,1=1 .NARCS) 

2740C  PRINT  2202, 'CAP  ' , (CAP(I) ,1=1 .NARCS) 

2750C  PRINT  2202, 'FLOW  ', (FLOW(I) ,1=1 .NARCS) 

2760C  PRINT  2202, 'BKPTR  ', (BKPTR(I) ,1=1 .NARCS) 

2770C  PRINT  2202  , ' FWDSP  ' ,  (FWDSP(I)  ,1=1  ,LNODE) 

2780C  PRINT  2202 , ' BACKSP' , (BACKSP (I ) , 1=1 ,LNODE) 

2790  2205  DO  2400  12400=1 .NARCS 
2800  J»SRCE(I2403) 

2810  K=DEST( 12400) 

2820  IF(J.NE.l)  GO  TO  2230 

2830  IF(K.GT.LCSE)  GO  TO  2210 

2840  SEARCP(K-1)=I2400 

2850  GO  TO  2300 

2860  2210  IF(K.GT.LCLPU)  GO  TO  2220 

2870  LDARC  (K— LCSF.)»I2400 

2880  GO  TO  2300 

2890  2220  S0ARC(K-LCLRU)»I2400 
2900  GO  TO  2300 

2910  2230  IF(J.LE.LCSE)  GO  TO  2300 
2920  IF(J.GT.LCLRU)  GO  TO  2250 

2930  IF(K.LE.LCSRU)  GO  TO  2240 

2940  LS ARC (J-LCSE) =12403 

2950  GO  TO  2300 

2950  2240  SBDARC(K-LCLRU)»I24O0 
2970  GO  TO  2300 

2980  2250  IF( J.GT.LCSRU)  GO  TO  2280 

2990  SSARC(J-LCLRU)=I2400 

3000  GO  TO  2300 

3010  2260  IF(J.GT.LNLRU)  GO  TO  2230 

3020  IF(K.NE.LMODE)  GO  TO  2270 

3030  LBARC(J-LCSRU)=T2400 

3040  GO  TO  2300 

3050  2270  IF(K.CT.LMSRU)  GO  TO  2300 

3060  SBSARC (K-LNLRU) =12403 

3070  GO  TO  2300 

3080  2280  IF(J.GT.LNSRU)  GO  TO  2290 

3090  IF(K.NE.LNODE)  GO  TO  2300 


BACKSP ( 12222 )  =  BACKSP( 1 2222-1  ) 


3100  SBARC(J-LNLRU)«I2400 

3110  GO  TO  2300 

3120  2290  SF.ARCP(J+NCSE-LNSRU)  *12400 

3130  2300  FLOW(I2400)=0 

3140  2400  CONTINUE 

3150C  IF(MSCL2.LT.2)  GO  TO  2440 

3160C  DO  2410  I* l ,LFMS 

3170C  PRINT, 'LPTRS  ' ,LDARC(I) ,LSARC(I) ,LBARC(T) 

3180C  2410  CONTINUE 
3190C  DO  2420  1-1 ,SFMS 

3200C  PRINT, 'SPTRS  * ,SDARC(I) ,SSARC(I) ,SBARC(t) ,SEDARC(I ) ,SBSARC(I) 

32 IOC  2420  CONTINUE 

3220C  DO  2430  1*1, NUMBER 

3230C  PRINT, 'SEPTRS  *,SEARCP(I) 

3240C  2430  CONTINUE 
3250C  2440  CONTINUE 
3260  RETURN 

3270  END 


* 
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Appendix  D:  MARLA  Subroutine 


subroutine  marla 

IMIS  PROGRAM  IMPLEMENTS  THE  CONCEPT  OF  ‘'AROI'LAL  * 
ANALYSIS  TO  DETERMINE  TI1E  OPTFTM  t'F.PMR  LEVEL  * 

FOR  AN  ITEM  (LRU  FAILURE  MODE  OR  SPU)  OF  A  WEAPONS  * 
SYSTEM.  IT  WILL  ACT  AS  A  PPEPPOCESSER  TO  THE  MAX  FLO  * 
SUBROUTINE  OF  THE  NRLA  PROGRAM.  * 

00  LOO?  10  CONDUCTS  THE  ANALYSIS  FOR  THE  FIRST  ITEM  * 
TO  THF.  LAST  ITEM  OF  THE  SYSTEM.  * 

TOTITU-LFMS+SF’S 
HO  10  I=*l ,TOTITM 
IF(I.CT.LFMS)  GO  TO  15 
T F ( L rMOO (I).NE.l)  CO  TO  12 
LFM0D(I)=2 
INDD»1000 

12  I F  (  L FMO B(I).NE.l)  GO  TO  14 
LFM03(I)a2 
INDB-10 

14  IF(LFMOC(I).NE.l)  GO  TO  15 

Lr'OC(I)=2 

INDC=*1 

16  IF(LFMOS  (I).NE.l)  GO  TO  IS 
LFl'OS  ( I  )*2 
INDS-100 

15  CONTINUE 

IF  (CIRF.EO.O)  LFU0C(I)=3 

00  LOOP  20  CHECKS  EACH  POSITION  OF  THF.  SE/LNU/SRU  * 
CROSS  REFERENCE  ARRAY  (SF.ARY)  FOR  THE  POINTER  VALUES  * 
OF  THE  SE  NECESSARY  TO  REPAIR  THF.  Ith  ITEM .  * 


►v* 

k 

Kt 

i 


DO  20  120*1 ,25 
PTR*SEARY ( I , 120 ) 

IF(PTR.EO.O)  CO  TO  20 
IO-SECOOE(PTR)/IO 

I F(LFUOD (I) .GT.O .OR. IP .GE . 5000 )  GO  TO  32 
3UY*AIMT ( 1 .0+(FSEU»n(I ) /OPHRS (PTR ) ) ) 

TSECD( I ) *TSECO( I )  -F(PUY*CADO (PTR)+FD3 ( PTR )  +  ( FSEUUD ( I ) / 
f,  OPHRS  (PTR)  )  *CODS( PTR) *PIUP) 

32  IF(LFMO3(I).GT.O.OR.ID.LT.5OO0.OR.Ir>.OE.9OnO)  GO  TO  34 
3UY-AINT(  l  .0+(FSF.U.lB(I )  /OPHRS  (PTR)  )  ) 

TSECP, (I  )*TSEC.R(T )+(BUY*C4DF.(PTR)+FPB(PTP )  +  (  ( FSFUHP (T )  / 

&  OPHRS (PTP ) ) *CODB(PTR) *P IUP ) ) 

34  IF(LF.MOC  (I).  GT.O.  OR. in.  LT  .9000 )  CO  TO  20 

SUY-P+AINTI  (FSELTC(I)*OS)  /OPirS(PTR)  )+AINT((Fr.F,'HC(I)*(I-OS)  ) 
4  /OPHRS (PTR) ) 

TSECC(  I)*TSECC(I)  +  <’n,"v*C\0B(PT!,)+Fnn(p7p)4.crqr,-ir  (  t  ^  / 
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•-  .  1 


510  S,  OPHRS(PTR)  >*Co"«Bn’T:'’)*Pr':M 

520  20  continue 

530C 

5'40C  NOW  THE  TOTAL  COST  FOR  E  \CM  ITT"  * 

550*  OETFP-'I'Tn  (TDCST,T3CST,AND  TCCST) .  FINALLY  USING  A  * 

560*  SERIES  OF  IF  STATEMENTS  A  MARL  l*.' At.  ANALYSIS  IS  * 

570*  PERFORMED  TO  IDENTIFY  IF  ONE  REP\tp  LEVEL  IS  Vo?r  * 

530*  ECONOMICAL  T’LYN  ANOTHER.  IF  THIS  IS  TRUE  A  MNIOYE  * 

590*  I  NO  I  CATER  IS  SET  ( INTLO  ,  INOf? ,  INOS  ,  IN'DO)  .  * 


600C 

610  IF(TSECO(I) .EO.fi)  CO  TO  33 

620  SF.O=TSr.CD(  I )  /M 

630  IF(TSECC( I ) .EO.O)  CO  TO  33 

640  SEC-TSF.CC(I)/M 

650  33  TDCST(I)-OEPOLC(I)+SED 

660  TBC5T  ( I )  -  5ASELC  ( I )  +TSF.C  0(1) 

670  IF(CIRF.EQ.O)  TCCST(I ) -JUMBO 

680  TCCST ( I )-CIRFLC ( T  )+S  F.C 

6DQ  TF(LP"'T)(I).CT.O.OR.LFMOS(T)  .CT.O)  CO  TO  80 

700  TF(SCPPLC(I) .LT.DEPOLO(I))  TNDD-1000 

710  IF(TDCST (I) .LT.SC1PLC (I ) )  INDS-l^O 

720  80  IF(LFMOS(I). CT.O. OR. LFMOB(I). CT.O)  CO  t o 

730  IF(T3CST(I) .LT.SCEPLC(I) )  I“D3«loo 

740  TF(SCRPLC(T)  .LT.BASELC(D)  INDR-IO 

750  32  IF  (LF'OB(I) .GT.O.OR.LFMnO(I) .CT.O)  CO  TO  34 

760  IF(T3CST(I) .LT.DEPOLC(T))  INDD-IORO 

770  I F ( TDCST ( I )  . LT . B AS F.LC ( I ) )  INDB-IO 

780  34  IF(LFMDC(I). CT.O. OR. LFMOD(I). CT.O)  CO  TO  74 

790  IF (TCCST ( I ) . LT.DEPOLC ( I ) )  INDD-IOOO 

800  IF(TOCST(I)  .LT .CtP.FLC(I) )  INDC-1 

810  74  IF(LFMOC(I).Gr.C.0R.LF‘,O3(l).CT.'))  CO  TO  77 

820  IF(SCRPLC(I) .LT.CIRFLC(I) )  INDC-l 

330  IF(TCCST(I)  .LT.SCP.PLC(I) )  IVPS-100 

840  77  IF(LF^OCd). CT.O. OR. LFMOR(I). CT.O)  GO  TO  37 

350  IF(TCCST(I) .LT . 8ASFLC ( I ) )  INOR-10 

860  IF(TBCSTd)  .LT.CIRFLC(I))  INOC-l 

870C 


330C  LINES  BELOW  CHECK  THE  VALDES  OF  THE  INDICATERS  * 
390*  IF  THE  I GDI CATER  HAS  BEEN  GIVEN  THE  RIGHT  VALUE,  THE  * 
900*  NRLA  FAILURE  MODE  EXCLUSION  ARRAY  IS  UPDATED  TO  * 
910*  ELIMINATE  THAT  REPAIR  OPTION  FROM  FURTHER  CONSIDER-  * 
920*  AT ION  FOR  THE  Ith  ITEM.  * 


930C 

940  37  IF(LFMOD(I)  .F.Q.2)  CO  TO  31 

950  I F ( I NDD . EO . 1 O00 )  LF*OD(I)-l 

960  31  IF(LFMOSd)  .E0.2)  CO  TO  3  3 

970  IFdNDS.EQ.  100  )  LFMOS(I)-l 

930  33  IF(LFMORd)  .E0.2)  CO  TO  35 

990  IFdNDB.EQ.lD)  LFM03(I)-l 

IO00  c 5  I F ( L FMOC (D.GF.2)  CO  TO  38 
1010  IFdNDC.ro. 1)  I.F’OC (I ) - 1 

1020  °S  TYP,JPR-INnn+INDS+I'TDB+TNDC 
1030C 


based  on  nr.  vat.'t  of  rv’iTK  r»:K  srArr"F'-n  * 

CHECK  TO  SEE  IF  3  JUT  OF  THE  4  POSSIBLE  OPTION'S  HAVE  * 

SEEN  r.LIMINATr.O.  IF  THIS  IS  TRUE  THE  OPTIMAL  DECISION* 

HAS  VEEN  MADE  AND  IS  STORF.D  IN  THE  ARRAY  'OPDECI/  * 

I F ( TYPRPR • EO . 0 1 1 1  )  OPDF.CL ( I ) -LOCAT ( 1 ) 

IF(TYPRPR.FO . 101 1  )  OPDECL(I)=LOCAT(2) 

IF(TYPRPR.FO.l  101 )  OPDECL ( I  )=L0CAT  ( 3 ) 

IF( TYPRPR.EQ.il 10)  OPD ECL ( T ) =L0C AT ( 4 ) 

TSEC0(I)=o 

TSEC3 ( I  )=0 
TEECC( T ) =0 

DO  LOOP  25  UPDATES  THE  *TUMPFR  OF  PIECES  OF  SE  AT  EACH  * 
LEVEL  EASED  ON  THE  REPAIR  LEVEL  DECISIONS  THAT  HAVE  * 
T.EEN  ’LADE  AS  CELL  AS  THE  SF.  HOURS  STILL  AVAILABLE.  * 

IT  THEN  COES  BACK  AND  CHECKS  TO  ‘.LAKE  SURE  THAT  ONLY  * 

ONE  ADDITIONAL  PIECE  OF  SE  IS  SUFFICIENT.  * 

DO  25  125=1  ,26 
PTR=SEARY (1,125) 

IF(PTR.EQ.O)  CO  TO  25 

HRAHSE(PTR)=NSECI (PTP)*(OPHRS (PTR)-RSYHRS  f PTR) ) 
ID-SECODF.(PTR)  / 1 0 
IF(LFMOD(T)  .C.T.O)  CO  TO  41 
IF(TYPRPR.EO. Dili. AND. TD.LT. 5000)  CO  TO  42 

41  IF(LFM03(I) .CT.O)  CO  TO  4S 

IF(TYPRPR.EO. 1101. AND. ID.CF. 5000. AND. TO. LT. 0000)  CO  TO  44 
4S  IF(LF>!OC (I )  .CT.O)  CO  TO  10 

IF(TYPr,PR.EO. 1110. AND. ID.r.E.QOOO)  CO  TO  46 
GO  TO  25 

42  IF(HRAVSE(PTR).0E.FSEUHD(I))  CO  TO  43 

REOMT(PTR)»REOMT(PTR)+l 
H RAVS  F ( PTR ) «H RAV5F.  ( PTR  )  +OPHRS  (PTR) 

GO  TO  42 

43  USEHR?(PTR)-USEHRS(PTR)+FSF.l:UD(I) 

FAC-FDB(PTR) 

IF(NSF.CI(PTR)  .CT.O)  FAC-0. 

TSECST  (PTR)  -TSF.CST  (PTR)  +(  (REQUT (PTR )  *CAD B ( PTR )+cAC+ 
f.  ((FSEUHD(I)/0?HRS  (PTR  )  )*OQOE  (  PTR)*?  TUP  )  )  /  PRO  AT  (  M)  ) 

HR AVS E (PTP ) -HRAVSE (PTR ) -^SEUID ( I ) 

CO  TO  25 

44  IF (HRAVSE (PTR)  .GE.FSE:’N3(I ) )  CO  Tn  45 

P.EO!  !T  ( PTR )  *R  EOMT  ( PTE )  + 1 

HRAVSE (PTR)-HPAVSF(PTR)+OPl,PS  (PTR) 

CO  TO  44 

45  USE HRS (PTR )- USE HRS (PTE)+FSEUHB ( T ) 

FAC-FDB(PTP) 

I F(NS  EC  I f nTR ) .CT.O)  *\\C«0 . 

If, FCST (PT9 ) -TSF.CST ( PTP )  +  (T,E0TT ( T’TR ) *C.\DP. ( n T • )  +FAO(  (’TE"I!r. 
6  /OP”RS(PTn) )*CODE(nTR)*?IUP) ) 

HRAVS  0  ( pTR  )  -HR  AVS  F.  (  PTR  )-FS  ECU  R  ( I ) 

CO  TO  25 


iWV.4 


15  70 
15S0 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 

16  300 
1690* 
1700* 
1 7 10C 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1330 
1840 
lp  50 
1360 
1870 
1880 
1390 
190^c 
191°* 
1920* 
1933* 
1940C 
1950 
I860 
1970 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2090 
7090 


44  IF(HRAVSE(PTR)  .GE.FSKUHG(  I)  )  G'i  TO  4  7 

RFAMT(PTR)-2+AINT (  ( FSPUHC ( I )*OS )  /OPHRS  (PTP.)  )  + 

&  AIN'T ( (  FSCl'HC  ( T )  *  ( l-OS ) )  /OPHRS  (PTP.) ) 

HRAVS  8  ( PTR )  =HRAVS  E  ( PIP )  +  ( P.EOMT  (  PT.R )  *OP HRS  ( PTR ) ) 

4  7  l'S  EH  P.3  (  PTR  ) =?  JS  EH  RS  ( PTP )  +FS  EU  HC  ( I ) 

FAC=FD3(PT.R) 

IF(NSECI(PTP.)  .GT.O)  FAC=0. 

T5ECST(PTR)=TSECST(PTR)+(  (?.EQMT(PTR)*CADB(PTR)+FAC+ 

&  (  (FSEUIIC ( I  ) /OPHRS  (PTR)  )  *CODR ( PTR)  *P  I'JP ) )  /  FLOAT!!1) ) 

II P  AVS  E  ( PTR )  =Ii  RAVS  F.  ( PTR )  -FS  E'JHC  ( I ) 

25  CONTINUE 

IF  THE  CURRENT  ITEM  IS  AN  LRU  FM  CO  TO  10  * 

BECAUSE  THE  FOLLOWING  STATEMENTS  RELATE  ONLY  TO  3R''S.  * 

00  TO  10 
15  CONTINUE 
J=I-LFMS 

IF(OPDECL  (LRUP  TR  ( J)  )  .F.Q  .LOCAT (15))  GO  TO  10 
IF(SOD(J).NE.l)  GO  TO  22 
SOD ( J) -2 
INDD-1000 

22  IF(SOB(J).NE.l)  GO  TO  24 

308( J)=2 
INDE-10 

24  IF(SOC(J).NE.l)  GO  TO  26 
SOC( J)=2 
INDC-l 

26  I F(S03 ( J ) .NE . 1 )  GO  TO  28 

SOS (J)=2 
INDS-100 

23  CONTINUE 

IF(CIRF.EQ.O)  S0C(I )«3 

DO  LOOP  20  CHECKS  FACM  POSITION  OF  THE  3E/LRU/SRU  * 

CROSS  REFERENCE  ARRAY  (SF.ARY)  FOR  THE  POINTER  VALUES  * 

OF  THE  SE  NECESSARY  TO  PEPAIP  THE  Ith  ITEM.  * 

DO  50  150*1 ,26 
PTR-S  EARY (1,150) 

TF(PTR.EO.O)  GO  TO  50 
ir*SECODF. (PTR)  / 1 0 

tF(SOD(J) .GT.0.OR.TD.GE.5000)  GO  TO  52 
3UY-AINT ( l .0+(SSEUHD (J) /ORHRS (PTR ) ) ) 
TSECD(J)»TSECD(J)+(BUY*CADB(PTR)+FDR(PTR)+(SSEUHD(J)/ 
f,  OPHRS  (PTR) )  *COOB  ( PTP )  *?T,TP) 

52  TP(SOB(J).GT .0  .OR  .ID.T.T . 500n .OR. ID .GE.9000)  GO  TO  54 
BUY- M NT (l .0+(SSFUHB( J) /OPUPS (PTR) ) ) 

TS  r.CR  ( J )  -TSECP  ( J  )  +  ( RUY* C  AD P' (  PTR )  +FD B ( ?TR  )  +  (  (SSF.U  IB  ( J  )  / 
o.  OPHRS  (nTR  ))*CODr,(  PTP. )*RIUP)) 

54  IF(SOC(J) .GT.O. OR.TD.LT. 9000)  GO  TO  50 

B"Y-2+MNT((SSEUHC( J )*OS) /OPHRS (PT” ) )+AINT( (FSF"HC(J)*( l-OS) ) 
&  /n?HRS(?TR)) 


tslcc(j)=tsecc(  j)+cm’y*c  \de ( p r* )  +"n s ( pt^ ) + (  n  / 

&  o?:i:<s(PTR))*co'i),5CPTK)*pir:j) 

50  CONTINUE 

SOU  THE  TOTAL  COST  FOR  RACE  ITEM  IS 
DF.TEP'!T\’F.P  (TOCST  ,TBCST  ,ANf)  TCCST).  FIS' ALLY  USING  A  * 

SERIFS  OF  IF  STATEMENTS  \  MARC I SAL  ANALYSIS  IS  * 

0 F.R FORMED  TO  IDENTIFY  IF  ONE  REPAIR  LEVEL  TS  wo?f  * 

ECONOMICAL  THAI'  ANOTHER .  IF  THIS  IS  TRUE  A  ULTO'JE  * 

INDICATE?  IS  SET  ( TMDr' ,  I  NOE  ,  I'PS  ,  INOC)  .  * 

IF(TSECD(J) . EO.O)  CO  TO  5S 
SED=TSr.CD(J)/M 
IF(TSF.CC(J)  .EQ.O)  CO  TO  58 
SEC*TSECC(.I)/M 
55  TOCST (.1  )  =  DCP03C(J)  rSF.D 

IF (OPDECL(LRuPTP. ( J ) )  •  EQ.LOCAT  ( 3 ) )  TOCST  ( J  )=TPCST  (J)  + 

&  BRLSS  (  J )+RRLr>S  ( J ) 

TOCST  (  I )  *  ?>  AS  ESC  (  J  )  +TS  EC  8  (  J  ) 

TCCST  ( J 1  *CIRFSC  (.D+SF.C 

I F  (OPPECL (LRUPTR ( J ) ) .EO.LOCAT ( A ) )  TCCST (J )»TCCST (J )+?RLCS (J ) 

I  F(OPr>ECL(I.R:J?TR(  J )  )  . EQ.LOCAT (3)  )  SCRPSC  ( J  )=SCRPSC (  J )+BRLSS  (.1) 
IF(CtRF.F.O.O)  TCCST (.I)=.IUMRO 

IF(S 00 (J).CT.O. OR. SOS (J).CT.O)  CO  TO  50 
IF(SCRPSC(J) .LT.DEPOSC(J))  IM0D-1000 
IF(THCST(J) .LT.SCRPSC(J) )  INOS-100 
50  IF(SOS(J).OT.O.OR.SOB(J).GT.O)  CO  TO  52 
IF(TRCST(J) .LT.SCRPSC(J))  INDS-IOQ 
IF  (SCRPSC  (. I)  .LT  .R.ASESC  (  J) )  INDB-10 
62  IF(SOB(.J)  .CT.O. OR. SOD(J)  .CT.O)  CO  TO  Si 
IF(T5C3T ( J) .LT . DEPOSC ( J) )  INDD=1000 
I F (TOCST (J ) .LT . LAS  ESC ( J ) )  IND3=10 
6i  I  F(SOC.(  J)  .CT.O.OR.SOD(J)  .CT.O)  CO  TO  57 
IF(TCCST ( J) .LT.nE?OSC(J) )  TNDD-IOOQ 
I F(TOCST(.I )  .LT.CIRFSC(J))  INOC-I 

67  IF(SOC(J) .CT.O. OR. SOS(J) .CT.O)  CO  TO  SQ 

IF(SCRPLC( J) .LT.CIPFSC(J))  IMDC-l 
IF(TCCST(J)  .LT.SCRPSC(J))  INOS=lOO 

68  LF(SOC(J) .CT.O. OR. S03(J) .CT.O)  GO  TO  59 

I  "(TCCST ( J I  .I.T  ."ASESC(.J)  )  IN0P.=  10 
TF(TRCST(.J)  .LT.CTrlFSC(J) )  IN1C=l 
Ir(0?r)F.CL(LRUPTR(J) ) . En .LOCAT (  J ) )  SCRPSC  (J ) -SCRPSC (J ) 

5  -OPLSS(J) 

LINES  PF.f.ov;  CHECK  THE  VALLES  OF  THE  INOICATERS  * 

IF  THE  INDICATED  HAS  REEN  CTVEN  THE  RIGHT  ,rM.UE,  THE  * 

MELA  FAILURE  NODE  EXCLUSION  ARRAY  IS  UPDATED  TO  * 

ELIMINATE  THAT  REPAIR  OPTION  FROM  FARTHER  CONSIDER-  * 

AT  ION  FOR  THE  Itn  ITEM.  * 

TF( SOD ( J ) . CO . 2 )  CO  TO  61 
IFdNDD.EO.lOOO)  SOD(.l)-l 
■1  I F ( JOS (  I ) . F0.2 )  GO  TO  53 


no 


2630  ir(ms.ro.inn  >  s-'S(j)  =  i 

2640  63  IF(S0D(  J)  .Fn.2)  GO  TO  66 

2650  I F  ( INDR,  .GO.  10)  SOR(j)«l 

2660  65  IF(S0C(.I)  .CE.2)  CO  TO  66 

2670  IF(INDC.EO.l)  SOC(J)=l 

2630  66  TVPRPR=r:Dn+r;js+i':np+rnr 

26900 

2700*  3  AS  TO  ON  TUP  VALUF.  OF  TtPRPR  THE  SELDU  STATEMENTS  * 

2710*  CHECK  TO  SEE  IF  3  DMT  OF  THE  4  POSSIBLE  OPTIONS  HAVE  * 

2720*  SEEN  ELIMINATED.  IF  THIS  IS  TRUE  TilE  OPTIMAL  DECISION* 

2730*  HAS  SEEN  MADE  AMD  IS  ST07FD  IM  THE  ARRAY  'OPDECS'  * 

2740C 

2750  IF(TYPRPR.F.O.OUl)  OPDFCS  ( J)=l.OC  NT  ( 1 ) 

2750  IFfTYPRPR.EO. 1°1 I )  OPDECS(J)*LOCAT(2) 

2770  IF(TYPRPR.EO.llOl)  OPDUCS(J)=LOC\T(3) 

2 7 SO  IF(TYPRPR.EQ.IIIO)  OPRECS ( J )-LOCAT(4) 

2790  TSECD(J)=o 

2800  TSF.CB  (J)=0 

2810  TSECC( J)=0 

2820C 


2S30*  DO  LOO?  75  UPDATES  THE  MUMPER  OF  PIECES  OF  SE  AT  SAC"  * 

2840*  LEVEL  BASED  ON  THE  REPAIR  LEVEL  DECISIONS  THAT  HAVE  * 

2850*  SEEM  MADE  AS  '.'ELL  \S  THE  SE  u.yjps  STILL  AVAILABLE.  * 

2860*  IT  THEN  GOES  BACK  AMn  CHECKS  TO  MAKE  SURE  THAT  ONLY  * 

2870*  ONE  ADDITIONAL  PIECE  OF  SE  IS  SUFFICIENT.  * 

2880C 

2890  DO  75  175*1,26 

2900  P7R*S  E  ARY ( T , 1 7  5 ) 

2910  IF(PTR.EO.O)  CO  TO  75 

2920  HRAYSE(PTR) *NSECI (PTE) * (OPHRS  (PTR)  -E.SVHRS  ( PTR ) ) 

2930  ID*»SECODE(?TP.)  / 1 0 

2940  IF(SOD(J).CT.O)  CO  TO  91 

2950  IF(TYPRPR. EQ.OIII.AND.ID.LT. 5000)  CO  TO  92 

2960  91  IF(S03(J).CT.0)  CO  TO  98 

2970  IF(TY?RPR.EQ.l 101 . AND. ID. CE. 5000. VvD.ID.LT.9009)  CO  TO  94 

2980  98  IF(SOC(.J)  .GT.O)  GO  TO  LO 

2990  IF(TYPRPR.EO. 1110. AND. ID. CE. 9000)  DO  TO  93 

3000  GO  TO  75 

3010  92  IF(HRAVSE. (PTR)  .CE.SSEUHD(J) )  GO  TO  93 

3020  REQ.MT  (  PTR )  *R  EQMT  ( PT  P ) +1 

3030  HP.  A  VS  F.  (PTR )  -HP.AVS  F.  ( PTR  )+*)PHRS  ( PTR  ) 

3040  CO  TO  92 

3050  93  USEHP.S  (PTR) “USEHRS  ( PTR)+SSEUHD (J ) 

3050  FAC-FDB(PTP) 

3070  I F (NSEC I (PTR) .CT.O)  FAC=0. 

3030  TSECST  (PTR  )“TSF.CST  (PTP.)  +  (  (REOv'T(FTR )  *CADP.  (PTR  )+FAC+ 

3090  &  ( (SSEUHD(J) /OPHRS (PTR ))*COD3 (PTR) *PIUP)) /FLOAT (M) ) 

3 1  DO  HRAVS  E  (  PTR )  -IPAVS  F.  (  PTR)  -SS  EU"D  (  J ) 

3110  CO  TO  75 

3120  94  I F (wraVSF. (PTR )  .CE.SSFUHB( J) )  GO  TO  95 

3130  REOMT(PTR )*9EQvir(°TR )+l 

3140  H  R  A  V  f>  E  (  P  T  R )  =  H  F.  A  V  S  E  (  P  T  R  )  +0  ?  M  R  S  (  P  T  R  ) 

3150  CO  TO  94 

3160  95  USEHRS  ( PTR) “USEHRS  (?TC )+S3EUHB( J) 

3170  FAC*FD3  ( PTR ) 


r vs i:c r  ( r tr ) . ft. o )  r aoo . 

TSECS r ( PTR ) -TSIICST  ( PTR )  +  (  ”.KO>! r ( PTR )  *CAr)3(PT!>. )  +F.V> ( ( 
&  /OPHRS (PTR )  )*COD!5(PTP. ) *?I"P ) ) 

HRAV R  R  (  PT R )  =:•■  RAV R  F  (  PTR  )  -5  5 F.l ’ : !  3  ( .!  ) 
no  TO  75 

96  IF(PRAVSF.(PrR)  .OF,  .SSF.l'PC  (  J  )  )  GO  TO  97 

RRQMT ( PTR )  -2+AINT ( (SSH'JHC(J) *OS )  /OP"  ’ S  (PTP  >  )  + 

S.  AINT  (  (SSFUHC  ( J )  *  ( L-OR ) )  /OPP.RS  ( PTP ) ) 

HRAVSE (PTR ) -"RAVST ( PTR)+(PEqMT ( PTR)*0?MP5 (PT& ) ) 

9  7  ITS E'!?.S  ( PTR )  «US F.HRS  ( PTR >  +S S E'.’HC  (  J ) 

FAOmr.(PTf:) 

irccr.ECKPTR) .o-t.o)  r\c=o. 

TS KCST  ( °TR )  -TSTCS T  (  PTR )  +  (  ( RR O’-rr  ( PTR )  *CAO 3  ( PTR )  +F AC+ 

1  ((SRF.UTC(J)  /0?P  AS(PTP)  )*COOF  (nTP.)*PIUP)  )/ FLOAT (M)  ) 
i RAV R r. ( PTR) RAV RE  (  PTP  )  -SSEL'HO  (  J  ) 

75  CONTINUE 
IND3-0 
IVJC-0 
IN99-0 
[T9S=n 
in  rn,:'rT‘"'f 


Appendix  E:  MSETNT  Subroutine 


10  SUBROUTINE  MSETNT (M5GL2 ) 

20C 

30C 

40C  CODE  FROM  STATEMENT  •«  2000  THROUGH  It  2210  UTILIZES  THE  ABOVE 

50C  COMPUTED  ITEM  &  SE  COSTS  TO  PRODUCE  AN  RLA  COST  NETWORK. 

60C 

70C  AS  A  PRELIMINARY  STEP  DETERMINE  NODE  NUMBERS  FOP,  KEY  NODES  IN 

80C  T!'E  NETWORK.  THE  VARIABLES  ARE  CODEO  WITH  THE  PREFIXES  'LD' 

90C  FOR  'LAST  DEPOT'  AND  'LB'  FOR  'LAST  BASE'. 

100C 

110  2000  LDSE-NDSE+1 
120  LDLRU-LDSE+LFMS 

130  LDSPU-LDLRU+SP’S 

140  LBLRU-LDSPU+LFMS 

150  LBSP.U-LBLRU+SFMS 

160  LBSE-LBSRU+NBSE 

170  LNODE-LBSE+l 

180  I F(LNODE .LE .MAXNOD )  CO  TO  2005 

190  WRITE(6 , 2002 ) 

200  2002  FORMAT ( ' 0  MODE  VECTORS  ARE  TOO  SMALL  --  STOP'’) 

210  STOP 

220C 

230C 

240C  SET  ft  OF  ARCS  TO  ZERO  &  It  OF  NODES  TO  1 

250  2005  NARCS-0 
260  NNODF.S-l 

270C 

280  IF(NDSE.LT. 1)  CO  TO  2015 

290C  CREATE  ARCS  FROM  SOURCE  TO  DEPOT  SE  —  DEPOT  SE  COST  ARCS 

300  DO  2010  12010*1 ,NDSE 

310  NARCS-NARCS+1 

320  NNODES-NNODES+1 

330  SRCE(NARCS)-! 

340  DEST (NARCS )-NNODES 

350  CAP (NARCS ) -S ECOST (12010) 

360  2010  CONTINUE 
370C 

380C  CREATE  ARCS  FROM  THE  SOURCE  TO  THE  DEPOT  LRU  FAILURE  MODE 

3R0C  MODES  —  ITEM  RELATED  COSTS  FOR  LRU  REPAIR  AT  DEPOT 

400  2015  DO  2020  12020-1, LFMS 

410  NARCS-NARCS+1 

420  NNODES-NNODES+l 

430  SRCE(NARCS)-1 

440  DEST (NARCS )-NN0DES 

450  CAP (NARCS ) -DEDOLC (12020) 

460  2020  CONTINUE 
470C 

4«0C  CREATE  ARCS  FROM  THE  SOURCE  TO  THE  DEPOT  SRU  NODES  —  ITEM 

490 C  RFEATED  COSTS  FOR  S“U  REPAIR  AT  DEPOT 

500  TF(RFMS.LT.t)  CO  TO  2040 


510  no  2030  12030=1 ,stms 

520  NARCS-NARCS+l 

530  NNODES-NNODES+l 

540  SRCE (NARCS )  =  l 

550  DEST (NARCS ) -NNODES 

550  CAP  (N'ARCS  )  =T)EPOSC  (12030) 

570  2030  CONTINUE 
580C 

590C  CREATE  ARCS  FROM  THE  DEpOT  SF.  NODES  TO  THE  DEPOT  LR1J/SRU  NODES 

600  2040  IF(NDSE.LT.l)  GO  TO  2070 
610C  FOR  EACH  DEPOT  SE  — 

620  DO  2060  12060=1, NOSE 

630C  SEPF  CONTAINS  A  POINTER  TO  THE  SE  CROSS  REFERENCE  (SEXREF) 

640  IPTR=SEPF(T2060) 

650  IF(IPTR.EQ.O)  CO  TO  2060 

660C  SEXREF  CONTAINS  THE  9  OF  AM  LRU  FAILURE  MODE  (ITEM  >  0)  OR 

670C  THE  NUMBER  OF  AN  SRU  (ITEM  <  0) 

680  2050  ITEM=SEXREF ( IPTR) 

690  IF(ITEM.LT.O)  ITEM-LFMS-ITEM 

700  NARCS =NARCS+1 

710  SRCE (NARCS) =12060+1 

720C  ITEM  IS  A  POINTER  TO  AN  LRU  OR  SRU  AND  CHEN  ADDED  TO  LOSE 

730C  GIVES  THE  APPROPRIATE  NODE  NUMBER 

740  DEST (NARCS )=LDSE+ITEM 

750  CAP (NARCS) -JUMBO 

760C 

770C  NXTITM  IS  A  POINTER  TO  THE  NEXT  ENTRY  IN  SEXREF  WHICH  IS  AN 

780C  LRU  OR  SRU  REQUIRING  THE  CURRENT  SE,  I.E.,  SE  #  12050. 

790  NPTR-NXTITM(IPTR) 

300  IF(NPTR.EO.O)  GO  TO  2060 

810  IPTR-NPTR 

820  GO  TO  2050 

830  2060  CONTINUE 

340C 

850C  CREATE  THE  ARCS  EMANATING  FROM  THE  DEPOT  LRU  FAILURE  MODE  NODES 

860  2070  DO  2030  12080=1, LFMS 
870C  CHECK  FOR  AN  SRU 

880  IPTR*SRUPTR(T2080) 

890  IF(IPTR.F.O.O)  GO  TO  2075 

900C  CREATE  AN  ARC  TO  THE  DEPOT  SRU  NODE  —  COSTS  UNIQUE  TO  BASE 

910C  REPAIR  OF  LRU  &  DEPOT  REPAIR  OF  SRU 

920  NARCS-NARCS+l 

930  SRCF.(NARCS)=LDSF.+I2080 

940  DEST (NARCS )«LQLRU+TPTR 

950  CAP (NARCS ) -BRLDS (IPTR) 

960C  CREATE  AN  ARC  TO  THE  BASE  LRU  NODES  —  COST  OF  SCRAPPING  THE  LRU 

970  2075  NARCS-NARCS+l 

980  S  PCE  ( NARCS  )=LDSF.+1 2080 

9°0  DEST (NARCS ) -LDSRU+I208O 

1000  CAP (NARCS ) -SCRPLC (12080) 

1010  2080  CONTINUE 
1020C 


114 


1  n 30 C  CREATE  ARCS  FROM  THE  DEPOT  SP"  NODES  TO  THE  PAS F.  SnU  NODES  — 

1040C  COST  OF  SCRAPPING  THE  SRC 

1050  IF(SFMS . LT. 1 )  GO  TO  2100 

1060  DO  2090  12090=1, SFMS 

1070  NARCS-NARCS+l 

1OS0  SRCE (NARCS )=LDLRU+I 2090 

1090  DEST(NARCS)=LBLRU+I2090 

1100  CAP  (N'ARCS )  =SCRPSC  ( 12090 ) 

1110  2090  CONTINUE 
1120C 

1130C  CREATE  ARCS  EMANATING  FROM  THE  RASE  LRU  NODES 

1140  2100  DO  2120  12120=1, LFMS 
1150C  CHECK  FOR  AN  SRU 

1160  IPTR-SRUPTR (12120) 

1170  IF(IPTR.EO.O)  CO  TO  2110 

1180C 

1190C  CREATE  THE  ARC  FROM  THERASE  LRU  NODE  TO  THE  BASE  SRU  NODE  ~ 

1200C  COSTS  UNIQUE  TO  BASE  REPAIR  OF  THE  LRU  AND  SCRAPPING  (OR  DEPOT 

1210C  REPAIR)  OF  THE  SRU 

1220  NARCS *MARCS+1 

1230  SRCE (NARCS ) =LDSRU+I2 1 20 

1240  DEST (NARCS ) -LBLRU+IPTR 

1250  CAP (NARCS )=BRLSS (IPTR) 

1260C  CREATE  AN  ARC  FROM  THE  BASE  LRU  NODE  TO  THE  LAST  NODE  (SINK)  — 

1270C  ITEM  RELATED  COST  OF  RASE  REPAIR  OF  THE  LRU 

1280  2110  MARCS »NAP.CS+1 

1290  SRCE(NARCS)=LDSRU+I2120 

1300  DEST (NARCS )=LNODE 

1310  CAP(NARCS)=BASELC(I2120) 

1320  2120  CONTINUE 
1330C 

1340C  CREATE  ARCS  EMANATING  FROM  THE  BASE  SP.U  NODES 

1350  IF(SFMS.LT.l)  GO  TO  2140 

1360  DO  2130  121 30=1, SFMS 

1370C  CREATE  AN  ARC  FPOM  THE  BASE  SRU  NODE  TO  THE  BASE  LRU  NODE  — 

138'):  THIS  ARC  PREVENTS  A  DECISION  TO  DO  BASE  LEVEL  SRU  REPAIR  UNLESS 

1390C  THE  LRU  IS  ALSO  BASE  REPAIRED 

1400  NARCS-NARCS+1 

1410  SRCE  (NARCS  )=LBT.RU+T2 130 

1420  DEST  (NARCS )  =LDS  P.U+LPUPTR  (12130) 

1430  CAP (NARCS) -JUMBO 

1440C 

1450C  CREATE  AN  ARC  FROM  THE  BASE  SRU  NODE  TO  THE  LAST  NODE  (SINK)  — 

1460C  ITEM  RELATED  COST  OF  BASE  REPAIR  FOR  THE  SRU 

1470  NARCS-NARCS+1 

1480  SRCE (NARCS )«LBLRU+I21 30 

1490  DEST (NARCS )=LNODE 

1500  CAP (NARCS )=BASESC(I2 130) 

1510  2130  CONTINUE 
1520C 

1530C  CREATE  ARCS  TO  AND  F°0"  THE  BASE  SE  NODES 

1540  2140  IFC1BSE.LT. 1)  GO  TO  2180 
155QC  FOP  EACH  BASF.  SF.  — 
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1560  NDSEPl-NDSE+l 

1570  NBSEPO-NDSF+NBSE 

1580  DO  2170  I2170-NDSEP1 .NBSEPO 

1590C  J2170  IS  A  POINTER  TO  THE  BASF  SF 

1600  J2170-I21  70-NDS  F. 

1610C  SFPF  CONTAINS  ’A  POINTER  TO  THE  SF  CROSS  REFERENCE  (SEXREF) 

1620C  (  SEE  COMMENTS  FOR  STATEMENTS  2040  TO  2060) 

1630  IPTR*SEPF( 12170) 

1640  IF(IPTR.EQ.O)  CO  TO  2160 

1650  2150  ITEM=SEXREF(IPTR) 

1660  IF( ITEM.LT .0)  ITFM=LFMS-ITEM 

1670  NARCS-N'ARCS+l 

1680  SRCE (NARCS ) =LDSRU+ITEM 

1690  DEST (NARCS )=LBSRU+J2 170 

1700  CAP (NARCS) -JUMBO 

1710C 

1720  NPTR-NXTITM(TPTR) 

1730  IF(NPTR.EO.O)  CO  TO  2160 

1740  IPTR-NPTR 

1750  GO  TO  2150 

1760C 

1770C  CREATE  AN  ARC  FROM  THE  BASE  SE  TO  THE  LAST  NODE  (SINK)  — 

1780C  BASE  SE  COST 

1790  2160  NARCS-NARCS+1 

1800  SRCE(NARCS )»L8SRU+J2 1 70 

1810  DEST (NARCS )-LNODE 

1820  CAP ( NARCS )  -S  ECOS  T(I2170) 

1830  2170  CONTINUE 

1840  2180  NNODES-LNODE 

1850C 

1860C  EXPLICITLY  SORT  THE  SRCF-DFST-CAP  TRIPLES  INTO  ASCENDING  ORDER 

1870C  BY  SRCE  (FLOW  IS  USED  BY  SORT  TO  STORE  ORDERING  POINTERS) 

1880  CALL  SORT(SRCE, DEST, CAP, FLOW, NARCS, MAXARC , 1 ) 

1890C  SAVE  POINTERS  FOR  FORWARD  SCAN  IN  FVnSP 

1900  DO  2190  12190*1 .NARCS 

1910C  UTILIZE  THE  VALUES  IN  SRCE  IN  REVERSE  ORDER.  (IF  NODE  '  V  IS  A 

1920C  SOURCE  FOR  'N'  ARCS  THEN  FWDSP(l)  WILL  SUCCESSIVELY  GET  THE 

1930C  VALUES:  N,  N-l ,  K-2 ,  ....  2,  1.  SIMILARLY,  FWDSP(K)  WILL 

1940C  HAVE  THE  VALUE  'J'  WHERE  J  CORRESPONDS  TO  THE  FIRST  OCCURENCE 

1950C  OF  NODE  K  IN  SRCE.) 

I960  J2190-NARCS+1-I2190 

1970  K2190*SRCE( J2190) 

1930  FWDSP (K2 1 90)»J2 190 

1990  2190  CONTINUE 

2000  FWDSP(LNODE)-NARCS+l 

2010  LNODM1-LNODE-1 

2020  DO  2195  12 195-1 .LHODMl 

2030  J2 1 95-LNODE-I2 195 

2040  IF(FWnSP(J2195)  .EO.O)  F'-DSP(J2195)  =  FWDSP(J21R5+1 ) 

2050  2195  CONTINUE 
2060C 

2970C  USE  SORT  TO  IMPLICITLY  SF.PUENCE  THE  DF.ST  ENTRIES  INTO  ASCENDING 

2080C  ORDER.  THE  ORDER  WILI.  BE  SPECIFIED  BY  THE  ENTRIES  TV  Ri’PTR. 


2O90C  THUS,  THE  ENTRIES  IN'  DACKSP  ARE  POINTERS  T«*  SRCE-DEST-GAr 

2101C  TRIPLES  RASED  ON  TMF.  ^EST  VALDES. 

2110  CALL  SORT(DEST,0,0,RKPTR,N\RCS,V4\ARC,0) 

2120  3ACKSP( 1  )-0 

2130  DO  2200  12200-1 .NARCS 

2140  J2200-BKPTR(T2200) 

2150  K2 200-DEST( J 2 200) 

2160C  BACKS  pa:)  IS  THE  POSITION  IN  3XPTR  WHICH  HOLDS  THE  LAST  POINTER 

2170C  TO  A  SRCE-DEST-CAP  TRIPLE  HAVINC  NODE  It  K  AS  THE  OEST.  THE  FIRST 

2180C  POINTER  TO  A  TRIPLE  HAVINC  NODE  ‘t  K  AS  THE  DEST  IS  BKPTR(K-1 )+l . 

2190  BACKS P(E2 20 0 ) *12200 

2200  2200  CONTINUE 


2210  DO  2222  12222-2 .LNODE 

2220  TF(BACKSP( 12222). EO.O)  BACKSP(I2222)= 

2230  2222  CONTINUE 

2240C 

2250C  IF(MSCL2.LT.1)  Co  TO  2205 

2260C  PRINT  2202, 'SRCE  ' , (SRCE(I) ,T-l ,NARCS) 

2270C  2202  FORMAT ( '  *,A6,('  ',2013)) 

2280C  PRINT  2202, 'DEST  ' , (DEST(I)  ,1-1  .NA’RCS) 
2290C  PRINT  2202, 'CAP  ' , (CAP(I) ,1=1 ,NARCS) 
2300C  PRINT  2202, 'FLOW  ' , (FLOU( I ) ,1-1 ,NARCS) 
2310C  PRINT  2202 , ' RKPTR  ' , (BKPTR ( I ) , 1-1 .MARCS ) 
2320C  PRINT  2202 ,' FWDSP  ' , (FWDSP(I) ,1-1 , LNODE) 
2330C  PRINT  2202 ,' BACKSP' , (RACKSP( I) , 1=1 .LNODE) 
2340  2205  DO  2400  12400-1 .NARCS 
2350  J-SRCE(I2400) 

2360  K=DEST(I2400) 

2370  IF(J.NE.l)  CO  TO  2230 

2380  IF(K.GT.LDSE)  CO  TO  2210 

2390  SF.ARCP  (K-l)-I2400 

2400  GO  TO  2300 

2410  2210  IF(K.GT.LDLRU)  GO  TO  2220 

2420  LDARC(K-LDSE) =12400 

2430  GO  TO  2300 

2440  2220  SDARC(K-LDLRU)-I2400 

2450  CO  TO  2300 

2460  2230  IF( J.LE.LDSE)  GO  TO  2300 

2470  IF(J.GT.LOLRU)  GO  TO  2250 

2480  IF(K.LE.LDSRU)  GO  TO  2240 

2490  LSARC(J-LDSE)«I2400 

2500  GO  TO  2300 

2510  2240  SBDARC(K-LDLRU) =12400 

2520  GO  TO  2300 

2530  2250  TF( J .GT.LDSRU)  GO  TO  2260 

2540  SSARC( J-LDLRU)-I2400 

2550  GO  TO  2300 

2560  2260  IF(J.GT.LBLRU)  GO  TO  2280 

2570  IF(K.NE. LNODE)  GO  TO  2270 

2580  LBARC  (J-LDSP.U )  -12400 

2590  GO  TO  2300 

2600  2270  IF(K.GT.LBSRU)  GO  TO  2300 

2610  SBSARC(K-LBLRU)-T2400 


BACKS P ( 12222 ) -PACKS P ( 12222-1) 


2350 

2360 

2370 

2380 

2390 

2400 

2410  2210 
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2620 

CO  TO  2300 

2630 

2280  IF(J.CT.LBSPU) 

CO 

TO 

2290 

2640 

IF(K.NE.LNODE) 

CO 

TO 

2300 

2650 

SBARC( J-LBLRU)- 

12400 

2660 

CO  TO  2300 

2670  2200  SEARCP(J+NDSE-LBSRU)-I2409 
2650  2300  FI/)W( 1 2400 )-0 
2690  2400  CONTINUE 
2700C 

27 IOC  THIS  ROUTINE  USEO  ONLY  IF  MARLA  IS  CALLED  IK  MAIN 

2720C 

2730  IF(MAR.EQ.O)  CO  TO  3000 

2740  DO  2450  12450-1 .MAXNOD 

2750  FWDSP(I2450)-o 

2760  RACKSP( 12450 )-0 

2770  2450  CONTINUE 

2780  DO  2460  12460-1 .MAXARC 

2790  BKPTR(I2460)-0 

2800  2460  CONTINUE 

2810C 

2820  DO  2181  12181=1, LFMS 

2830  I F (OPDECL (12181) .  EO . LOCAT (15))  CO  TO  2181 

2840  IPTR-LDARC(I218l ) 

2850  SRCE ( IPTR)-LN0DE+10 

2860  IPTR-LBARC(I218 1 ) 

2870  SRCE(IPTR)-LN0DE+10 

2880  IPTR-LSARC (12131) 

2890  SRCE ( IPTR)«LNODE+10 

2900  IF(OPDECL(I2 1 8 1 ) .NE.L0CAT(3) )  GO  TO  2191 

2910  CPTR-SRUPTR(I2181) 

2920  TC1-3RLDS  (CPTP. ) 

2930  TC2»BRLSS(CPTR) 

2940  L’=SDARC(CPTR) 

2950  C4P(W)=CAP(W)+TC1+TC2 

2960  X-SSARC(CPTR) 

2970  C AP ( X ) -CAP ( X ) +TC2 

2930  IPTPv-S  BS ARC  (CPTP ) 

2990  SRCE ( IPTR)*LM0DE+10 

3000  I PTR-S  BDARC ( CPTR ) 

3010  SRCE (IPTR)-LNODE+10 

3020  2191  DO  2182  12182-1 ,NARCS 

3030  I F(DEST (12182)  .EQ.  (LDSE+I2 18 1 ) )  SRCF.(I2l  82)-LNODE+10 

3040  IF (SRCE (12182 ).E0. (LDSE+I218 1 ) )  SRCE(T2182)=LN0DF>10 

3050  IF(SRCE(T2lS2).EO.(LDSRU+I218l))  SRCE(I2182)-LN0DE+10 


3050  IF(DEST(I2182)  .EO.  (LDSRU+I2 181 ) )  SRCF.(I2182)-LNODE+lO 

3070  2182  CONTINUE 

3080  2181  CONTINUE 

3090  DO  2184  12184-1, SFMS 

3100  IF(OPDECS (12184) .EO. LOCAT (15))  CO  TO  2134 

3110  I PTR-S D ARC (12184) 

3120  SRCE(  IPT?>  )-LNODE+10 

3130  I PTR-S E ARC (12184) 

3140  SRCF.(  IPTR  )-LNODE+l  0 


=w 


"k 


Xi2 


3150 

3160 

3170 

3180 

3190 

3200 

3210  2185 

3220  2184 

3230C 

3240C 

3250 

3260 

3270 

3280 

3290  2487 

3300 

3310C 

3320 

3330C 

3340C 

3350C 

3360C 

3370C 

3380 

3390 

3400 

3410  2490 

3420 

3430 

3440 

3450 

3460 

3470  2495 

3480C 

3490C 

3500C 

3510C 

3520C 

3530 

3540 

3550 

3560 

3570 

3580C 

3590C 

3600C 

3610 

3620  2500 

3630 

3640 

3650  2525 

3650C 

3670 


I  PTR=SSARC  (121-84) 

S8CF(IPTR)-I.V00c>in 

no  2185  12185=1 .MARCS 

IF(DEST(I2 1 85 ) .  F!Q .  (LDLRU+I2 l 84  )  )  5RCE(I2 1  85  )=L\TODE+lO 
IF(RRCE(I2185).EQ. (LRLRU+I2 134 ) )  SRCE ( T2 1 85 ) -LMODE+1 0 
I  F(0EST (12135)  .  EO .  (L3LRU+I2 184 ) )  SRC.G(I2185)=LNODE+10 
COM! I MU E 
CONTINUE 

EXPLICITLY  SORT  THE  SRCE-DEST-CAP  TRIPLES  INTO  AS  CENT  I  NO  ORDER 
BY  SRCE  (FLOW  IS  USED  BY  SORT  TO  STORE  ORDERING  POINTERS) 

CALL  SORT(SRCE,DF.ST, CAP  .FLOW, MARCS  .MAX ARC  ,  1 ) 

ITER=0 

PC  2487  12487=1 .NARCS 
IF(SRCE(I2487) .LT.LNODE)  ITER=ITER+1 
CONTINUE 
NARCS-ITER 

SAVE  POINTERS  FOR  FORWARD  SCAN  IM  FWOSP 
00  2490  12490=1 .MARCS 

UTILIZE  THE  VALUES  IM  SRCE  IN  REVERSE  ORDER.  (IF  NODE  ' 1 '  IS  A 
SOURCE  FOR  'N'  ARCS  THEN  FWOSP(l)  WILL  SUCCESSIVELY  GET  THE 
VALUES:  M,  N-l ,  N-2,  .  .  .  ,  2,  I.  SIMILARLY,  FWDSP(K)  WILL 
HAVE  THE  VALUE  *J*  WHERE  J  CORPUS PONOS  TO  THE  FIRST  OCCURENCE 
OF  NODE  *  K  IN  SRCE.) 

J2490=NARCS+l— 12490 
X2490=SRCE(J2490) 

FWDSP(K2490)«J2490 

CONTINUE 

FWOSP (LNODE)=NARCS+l 
LNODMl*LMODE-l 
DO  2495  12495=1, LNODM1 
J2495-LNODE-I2495 

IF(FWDSP(J2495)  .EO.O)  FWOSP ( J2495 )  =  FWDSP ( J2495+1 ) 

CONTINUE 


USE  SORT  TO  IMPLICITLY  SEQUENCE  THE  BEST  ENTRIES  INTO  ASCENDING 
ORDER.  THE  ORDER  WILL  BE  SPECIFIED  BY  THE  ENTRIES  IN  BKPTR. 

THUS,  THE  ENTRIES  IV  BACKS?  ARE  POINTERS  TO  SRCE-DF.ST-CAP 
TRIPLES  BASED  ON  THE  DEST  VALUES. 

CALL  SORT (DEST, O.O.BKPTR, NARCS ,MAXARC,0) 

RACKSP(1)=Q 

DO  2500  12500=1 .NARCS 

J2500=BKPTR( 12500) 

K2  50D=r'EST(  J2  500) 

BACXSP(K)  IS  THE  POSITION  IN  BKPTR  WHICH  HOLDS  THE  LAST  POINTER 
TO  A  SRCE-DEST-CAP  TRIPLE  "AVINC  NODE  #  K  AS  THE  DEST.  THE  FIRST 
POINTER  TO  A  TRIPLE  HAVING  NODE  ?  K  AS  THE  DEST  IS  BKPTR(K-l )+l . 
BACKS P(K2 500) =12 500 
CONTINUE 

DO  2525  12525=2, LNODE 

TF(BACKSP( 12525) .EO.O)  BACKSP( 12525 )=BACESP(T 2525-1 ) 

CONTINUE 

DO  2600  12600=1 .NARCS 


3630  J«SPCE(I2600) 

3690  K**f)F.ST  (12600 ) 

3700  IF(J.NE.t)  00  TO  2530 

3710  IF(K.CT.LDSE)  GO  TO  2510 

3720  SEARCP(K-l)=I2600 

3730  GO  TO  2595 

3740  2510  IF(K.GT.LDLRU)  GO  TO  2520 

3750  LnARC(K-LDSE)«I2603 

3760  GO  TO  2595 

3770  2520  SDARC(K-LDLRU)=I2600 

3780  GO  TO  2595 

3790  2530  IF( J .LE .LDSE)  GO  TO  2595 

3300  IF(J.CT.LQLRU)  GO  TO  2550 

3810  IF(K.LE.LOSRU)  CO  TO  2540 

3820  LSARC(J-LOSE)-I2600 

3330  GO  TO  2595 

3840  2540  3RDARC(K-LOLRU)*I26GO 

3850  CO  TO  2595 

3860  2550  IF( J .GT.LnSRU)  GO  TO  2560 

3870  SSAPC(J-LOLRU) *12600 

3880  GO  TO  2595 

3890  2560  IF(J.GT.LBLRU)  GO  TO  2580 

3900  IF(K.NE.LNODE)  GO  TO  2570 

3910  L3ARC( J“LDSRU)“T2600 

3920  GO  TO  2595 

3930  2570  IF(K.GT.LBSRU)  CO  TO  2595 

3940  SBSARC (K-L3LRU)*I2600 

3950  GO  TO  2595 

3960  2530  IF(J.GT.LBSRU)  CO  TO  2590 

3970  IF(K.NE.LNODE)  GO  TO  2595 

3980  SBARC(J-LBLRU)=I2600 

3990  GO  TO  2595 

4000  2590  SEARCP(J+NDSE-LBSRU) =12600 

4010  2595  FLO W(1 2600 )*0 

4020  2600  CONTI MU E 

4030  3000  CONTINUE 

4040C  IF(MSGL2.LT.2)  GO  TO  2440 

4050C  DO  2410  I«l,LFMS 

4060C  PRINT, 'LPTRS  ' ,LDARC(I) ,LSARC(I) ,LBARC(I) 

4070C  2410  CONTINUE 
40S0C  00  2420  I=1,SFM3 

4090C  PRINT, 'SPTRS  ' ,S0ARC ( I ) ,SS ARC ( I ) ,SBARC(I) ,SBDARC(I) ,SBSARC(I) 

4100C  2420  CONTINUE 

41  IOC  00  2430  1=1 ,NUMSER 

4120C  PRINT, 'SEPTRS  ',SF.ARCP(I) 

4130C  2430  CONTINUE 
4140C  2440  CONTINUE 
4150  RETURN 

4160  END 


SUBROUTINE  NMXFLO(MSGL2 ) 


10 
20C 
30C 

40C  TEMPORARY  DATA  INPUT  AREAS 

50  INTEGER  MIRROR, NSTART ,  IAPRF.V  .NAPATH  ( 300) 

60C 

70  IF  (ICIRF.EQ.2)  GOTO  3012 

80  DO  3005  13005*1 ,LFMS 

90  IF  (OPnECL(I3005) .NE.LOCAT( 1 5) )  GOTO  3005 

100  IF(LFMOD(I3005) .EO.O)  GO  TO  3001 

110  J  300  5  *LU ARC (13005) 

120  CAP (J3005)* JUMBO 

130  3001  IF(LFMOS(I3005) .EO.O)  GO  TO  3002 

140  J3005-LSARC(I30Q5) 

150  CAP(J3005)- JUMBO 

160  3002  IF(LFMOB(I3005) .EQ.O)  GO  TO  3005 

170  J  3005-LBARC ( 1 3005 ) 

180  CAP (J3005) -JUMBO 

190  3005  CONTINUE 

200C 

210  IF(SFMS.LT.l)  GO  TO  3012 

220  DO  3010  T3010-1.SFMS 

230  IF  (OPDF.CS  (13010)  .NE.L0CAT(15)  )  GOTO  3010 

240  IF(SOD(I3010) .EQ.O)  GO  TO  3006 

250  J3010-SDARCU3010) 

260  CAP (J30 10) -JUMBO 

270  3006  I F( SOS (I 30 10) .EO.O)  GO  TO  3007 

280  J3010»SSARC( 13010) 

290  CAP (J30 10) -JUMBO 

300  3007  IF(SOB(I3010) .EO.O)  GO  TO  3010 

310  J30 1 0-SBARC (13010) 

320  CAP (J30 10) -JUMBO 

330  3010  CONTINUE 
340  3012  CONTINUE 
350  CALL  PTIME(X) 

360C 

370C  DEPTH  FIRST  MAX  FLOW  ALCOTITHM 

380C 

390C  DATA  USED  IS  CONTAINED  IN  THE  VECTORS  SRCE,  DEST  ,CAP,  FLOW 

400C  AND  BKPTR  WHICH  HAVE  AN  ENTRY  FOR  EACH  ARC;  AND  IN  VECTORS 

410C  N'PATH,  DLTAFL ,  STATE,  FWDSP,  AND  BACKSP  WHICH  HAVE  AN  ENTRY  FOR 
420C  EACH  NODE.  STATE  IS  A  'l"  FOR  AN  UNLABELED  NODE  AND  A  '2'  FOR  A 
430C  LABELED  NODE. 

440C 

450C  LABEL  THE  SOURCE 

460C 

470  DLTAFL ( 1 )  -  JUMBO 

480  MIRROR  -  0 

490  NPATH(l)  -  -1 

500  5015  STATE( I )  -  2 


5  IOC 

520C  LABEL  ALL  OTHER  NODE*?  AS  UNLABELEO  '1' 

530C 

540  DO  5020  15020  -  2.LN0DF. 

550  STATE(T5020)  -  l 

560  5020  CONTINUE 

570C 

580C  FIND  AN  ARC  FROM  NODE  1  WITH  CAPACITY  REMAINING 
590C 

600  J5050  -  FWDSP(t) 

610  K5050  -  FWDSP(2)  -  1 

620  DO  5060  15060  -  J5050.K5050 

630  NODE  -  DEST( 15060) 

640  IF  ( FLO W( 15060) .LT.CAP ( 1 5060 ) .AND .STATE (NODE) . EO . 1 )  GOTO  5065 

653  GOTO  5060 

660C 

670C  MARK  FORWARD  OPATH  AND  CALC  MAX  FLOW  ON  ARC  15060 
680C 

690  5065  NPATH (NODE)  *  1 

700  STATE (NODE)  -  2 

710  NAP ATH (NODE)  -  15060 

720  JDELTA  *  CAP(I5060)  -  FLOW(I5060) 

730  DLTAFL(NODE)  -  MT NO  (JDELTA, DLTAFLd  )  ) 

740C 

750C  CONTINUE  FORWARD  DEPTH  FIRST  FLOW  AUGMENTATION 
760C 

770  5070  NNODE  -  NODE 

780  J5030  -  FWDSP (NNODE) 

790  K5030  -  FWDSP (NNODE  +  I)  -  1 

800  DO  5090  15090  -  J5080.K5080 

810  NODE  -  DEST(I5090) 

820  IF  (FLOW(I5090)  .LT.CAP (15090)  .AND  .STATE (NODE)  .  F,0 . 1 )  GOTO  5085 

830  GOTO  5090 

840  5085  NPATH (NODE)  -  NNODE 

850  ST,\TE(NODE)  -  2 

860  NAPATH(NODE)  -  15090 

870  JDELTA  -  CAP(I5090)  -  FLOW(I5090) 

880  DLTAFL(NODE)  -  MINO(JDELTA,DLTAFL(NNODE) ) 

890C 

900C  IF  NODE  EQUALS  LAST  NODE  INCREMENT  FLOW  ON  PATH 
910C 

920  IF  (NODE.EQ.LNODE)  GOTO  5150 

930  GOTO  5070 

940  5090  CONTINUE 
950C 

960C  FIRST  TIME  THROUGH  SKIP  MIRROR  ARC  SECTION 
970C 

930  IF  (MIRROR. EQ.O)  GOTO  5140 

990C 

1000C  MIRROR  ARC  AUGMENTING  PATH  SECTION 
I010C 

1020  5100  J510Q  «  RACKSP(NNODF.  -  l)  +  1 

1030  15110  -  RACKSP(NNOOF) 


1040  no  5130  15130  -  J5100.J5110 

1050  K5100  -  BKPTR(I5130) 

1060  NODE  -  SRCE(K5 100) 

1070C 

1030C  CHECK  FOR  SOURCE  NODE  OR  PATH  JUST  TRAVELED 
1090C 

1100  IF  (NODE. E'3.1. OR. NODE. EO.NPATH(NNODE) )  GOTO  5130 

1110  IF  (FLOW(K5 100) .GT .0 . AND .STATE(NODE) . EQ . 1 )  GOTO  5120 

1120  GOTO  5130 

1 1 30C 

1140C  MARK  REVERSE  PATH  AND  CALC  MAX  FLOW  ON  ARC  K.5100 
115DC 

1160  5120  NPATH(NODE)  =  0  -  NNODE 
1170  STATE (NODE)  -  2 

1180  NAPATH(NODE)  -  K5100 

1190  DLTAFL(NODE)  -  MIN0(FLOW(K5100) .DLTAFL(NNODE) ) 

1200  GOTO  5070 

1210  5130  CONTINUE 
1220C 

1230C  IF  UNABLE  TO  AUGMENT  PATH  TO  SINK  FROM  THIS  NODE 

1240C  LABEL  IT  ' 2 '  AND  BACKUP  TO  PREVIOUS  MODE  ON  PATH 

1250C 

1260  5140  NODE  =  IABS (NPATH (NNODE) ) 

1270  IF  (NODE.E0.1)  GOTO  5060 

1280  GOTO  5070 

1290C 

1300C  FLOW  AUGMENTATION  ON  PATH  FROM  SOURCE  TO  SINK 
1310C 

1320  5150  LAST  -  LNODF. 

1330  STATE (LAST)  *  l 

1340  NSTART  *  0 

1350  INC  *  DLTAFL  (LNODF.) 

1360  5160  IPREV  -  NPATH (LAST) 

1370  IAPP.EV  -  NAPATH(LAST) 

1380  IF  (IPREV. GT.O)  GOTO  5170 

1390  FLOW(IAPREV)  -  FLOW(IAPREV)  -  INC 

1400  GOTO  5180 

1410  5170  FLOW(IAPREV)  -  FLOW(IAPREV)  INC 

1420  5180  LAST  -  IABS(IPRF.V) 

1430  IF  (LAST.EQ. I )  GOTO  5200 

1440  DLTAFL(LAST)  -  DLTAFL (LAST)  -  INC 

1450C 

1460C  IDENTIFY  AND  MARK  THE  LAST  NODE  IN  THE  PATH  WITH  EXCESS 

1470C  FLOW  REMAINING  FOR  AUGMENTATION  TO  THE  SINK 

1480C 

1490  IF  (FLOW(IAPREV) .LT.CAP(IAPREV) .AND. DLTAFL (LAST) .E0.0) 

1500  &  STATE(LAST)  -  l 

1510  IF  (NSTART. EO.l. OR. DLTAFI.(LAST).EQ.O)  GOTO  5190 

1520  NSTART  -  1 

1530  NODE  «  L4ST 

1540  5190  GOTO  5160 


1550  5200 
1560  5060 
1570C 
1530C 
1590C 
1600 
1 6 10C 
1620C 
1630C 
1640C 
1650 
1660 
1670C 
1680C 
1690C 
1700C 


IF  (NSTART.EQ. 1 )  GOTO  5070 
CONTIMUE 


AFTFR  SECOND  TIME  THROUGH  SOLUTION  COMPLETE 


IF  (MIRROR. EO. 2)  COTO  3300 


SET  SWITCH  TO  ALLOW  FOR  FLOW  AUGMENTATION  ON  MIRROR  ARCS 
AND  START  OVER 


MIRROR  -  MIRROR 
GOTO  5015 


+  1 


**** 


****  NETWORK  SOLUTION  COMPLETE 
COMPUTE  MAX  FLOW 
1710  3300  J3310-1 
1720  K3310-FWDSP(2)-1 

1730  ITFLOW-O 

1740  DO  3310  I3310-J3310.K3310 

1750  ITFLOW-ITFLOWfFLO  W(  I 3310) 

1760  3310  CONTINUE 
1770  FLOW(MAXARC+l)-ITFLOW 

1780C  IF(MSGL2.LT.l)  GO  TO  3400 
1790C  PRINT,'  TOTAL  FLOW  '.ITFLOW 
1800C  PRINT  2220, 'NODE  #',(1,1*1 ,LMODE) 

18 IOC  PRINT  2220, 'SRCE  ' , (SRCE(I) , r-l ,NARCS) 
1820C  PRINT  2220, 'DEST  ' , (DEST(I) ,1*1 ,NARCS) 
1830C  PRINT  2220,'NPATH  ' , (NPATH(I) ,1*1 ,LNODE) 
1840C  PRINT  2220,'DLTAFL' ,  (DLTAFT.(I)  ,1-1  ,LNOPE) 
1850C  PRINT  2220, 'STATE  ' , (STATE(I) ,1-1 ,LNODE) 
1860C  PRINT  2220, 'FLOW  ', (FLOW(I) ,1-1 .MARCS) 


1870C  PRINT  2220, 'CAP 
1880C  3400  CONTINUE 


' , (CAP(I) ,1-1 ,NARCS) 


1890 

1900 

1910 

1920 

1930 

1940 

1950 


GO  TO  3330 


GO  TO  3320 


LCSEC-0 
NDSEP1-NDSE+1 
IF(NDSE.EQ.O) 

DO  3320  13320-2, NDSEP1 
IF (STATE (1 3320) .NE. 1 ) 

S  EPTR-S  EARC? (13320-1) 

LCS EC-LCS  EC+CAP (SEPTR) 

I960  3320  CONTINUE 

1970  3330  IF(NDSE.EO.NUMSER)  GO  TO  3350 
DO  3340  I3340-NDSEP1 .NUMSER 
J3340-L3SRU+I3340-NDSE 
IF(STATE(J3340) .EQ.i)  GO  TO  3340 
SEPTR-SEARCP(I3340) 
LCSEC-LCSEC+CAP (SEPTR) 

2030  3340  CONTINUE 

2040  3350  FL0W(MAXARC+2)-LCSEC 

2050C  J-0 


1930 

1990 

2000 

2010 

2020 


2 


2060C  00  3360  I-l.NARCS 

2070C  IF(0EST(T) .  E0.LN0DE)  J-.J+FLOW(I) 

2080C  3360  CONTINUE 

2090C  IF(J.NE.ITFLOU)  PRINT, 'NOT  CONSISTENT  FLOW  '.ITFLOW.J 

2100  Y  -  0 

2110  CALL  PTIME(Y) 

2120  Z-Y-X 


2130  TT-TT+Z 

2140  WRITE (6 ,400)  X*36OO,Y*36OO,Z*36OO,TT*36O0 

2150  400  FOP.MAT ( '  TIME-' ,4F10.5) 

2160  RETURN 

2170  END 


Appendix  G;  OUTPUT  Subroutine 


10  SUBROUTINE  OUTPUT 

20C 

30C 

40C 

50C  OUTPUT  RESULTS  FOR  DEPOT  SE 

60  3400  I P AG E= I PAG E+ 1 

70  '..'RITE (6 ,3401 )  (SYSNAM(ISUB)  ,ISUB=1 ,3)  ,DATE,IPAGE 

80  3401  FORMAT ('l' , 2K , 3A4 , R6X , AS , 5X , ' PACE'  ,14) 

90  WRIT£(6 , 3402 ) 

100  3402  FORMAT ( '  ' ,24X , 'NETWORK  RLA  RESULTS') 

110  WRITE (6 , 3403 ) 

120  3403  FORMAT ('  ' ,90X,' WHOLESALE  CHANGE  FACTORS') 

130  WRITE (6, 3404)  WCFUC ,WCFF ,WCFSE 

140  3404  FORMAT  ( '  ',96X,'COST  -'  ,F6.2/97X,'MTBF  ='  ,  F6 . 2 /99X , '  SF.  =',F6.2) 
150  WRITE(6,3410) 

160  3410  FORMAT!'  ' ,22" ,' SUPPORT  EQUIPMENT  REQUIPEMENTS ' ) 

170  WRITE(6 ,3412) 

130  3412  FORMAT! '0' ,4X, 'DEPOT' ,2*X, 'TOTAL  SE') 

190  WRITE (6 ,3415) 

200  3415  FORMAT!'  ','SE  CODE'  ,  5X  ,  'NAME'  ,  1  OX ,'  RF.OMT'  ,  5X  ,'  LC  S/BASE', 3X, 
210  &  'USE  HRS' ,7X, 'USE  %') 

220  TOTSEn-O. 

230  IF1NDSE.LT. 1)  GO  TO  3431 

240  DO  3430  13430*1, NDSE 

250  IF (STATE (1 34 30+1) ,GT. 1 )  GO  TO  3430 

260  WRITE (6 , 3420 )  SECODE! 1 3430 ) , (XSF ( T 3430 , TSUB) , ISUB-1 , 3 ) , 

270  h  REQMT  (13430)  ,TSF.CST  (13430)  ,USFHRS (13430)  ,SEUP(I3430) 

280  3420  FORMAT!'  ' ,2X,I5 ,5X , 3A4 ,3X , 13 , IX, FI 2 .0 ,F1 2 .0 ,2PF1 2 . 1 ) 

290  TOTSEO*TOTSED+TSECST( 13430) 

300  3430  CONTINUE 

310  3431  CONTINUE 

320  WRITE (6 , 3435 )  TOTSED 

330  3435  FORMAT ('0' ,2X, 'TOTAL' ,24X,F12.0) 

340C 

350C  OUTPUT  RESULTS  FOR  BASE  SUPPORT  EQUIPMENT 

360  IF(NDSE.EQ.NUMSFR)  CO  TO  3451 

370  WRITE (6, 34 40) 

380  3440  FORMAT! '0' ,4X, 'BASE' ,29X, 'TOTAL  SE') 

390  WRITE (6,3415) 

400  TOTSEB=0. 

410  NDSFP1-NDSE+1 

420  MBS  F.  POUNDS  E+NBSF 

430  DO  3450  I3450-NDSEP1 ,NBSEP0 

440  J3450*LBSPU+I3450-NDSE 

450  I F( STATE (.134 50)  . EO. I )  GO  TO  3450 

460  3445  WRITE(6 , 3420)  SECODE! 1 3450 ) , (XSE ( 1 3450 , ISU3) , ISU3- 1  , 3 ) , 

470  &  PEQMT(I3450) ,SECOST(T3450) .USEURS (13450) ,SEUR(I3450) 

480  TOTS  EB*TOTSEB+TS ECS T (13450) 

490  3450  CONTINUE 

500  WRITE (6,3435)  TOTSEP, 


126 


510r  OUTPUT  RES'U.TS  FOR  CTRp  S'TPOPT  F.ouir“ir‘:r 
520  3442  TF  (NBSEPO.EO. NUMBER)  GOTO  3451 
530  WRITE(6,3SOO) 

540  3900  FORMAT ( '0' ,4X, 'CIRF' ,29X , 'TOTAL  SF') 

550  WRITE  (6,3415) 

560  •  IF  (ICIRF.EO.l)  WRITE  (6,3810) 

570  3310  FORMAT  ( '  NO  CIRF  SF.  CALCULATE")  THIS  PASS') 

530  IF  (ICIRF.EO.l)  GOTO  3451 

590  TOTSEC=0 

600  NBSEP1=NBSEP0+1 

610  DO  3920  I3820=MRSFP1 .NUMSER 

620  J3920=I3^20-NBSEP0 

630  I*  (STATF(J3320) .CT.l)  GOTO  3820 

640  WRITE  (6,3420)  SEC0DE(  1 38  20 )  ,  (XSF.  ( 13920  ,  ISUB)  ,  TSUB=  1  ,  3 )  , 

650  &  P.EQMT (13820)  ,TSECST(I3320)  ,USEHRS  ( 13820 )  ,SF.!’R(T 3820 ) 

660  TOTSEC=TOTS  F.C+TSECST (13820) 

670  3820  CONTINUE 

680  WRITE  (6,3435)  TOTSEC 

690C 

700C  OUTPUT  RESULTS  FOR  LRU'S  &  SRU'S 

710  3451  TPAC,F.=I PACE+1 

720  WRITE ( 6,3401)  ( SYSNAM(ISUB) ,  ISUB=1 ,3) .DATE , I PACE 

730  WRITE(6 , 3452 )  CHARS (2 ) .CHARS (2 ) 

740  3452  FORMAT ( '  ',40X,'LRU  &  SRU  REPAIR  LEVEL  DECISIONS' ,2X,2A4) 
750  WRITE(6,3453) 

760  3453  FORMAT  ('0'  ,  29X , '  — —  LRU  (LC  S/BASE)  — ',7X,'—  SRU 
770  &  '(LC  S/BASF.)  — '  ,  7X , 'LRU/SRII  LRU' ,7X,  '  INPUT  MTE 

780  &  'LC  REP'/'  ' ,4X,'LRU/SRU  NAME  IDENT  MO  DEPOT  SCF 

790  &  '  BASF.','  CIRF', 'DEPOT  SCRAP  BASE','  CIRF 

800  &  'COST  FAIL'S' ,7X,'MTBF  REPAIR  DEM/BASE') 

810  TOTLO=0 . 

820  T0TLS=0 . 

830  TOTLC=0. 

840  TOTLB-O. 

850  TOTSD=0, 

860  TOTSS-O. 

870  TOTS  8=0 . 

880  LINES =3 

890  UUC ( l ) “LOCAT (15) 

900  WUC ( 2 ) =L0CAT (15) 

910  NOSEP2=NDSE+2 

920  MSE=MDSEP2 

930  ICTRF-l 

940  LLRU=LDLRU 

950  DO  3600  1 3600-MS  F..LLRU 

960  ITF.M=I3600-(MSE-1 ) 

970  I F (LFMVUC ( ITEM,  1 ) .  EQ.L1,C(  l ) .  AND .  LFMWUC ( ITEM , 2 ) . EO .  WUC  ( 2 ) ) 

930  &  CO  TO  3457 

9°0  DO  3454  13454=1, MLRU 

1000  TE(T.FMWT’C(TTE*',1  )  .EO  .T.WUC  (13454,1  )  .  A’D  . 

1010  &  LFMWl'C  ( TTEM ,  2  )  .  EO  *  LWUC  (13454,2)'  LPTD  =  I3454 

1020  3454  CONTINUE 

1030  MTBCT=MTBF  (LPTP )  /  (UF(T.D7R)*(  l  .-RIP(LPTR) ) ) 


1040 

1050 

1060  3457 
1070 
1080 
1090 
1100 
1110 

1129  3458 

1130 
1140 
1150 
1160 
1170 
1130 
1190 
1200C 
1210 
1220C 
1230 
1240C 
1250 
1260 

1270  3459 

1230 

1290C 

1300  3460 

1310 

1320 

1330 

13  40 
1350C 
1350  3470 
1370 
1330 
1390 
1400 

14  IOC 
1420  3480 
1430 
1440 
1450 
1460 

1470  34«3 

1430 

1490 

1509 

1510  3490 


TOCTOM=PGMc<*'M,A  ( !.n  TR  )  /MT  3CT 
TLCn=T0CTCM*PTUP*12. 

I  FO.INES  .LT.  54 )  GO  TO  3458 
I P  AO  F.=  I  PAG  E+ 1 


,.’RITE(6,3401  )  (SYSNAM(ISUB)  .TSUB-l  ,3)  .DATE  ,  IPAGF. 
WRITE  (6,345  2)  CHARS ( 5 ) .CHARS (6 ) 

WRITE(S,3453) 

LINF.S-4 

TLCDF=TLCD*FAILP (ITEM) 

FMTBCT-MTBCT/FAILP (ITEM) 

NODEL-I360G 
NOD  F.  R = NOD  EL+LFMS +S  FMS 

IF  (OPDF.CL(  ITEM)  . EO .LOCAT (4)  )  GOTO  3483 
IF  (OPDECLC ITEM) . EO .LOCAT ( 1 ) )  GOTO  3460 
IF  (OPDECL(ITEM) .EQ.LOCAT(2) )  GOTO  3470 
IF  (OPDECL (ITEM) . EO .LOCAT ( 3 ) )  GOTO  34pO 
CHECK  FOR  LRU  DEPOT  REPAIR 
I F( STATE (MODEL) . EO . 1 .AND .STATE (MODER ) . EO . 1 ) 

CHECK  FOR  LRU  SCRAPPED 
IF (STATE (MODEL) .NE.l .AND.STATE(NODER) .FQ. 1 ) 

CHECK  FOR  LRU  RASE  REPAIRED 
I F (STATE (MODEL ) .ME . 1 . AN  D . STATE (NOD  ER ) .NE.l) 
URITE(6,3459) 

FORMAT ('  TROUBLE  ON  LRU  LABELS') 

STOP 

COSTL=DEPOLC ( ITEM) 

TOTLD-TOTLD+COSTL 
LEVEL- 1 

0?9ECL( ITEM) -LOCAT (1 ) 

GO  TO  3490 


COSTL-SCRPLC ( ITEM) 

TOTLS-T0TLS+COSTL 

LEVEL-2 

OPDECL ( ITEM) -LOCAT ( 2 ) 
GO  TO  3490 

COSTL-BASELC( ITEM) 
TOTLB-TOTLB+COSTL 
LEV EL- 3 

OPDEi  .. ( ITEM) -LOCAT (3) 
GOTO  3490 

COSTL-CI RFLC ( ITEM) 

TOTLC-TOTLC+GOSTL 

LEVEL-4 

OPDECL (ITEM) -LOCAT (4) 
FPCT-FAILP( ITEM) *100, 
LOD-CHARS (2) 

1.9 S -CHARS  (2  ) 

LOD-CHARS (2) 

LOC-CMAES ( 2 ) 

I E ( LFMOO ( ITEM) . EO . 2 ) 


GO  TO  3460 


GO  TO  3470 


GO  TO  3480 


•CHARS (14) 


x/’-VS.O 


,  •  . 


1  570  l F(LFMOS (ITEM) .  T.  ? .  2  )  i.or=CUAR3  ( 1 4  ) 

1530  IF(LFMOB(ITEM) .EQ.2)  L03=CUARS  ( 1  4 ) 

1590  IF  (LFMOC(ITFM)  .EQ.2)  LOOCFARS  ( 1 4 ) 

1600  IF(SRUPTR(  ITEM)  .UE.O)  GO  TO  3520 

1610  GO  TO  (3493, 3500, 3510, 3580), LEVEL 

1620  3493  IF(VTCO)  .ME.  I.  FM  WUC  ( I T  FM ,  1 )  .OR .  WUC  (2  ) .  N'F  .LFMM’C  ( ITEM ,  2  ) ) 

1630  &  GO  TO  3496 

1640  3494  WRITE (6, 349 5)  COSTL.LOS  ,LOP,,LOC,FPCT,FMT3CT,TLC9F 

1650  3495  FORMAT ('  ' ,27X,F8 .0 , 3X ,A4 ,4X ,A4 ,2X , A4 ,38X,F5 . 1 , 1 3X ,F8 .0 ,2X ,FS . 1 ) 

1660  LIMES=LINF.S-K 

1670  GO  TO  3600 

1630  3496  IF(FRSTFM(LPTR) .EQ.LASTFM(LPTR) )  GO  TO  3498 

1690  WR*ITE(6  ,3497)  LPTR  .LRUNA'l  (LPTR ,  1 )  ,LRUNAM(LPTR  ,2)  ,LRUNAM(LPTR  ,3) 

1700  &  LWUC(I,?TR,1)  ,LVUC(LPTR,2)  .UCL(LPTR)  .MTBF(LPTR)  ,MTBCT ,TLC0 

1710  3497  FORMAT ('0' , 12 ,2X , 3A4 , 2X, A4 ,A3 ,59X ,F8 .0 , 10X , F8 .0 , 2X , F3 .0 , 

1720  &  2X.F8.1) 

1730  LINF.S-LINF.S+2 

1740  WUC ( 1 )=LFMVUC (ITEM, l ) 

17  50  WUC  ( 2 )  =LFMWUC  (ITEM,  2) 

1760  GO  TO  3494 

1770  3498  WRITE(6 ,3499)  LPTR,LRUNAM(LPTR,  1 )  ,L?.UNAM(L?TR  ,2 )  ,LRUM AM(LPTR ,3 ) 
1780  &  L WUC (LPTR, 1 ) ,L WUC (LPTR, 2) , COSTL.LOS ,LOB ,LOC ,UCL (LDTP ) , 

1790  S,  FPCT.MTBF(LPTR) , FMTRCT ,TLCDF 

1800  3499  FORMAT (*0',I2,2X,3A4,2X,A4,A3,2K,F8.0,3X,A4,4X,A4,2X,A4,28X,F8.0 
1810  &  2X.F5.1 ,3X,F8.0,2X,F8.0,2X,F8.1) 

1820  LIMES-LINES+2 

1830  WUC(  1)»LFMLT1C(  ITEM,  1) 

1840  WUC  ( 2 )  “LFMVUC  ( I TFM ,  2 ) 

1850  GO  TO  3600 

1860C 

18  70  3500  I F ( WUC ( l ) . ME .LFMWUC (ITEM, 1 ) . OR . WUC (  2 ) . WE . LFM WUC (ITEM , 2 ) ) 

1830  &  GO  TO  3503 

1890  3501  VRITE(6,3502)  LOO ,COSTL ,LOB ,LnC ,FPCT .FMTBCT ,TLCOF 

1900  3502  FORMAT ( '  ' , 30X ,A4 , IX ,F8 .0 ,3X ,A4 ,2X, A4 , 38X , F5 . I , 1 3X , F8 . 0 , 2X , F8 . 1 ) 

1910  LINES-LINES+l 

1920  GO  TO  3600 

1930  3503  IF( FPSTFM(LPTR)  .F.O.LASTFM(LPTR) )  GO  TO  3504 

1940  WRITE(6,3497)  LPTR,LRUNAM(LPTR, 1 ) ,LRUWAM(LPTR , 2 ) ,LRUUAM(LPTR , 3 ) 

1950  &  L WUC (LPTR, 1) ,LWUC(LPTR,2) .UCL(LPTR) .MTBF(LPTR) .MTBCT.TLCD 

1960  LIMES*!, IKES+2 

1970  WUC ( 1 )=LFMWUC( ITEM, l ) 

1980  WUC  ( 2 )  -LFMWUC  ( I T  EM ,  2 ) 

1990  GO  TO  3501 

2000  3504  '.TRITE (6 ,3505)  LPTR, LRl'NAM( LPTR,  L  )  ,LRUMAM(LPTR  ,2)  ,  LRUS AM (LPTR  ,3) 
2010  &  L  WUC  (LPTR,  1 )  ,LWUC(LPTR, 2 )  ,L0D,C0S7L  ,L03  ,L0C , 

2020  &  'JCL(LPTR)  ,FPCT,MTBF(LPTP)  .FMTBCT .TLCUF 

2030  3505  FORMAT ( ' O' , 12 ,2X , 3A4 , 2X , A4 ,A3 , 5X , A4 , IX , F8 .0 , 3X , A4 ,2X ,A4 , 28X , 

2040  &  F8.0,2X,F5.1,3X,F8.0,2X,F3.n,2X,F8.1) 

2050  LINES-LIUES+2 

2060  WUC(  1  )*LFMWUC  (ITEM,  1 ) 

20  70  'VC ( 2  )  =L FMWUC ( ITEM , 2  ) 

2030  GO  TO  3600 

2090C 


2100  3510  IF(LUC(l)  .MF.LFMWC(TTE'1,  l )  .  OR  .  LVC  (  2  )  .ME  .I.FMU'C  (  T  Tr  “ ,  2  )  ^ 

2110  &  00  TO  3513 

2120  3511  UP.ITE(8,3512)  LCD  ,LOS  ,OOSTL  ,LOC  , FPCT , FMTBCT.TLCDF 

2130  3512  FORMAT ( '  ' , 30X ,A4 ,4X, A4 , 1 X ,F8 .0 , IX ,A4 , 38X, F5 . 1 , 1 3X , FS .0 , 2K,FS . 1 ) 

2140  LIUES=LIN’ES  +  1 

2150  .  00  TO  3500 

2150  3513  IF(FRSTFM(LPTR) . EQ.LASTFM(LPTR) )  00  TO  3514 

2170  KRITF.(6 ,3497  )  LPTR ,LRUMAM(LPTR , 1 )  .LRUNAM (LPTR ,  2 )  , LRIJUAM (LPTR ,  3  )  , 

2180  &  L’.7»2C(LPTR,1)  ,L'.vi;C(LPTR,2)  .UCL(LPTR)  , 

2190  &  MTSF(LPTR)  ,MTFCT,TLCr> 

2200  LINES-LINES+2 

2210  WITC  ( l )  “I.FMWUC ( ITEM ,  1 ) 

2220  WUC ( 2  )=LF>U.T'C ( ITEM, 2  ) 

2230  00  TO  3511 

2240  3514  WRITF(6,3515)  LPTR ,LRU>iAM(LPTR , 1 ) ,LRUMAM(LPTR , 2 ) ,LRUNAM(LPTR , 3) , 
2250  &  LWl'C  (LPTR ,  1 )  ,LWUC (LPTR ,  2 )  ,LOD  ,LOS  ,COSTL  ,L0C  , 

2260  &  'JCL(LPTR)  ,FPCT,MT3F(LPTR) ,  FMTBCT , TLCOF 

2270  3515  FORMAT ( '0' ,I2,2X,3A4,2X,A4,A3,5X,A4,4X,A4,1X,F8.0,1X,A4,28X,F8.0,2X 
2280  &  F5.1 ,3X,F8.0,22<,F8.0,2X,F8.1) 

2290  LIN’ES-LINES+2 

2300  lv’i,C(l)»LFM,..rC(ITEM,l> 

2310  LUC  ( 2 )  “LFMWUC  ( ITEM ,  2 ) 

2320  GO  TO  3600 

2330C 

2340  3580  IF( WUC(  1 )  .NE.LFMWUC(ITF.M,1 )  .OR.  WUC(2)  .NE.LFMWUC (ITEM, 2 ) ) 

2350  &  GO  TO  3583 

2360  3581  WRITE (6 , 3582 )  LOO , LOS , LOB, COSTL, FPCT .FMT3CT , TLCOF 

2370  3582  FORMAT('  ' , 30X.A4 ,3X,A4 ,3X,A4 , IX, F7 .0,3SX,F5 . 1 , 1 3X , F8 .0 , 2X,F8 . 1 ) 

2380  LINES-LINES+1 

2390  GO  TO  3600 

2400  3583  IF(FRSTFM(LPTR) . EQ .LASTFM(LPTR) )  GO  TO  3534 

2410  WRITE(6 , 3497 )  LPTR , LRUN’AM(L?TR,  1 )  ,LR1JN’AM(LPTR ,  2  )  ,LRU>JAM(L?TR ,  3 ) , 

2420  &  LUUC(LPTR, l) ,LWUC(LPTR,2) ,UCL(LPTR) , 

2430  &  MTBF(LPTR) ,MT3CT,TLCD 

2440  LIMES-LINES+2 

2450  WUC  ( 1 )  *LFMWTJC  ( ITEM ,  1 ) 

2460  WC  ( 2 )  “LFMWUC  ( ITEM ,  2 ) 

2470  GO  TO  3581 

2480  3584  WRITE(6,3585)  LPTR,LRUXAM(LPTR , 1 ) ,LRUNAM(LPTR,2) ,LRUNAM(LPTR,3) , 
2490  &  LWnC(LPTR.l) , LWUC(LPTR , 2 ) ,LOD ,LOS ,LOB .COSTL , 

2500  &  UCL(LPTR) , FPCT, MTBF(LPTR) .FMTBCT, TLCDF 

2510  3585  FORMAT ('0' ,12 ,2X ,3A4 ,2X,A4 ,A3,5X,A4 , 3X,A4 , 3X,A4 , 1X,F7 .0,29X, 

2520  &  F8.0,2X,F5.1,3X,F8.0,2X,F8.0,2X,F8.1) 

2530  LIt!ES“LIMES+2 

2540  W17C(  1  )=LFMWUC(ITEM,  1 ) 

2550  LUC  (2 )  “LFMWUC  (ITEM ,  2 ) 

2560  GO  TO  3600 

2570C 

2530  3520  IF(WUC(  1 )  •  F.O.LFMLUCdTEM, l )  .AND. WUC ( 2 ) . FO .LFMtUC ( ITEM, 2 ) ) 

2590  &  GO  TO  3525 

2600  WRITE  (6, 3497)  LPTR,  LRUXAM  (LPTR,  1 )  ,LRUNAM(LPTR,2)  ,LRUNAM(LPTR ,3)  , 

2610  &  LWUC(LPTR, 1 ) ,LWUC(LPTR ,2) ,UCL(LPTR) , 

2620  &  MTRF(LPTR) ,MT8CT,TLCO 
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2630  LINES-LINF.S+2 

2640  MTC(1  )=LFV "UC  ( ITEM  ,  1 ) 

2650  WC(2)«LFMWnc(lTEM,2) 

2660C 

2670  3525  NODEL-LDLRU+S RUPTR( ITEM) 

2630  S'ODER-LBLRU+SRUPTR  (ITEM) 

2690  ITFM=SRUPTR( ITEM) 

2700  ISOD-CUARS  (2 ) 

2710  IS0S-C!1ARS(2) 

2720  ISOB-CHARS(2) 

2730  ISOC«CUARS(2) 

2740  I F(SOD ( ITEM) • EO . 2 )  1300-CHARS ( 14 ) 

2750  T  F(SOS (ITEM) . EO. 2 )  ISOS-CHARS ( 1 4 ) 

2760  I F(SOB( ITEM) . EO. 2 )  ISOS-CHARS (14) 

2770  IF(S0C(ITEM).EO.2)  TS0C-C'IARS(14) 

2730  IF  (OPDECS ( ITEM) . FQ. LOCAT (4 ) )  GOTO  3630 

2790  IF  (ICIRF.FQ.2  .AND. STATE (MODEL)  .  FQ .  1 )  GOTO  3630 

2800  IF  (OPDECS(ITEM)  . EQ. LOCAT (1 ) )  COTO  3530 

2810  IF  (OPDECS ( ITEM) .EQ.LOCAT(2) )  GOTO  3550 

2820  IF  (OPDECS ( ITEM) . FQ . LOCAT ( 3) )  GOTO  3570 

2830C  CHECK  FOR  SRU  DEPOT  REPAIR 

2840  IF(STATE(NODEL)  .EQ.  1  .AND.STATE(NODER)  .EO.l )  CO  TO  3530 

2850C  CHECK  FOR  SRU  SCRAPPED 

2860  I  F( STATE (MODEL)  .N’E.  1  .AND.  STATE  (MODER)  .EO.l)  CO  TO  3550 

2870C  CHECK  FOR  SRU  BASE  REPAIRED 

2880  IF(STATE(NODEL).NE.l.AND.STATE(N'ODER).NF,.l)  GO  TO  3570 

2390  '.TRITE  (6 ,3526) 

2900  3526  FORMAT (*  TROUBLE  OM  SRU  LABELS') 

2910  STOP 

2920C  SRU  IS  DEPOT  REPAIRED 

2930  3530  COST-DEPOSC(ITEM) 

294CC  IF  LRU  IS  BASE  REPAIRED  ADD  EXTRA  COSTS 

2950  IF(LEVEL.E0.3)  C05T-C03T+BRLDS(ITEM)+RRLSS(ITEM) 

2960C 

2970  TOTSD-TOTSD+COST 

2930  OPDECS (ITEM) -LOCAT (1 ) 

2990  GO  TO  (3532, 3534, 3536, 3533), LEVEL 

3000  3532  VRITE(6, 3533)  S  RUN  AM  ( ITEM,  1 )  , SRUNAM  ( ITEM ,  2 )  ,SV/UC  (ITEM,  1 ) , 

3010  &  SWC  (ITEM,  2),  COSTL,  LOS,  LOB,  LOC,  COST,  ISOS,  IS03,TS0C,UCS  (ITEM) , 

3020  &  FPCT,FMTBCT,TLCDF 

3030  3533  FORMAT ( '  ' ,8X ,2A4 ,2X ,A4 ,A3 ,2X,FS .0 , 3X.A4 ,4X ,A4 , IX ,A4 ,FS .0 , 

3040  &  3X,A4 ,4X ,A4 ,2X,A4 ,F?.0 ,2X,F5 . 1 , l 3X ,F8.0 ,2X ,F3. I ) 

3050  CO  TO  3590 

3060  3534  WP.ITE(6,3535)  SRUNAM  (ITEM ,  I )  .SRUNAM  (ITEM,  2)  ,SUVC(ITEM,  l ) , 

3070  &  SUICdTEM, 2)  ,L0D, COSTL, LOB, LOC, COST, ISOS  , IS03 , ISOC ,UCS ( ITEM)  , 

3080  6  FPCT , FMT3CT ,TLCUF 

3090  3535  FORMAT ( '  ' ,3X ,2A4 ,2X ,A4 ,A3 , 5X , A4 , IX ,F8 .0 ,3X ,A4 , IX , A4 ,F3 .0 , 

3100  5.  3X,A4,4X,A4,2X,A4,F8.0,2X,F5.1 , 

3110  &  l 3X , F8 . 0 , 2X , F8 . I /40X , ' CASE  7  ERROR  IN  ABOVE  DECISION') 

3120  GO  TO  3590 

3130  3536  'TRITE  (6 , 3537  )  SRUNAM  (ITEM,  l  )  ,SRUNAM(ITF.M,2  )  ,SU'C(ITF.M,  1 )  , 

3140  6  SUDC ( ITEM , 2  )  ,L00  ,L0S  .COSTL  ,I.OC  .COST  ,  ISOS ,  ISOS,  ISOC  ,UCS  (TTEM)  , 

3150  &  FPCT , FMT3CT ,TLCDF 
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3160  3537  FORMAT ('  '  ,HX  ,2.A4 ,2X  ,A4  ,  A3  ,5\  ,A4 ,4V  ,A4  ,  IX  ,F 7 .0 , 1 X  ,A4  , F8  .0  , 

3170  &  3X,A4,4X,A4,2X,A4,F8.0,2X,F5.l, 

3180  &  13X,F8.0,2X,F3.1) 

3190  00  TO  3590 

3200  3538  LT.ITE(6 ,3539)  SRUNAM( ITEM , 1 )  ,SRTJ5iAM(ITEM,2)  ,SLT’C(  ITEM,  1 )  , 

3210  &  SWUC( ITEM, 2), LOD, LOS, LOB, COSTL, COST, ISOS, ISOB, ISOC, UCS (ITEM) , 

3220  &  FPCT , FMTBCT ,TLCDF 

3230  3539  FORMAT ( '  ' ,8X,2A4,2X,A4,A3,5X,A4,4X,A4,2X,A4,F7.0,F8.0, 

3240  &  3X,A4,4X,A4,2X,A4,F8.0,2X,F5.1,13X,F8.0,2X,F8.1) 

3250  GO  TO  3590 

3260C 

3270C  SRU  IS  SCRAPPED 

3280  3550  COST-SCRPSC(ITEM) 

3290C  IF  LRU  IS  BASE  REPAIRED  ADD  EXTRA  COST 

3300  IF(LEVEL.EQ. 3)  COST-COST+RRLSS (ITEM) 

3310  TOTS S -TOTS S+COST 

3320  OPDECS (ITEM) -LOCAT (2 ) 

3330  GO  TO  (3552,3554,3556,3558) .LEVEL 

3340  3552  WRITE(6 , 3553)  SRUMAM( ITEM, 1 ) ,S RUN AM (ITEM, 2) ,SWUC(ITE‘1, 1 ) , 

3350  &  SWUC(ITEM,2) ,COSTL ,LOS , LOB, LOC.ISOD, COST, ISOB , ISOC ,UCS ( ITEM) , 

3360  &  FPCT .FMTBCT ,TLCDF 

3370  3553  FORMAT ( *  ' ,SX,2A4 ,2X ,A4 ,A3 ,2X ,F8 .0 , 3X ,A4 ,4X , A4 , 2X ,A4 ,2X ,A4 , IX , 
3380  &  F8.0,3X,A4,2X,A4,F8.0,2X,F5.1 , 

3390  &  13X,F3.0,2X,F8.l) 

3400  GO  TO  3590 

3410  3554  WRITE(6,3555)  SRUMAM(ITEM,  1 )  ,SRl'MAM(ITEM,2)  ,S  LUC  (IT  EM ,  l ) , 

3420  S,  SWUC(ITEM, 2), LOD.COSTL.LOB, LOC.ISOD, COST, ISOB, ISOC, UCS (ITEM) , 

3430  &  FPCT, FMTBCT, TLCDF 

3440  3555  FORMAT ( '  ' ,8X,2A4,2X,A4,A3,5X,A4,1X,F8.0,3X,A4,2X,A4,2X,A4,1X, 
3450  &  F8 .0 ,3X ,A4 ,2X ,A4 ,F8 .0 ,2X  ,F5 . 1 , 

3460  &  13X,F8.0,2X,F8. 1) 

3470  GO  TO  3590 

3480  3556  WRITE(6 , 3557 )  SRUNAM(ITEM, 1 )  ,SRUNAM(ITEM,2)  ,SWUC(ITEM, 1 ) , 

3490  &  S  MIC (ITEM, 2), LOD, LOS, COSTL, LOC.ISOD, COST, IS03, ISOC, UCS (ITEM) , 

3500  &  FPCT, FMT3CT, TLCDF 

3510  3557  FORMAT ('  ' ,8X,2A4 ,2X,A4 , A3 ,5X ,A4 ,4X,A4 , IX, F8 .0 , IX ,A4 ,2X,A4 , IX , 
3520  &  F8.0 ,3X ,A4 ,2X,A4 ,F8 .0 ,3X ,F4 .0 , 

3530  &  13X,F8.0,2X,F8.1) 

3540  GO  TO  3590 

3550C 

3560  3558  WRITE(6 , 3559 )  SRUXAM(ITEM,  1 )  ,SRU.\’AM(ITEM,2)  ,SLLTC(ITEM,  1 ) , 

3570  &  SWUC(ITEM, 2), LOD, LOS, LOB, COSTL, ISOD, COST, ISOB, ISOC, UCS (ITEM) , 

3580  &  FPCT, FMTBCT, TLCDF 

3590  3559  FORMAT ( '  '  ,8X  ,2A4 ,2X , A4  ,A3  ,5X  ,A4 ,4X  ,A4 ,2X  ,A4  ,  IX  ,F7  .0  ,2X  ,A4 ,1X , 
3600  &  F8.0,3X,A4,2X,A4,F8.0,3X,F4.0, 

3610  &  13X ,F8 .0 ,2X ,F8 . 1 ) 

3620  GO  TO  3590 

3630C 

3640C  SRU  IS  BASE  REPAIRED 

3650  3570  COST-BASESC( ITEM) 

3660  TOTS5-TOTSR+COST 

3670  OPDECS (ITEM) -LOCAT ( 3) 

3680  GOTO  (3572,3574,3576,3573) .LEVEL 


3690  3372  VRTTK(6,3573)  SRI'NAM  ( ITEM ,  1  )  .SRUSAM (ITEM, 2)  .SV'Cl  ITEM ,  1  )  , 

3700  &  S'..>VC(  ITEM, 2)  .COSTL  ,LOS  ,  LOB ,  LOC ,  I  SOD  ,  ISOS  .COST  .ISOC  ,UCS  ( ITEM) 

3710  &,  FPCT, FMTBCT, TLCDF 

3720  3573  FORMATS  ' ,3X ,2A4 ,2X,A4 ,A3 ,2X ,FS .0 ,3X,A4 ,4X,A4 ,2X ,A4 ,2X ,A4 ,4X , 
3730  &  A4, IX, F3.0, 'ERROR' , F8 .0 , 2X , F5 . 1 , 

3740  &  13X.FR. 0.2X.F8.1) 

3750  GO  TO  3590 

3760  3574  URITF.(6,3575)  SRUMAM( ITEM,  1 )  .SRUSAM ( ITEM , 2 )  .SW.'CI ITEM.  1 )  , 

3770  &  SWi;C(ITEM,2) , LOO, COSTL, LOB, LOC, ISOO, ISOS  .COST, ISOC .KCS  ( ITEM) 

3730  &  FPCT.FMT3CT.TLCDF 

3790  3575  FORMAT ( '  ' ,3X ,2A4 .2X.A4 ,A3 ,5X ,A4 , IX ,F3 .0 , 3X.A4 ,2X , A4 ,2X ,A4  ,4X , 
3800  &  A4, IX, F8.0, 'ERROR' ,  F8.0 , 2X ,  F5 . 1 , 

3810  &  13X,F8.0,2X,F3.1) 

3820  GO  TO  3590 

3830  3576  WRITE(6,3577)  SRUNAM( ITEM, I ) .SRUNAM ( ITEM , 2) ,SWUC( ITEM, 1 ) , 

3840  &  SWUC(ITEM,2) .LOD.LOS, COSTL, LOC, ISOO, ISOS, COST, ISOC, UCS(ITEM) 

3850  &  FPCT , FMTBCT ,TLCD  F 

3860  3577  FORMAT ( '  ' ,8X,2A4,2X,A4,A3,5X,A4,4X,A4,1X,F8.0,1X,A4,2X,A4,4X, 
3870  &  A4,1X,F8.0,1X,A4,F8.0,2X,F5.1, 

3830  &  13X,F8.0,2X,F8.1) 

3890  GOTO  3590 

3900  3578  WRITE(6,3579)  SRUSAM ( ITEM, 1 ) .SRUMAM (ITEM, 2) ,SWTC(ITEM, 1 ) , 

3910  &  S WUC(ITEM, 2), LOD.LOS, LOB, COSTL, ISOD, ISOS, COST, ISOC, UCS( ITEM) 

3920  &  FPCT , FMT3CT ,TLC0F 

3930  3579  FORMAT ( '  ' ,8X, 2A4 ,2X ,A4 ,A3 , 5X ,A4 ,4X , A4 , 2X ,A4 ,F7 .0 ,2X ,A4 ,4X , 

3940  &  A4,1X,F8.0, 'ERROR' ,F3.0,2X,F5. 1 , 

3950  6,  13X,F8.0,2X,F8. 1 ) 

3960  GO  TO  3590 

3970C 

3980C  SRU  IS  CIRF  REPAIRED 

3990C 

4000  3680  COST*CIRFSC( ITEM) 

4010  TOTSC-TOTSC+COST 

4020  0PDECS(ITEM)«L0CAT(4) 

4030  GO  TO  (3682,3684,3636,3683) .LEVEL 

4040  3682  ’.TRITE  (6 ,3683)  SRUNAM(  ITEM,  l )  , SRUMAM  (ITEM, 2)  ,SWUC(ITEM,  l ) , 

4050  &  SWUCdTEM,  2),  COSTL,  LOS,  LOB,  LOC,  ISOD,  ISOS,  ISOB,  COST,  UCS(ITEM) 

4060  &  FPCT, FMTBCT, TLCOF 

4070  3683  FORMAT ('  ',RX,2A4,2X,A4,A3,2X,F8.0,3X,A4,4X,A4,2X,A4,2X,A4,4X, 
4030  &  A4 , IX, 'ERROR' ,F8.0,F8.0,2X,F5.l, 

4090  &  1 3X.F8.0 ,2X,F8 . I ) 

4100  GO  TO  3590 

4110  3684  WRITE(6,3685)  SRUNAM(ITEM,  l  )  ,  SRI'NAM  (ITEM,  2 )  ,SWUC(ITEM, 1 )  , 

4120  &  SWUCdTEIl,  2),  LOT),  COSTL,  LOB,  LOC,  ISOD,  ISOS,  ISOB,  COST,  UCS(ITRM) 

4130  &  FPCT, FMTBCT, TLCDF 

4140  3685  FORMAT ( '  ' ,SX,2A4 ,2X,A4 ,A3,5X ,A4 , IX, F8 .0 ,3X,A4 , 2X,A4,2X ,A4 ,4X , 
4150  &  A4, IX, 'ERROR', F8.0,F8.O,2X,F5.1, 

4160  &  1 3X ,F8.0 ,2X,F8 . I ) 

4170  GO  TO  3590 

4180  3686  VRTTE(6 ,3687)  SRUSAM  (ITEM,  1 )  ,  SRI’NAM  (ITEM,  2  )  ,S WUC (ITEM,  1 )  , 

4190  S.  SLTTC(TTEM,2) ,  LOD.LOS,  COSTL,  LOC,  I  SOD,  ISOS,  I  SOB,  COS  T.UCS  (ITEM) 

4200  &  FPCT,  FMTBCT, TLCDF 

4210  3637  FORMAT ( '  ' ,8X,2A4 ,2X ,A4 ,A3 ,5X ,A4 ,4X , A4 , IX ,F8.0 , IX ,A4 ,2X ,A4 ,4X , 
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4220  &  A4 , 2 X  ,  A4  ,  F 7  •  0 , 1 X  ,  F8 . 0 , 2 X  ,  F 5  . 1  , 

4230  &  1)X,F8.0,2X,F8.1) 

4240  GOTO  3590 

4250  3688  VRITF.(6,3639)  RRUNAM(ITEM.l)  ,SRUNA>1(ITF.M,2)  ,S WUC(ITEM , l )  , 

4260  f.  SUnC(ITEM,2)  ,LOD , LOS  , LOB  ,COSTL ,  ISOD  , ISOS  ,IS03, COST  ,l'CS  ( ITEM) 

4270  &  fpct,fmtbct,tlcdf 

4280  3639  FOR'lAT('  ' ,8X,2A4 ,2X ,A4 ,A3 , 5X ,A4 , 4X , A4 , 2X ,A4 , F7 .0 ,2X , A4 ,4X , 

4290  &  A4, IX, 'ERROR', FS .0 , F8 .0 , 2X , F5 . 1 , 

4300  6,  13X,F8.0,2X,F8.1) 

43  IOC 

4320  3590  LIMES=LINES+1 
4330  3600  CONTINUE 
4340C 

4350C  WRITE  COST  TOTALS 

4360  WRITE (6, 3604) 

4370  3604  FORMAT ( '  ',39X,' -  LC  S/BASE  - ') 

4380  WRITE(6,3605) 

4390  3605  FORMAT('0' , 35X ,7X, 'DEPOT' , 7X, ' SCRAP' ,8X,'BASF/  ,7X,'CIRF' ,14X, 
4400  &  'WHOLESALE  CHANGE  FACTORS') 

4410  WRITE(6,3610)  TOTLO.TOTLS  ,TOTLR ,TOTLC , WCFUC 

4420  3610  FORMAT ('  TOTAL  LRU  COSTS ', 20X ,4F1 2 .0 , 19X ,' COST  «',F6.2) 

4430  WRITE(6,3620)  TOTSD ,TOTS5 ,TOTSB ,TOTSC , WCFF 

4440  3620  FORMAT ('  TOTAL  SRU  COSTS ' ,20X ,4F1 2 .0 , 19X, 'MTBF  »',F6.2) 

4450  WRIT£(6f 3630)  TOTSED,TOTSEB,TOTSEC,UCFSE 

4460  3630  FORMAT ( '  ', 'TOTAL  SE  COSTS' , 2 IX ,Fl 2 .0 , 1 2X ,F1 2 .0 , 1 2X , FI 2 .0 , 

4470  &  9X,'SE  *' , F6 . 2 ) 

4480  DEV-DEV/M 

4490  '.TRITE (6, 3635)  DEV 

4500  3635  FORMAT ('  ','SE  DEVELOPMENT  COST' ,24X,F12.0) 

4510  GTOT-TOTLD+TOTLS+TOTL3+TOTSD+TOTSS+TOTSB+TOTSED+TOTSF.B+DEV 

4520  &  +TOTLC+TOTSEC+TOTSC 

4530  WRITE(6,3640)  CTOT 

4540  3640  FORMAT (' 0' ,5X,' TOTAL  COST  OF  REPAIR  LEVEL  DECISIONS  «' , 

4550  S,  F12 .0/20  ('  *' ) ) 

4550  RETURN 
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Appendix  H:  SECMP  Subroutine 


10  subroutine  secmp(ind.mscli) 

20C 

30C 

40  I F  ( I ND .  EO .  0 )  CO  TO  306 

50  IPAGF.=  IPAGE+l 

60  800  WRITF,(6,801)  (  SYSNAM( ISUB)  ,ISUB=1 ,3)  ,DATE  ,  I  PACE 

70  301  FORMAT  ('l',2X,3A4,86X,A8,5X,'PAOF.',I4) 

80  WRITE(6 ,802)  CHARS (2) , CHARS (2) 

90  802  FORMAT ( *  ' ,43X , 'COMPUTED  SE  COSTS ', 2X , 2A4 ) 

100  WRITE (6 ,803)  PIUP 

110  803  FORMAT ('O' ,4X,'SE',9X,'SE' ,9X,'***  NO.  ***  ***  SE  USE  *** 

120  &  '  '  ,3X,  'CODE'  ,7X,  'NAME',  8X/CURR  NEW'  ,8X,'HR3/MO  ',4X 

130  &  '%USE'  ,15X,F4.0,'  YR' ,  14X, 'TOTAL  SF/) 

140  WRITE(6 ,804) 

150  804  FORMAT ('  '  , 26X , '  IN'V  REQD  CURR  NEW'  ,  12X  , '  ACQ .  S' ,4X 

160  &  'OPER.  $  '  ,  5X, '  FAC.  $',2X,'LC  S/BASF.') 

170  W?ITE(6,805) 

180  805  FORMAT ( '  DEPOT') 

190  LtNES-7 

200  IF(NDSE.CT.O)  CO  TO  S06 

210  WRITE (6 ,8050) 

220  8050  FORMAT ( '  BASE') 

230  LINES=LINES+l 

240  IF  (NBSE.GT.O)  GOTO  306 

250  WRITE ( 6 ,8051 ) 

260  3051  FORMAT ('  CIRF') 

270C 

280  806  IF(N’UMSER.LT.l)  RETURN 

290  DO  900  1900-1 , MUMS ER 

300  SAVCST«SECOST(I900) 

310C  COMPUTE  TOTAL  USAGE  HOURS  FOR  THE  SE 

320  ITYPE*SECODE(I900) /10000 

330  JTYPE«(SECODE(I900)/1000)-(ITYPE*10) 

340  TMURS*U3EHRS ( 1900) 

350  TMROT*REQMT ( 1900 ) 

360  USEHRS ( I 900 ) *0 . 

370  NEXT“SEPF(I900) 

380  IF(NEXT.EO.O)  CO  TO  815 

390  LAS  T=S  EPL (1900) 

400  807  ITMPTR-SEXREF(NF.XT) 

410  I F ( ITMPTR.CT .0 . AND .OPDECL( ITMPTR) .NE.LOCAT (15))  GOTO  311 

420  IF( ITMPTR. LT.O. AND .OPD ECS (0- ITMPTR) .NE.LOCAT (15))  GOTO  811 

430  IF(ITYPE.GT. 4. AND.ITYPE.LT. 9)  GO  TO  809 

440  IF(ITYPE.EQ.9)  GOTO  8090 

450  IF (ITMPTR. LT.O)  GO  TO  808 

460  US  EHRS ( 1 900 ) *US  EHRS ( 1900 ) -*TS  EUHD ( ITM°TR ) 

470  GO  TO  811 

480  803  ITMPTR-O- ITMPTR 

490  USEHRS (1900) *US EHRS (1900) +SSEUHD( ITMPTR) 

500  GO  TO  811 

510  809  I F( ITMPTR.LT .0)  GO  TO  810 
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520  USEHRS  ( I90Q)=USEHRS  ( I?03)+FSF"MB(  I  IMPTR) 

530  GO  TO  Sll 

540  3090  IF(ITMPTR.LT.O)  GOTO  8095 

550  USEHRS  (I900)=!JSE!IRS  ( I900)+FSEUHC ( ITMPTR ) 

560  GOTO  Sll 

570  3095  ITMnTR*0— ITMPTR 

5S0  USEHRS  ( 1900  )=U5EHRS(  1900  )+SSEUHC(TTrPTP.) 

590  GOTO  811 

600  810  ITMPTR*0-ITMPTR 

610  USEHRS  ( 1 900)  =USEHRS  ( 1900) +S  SFl'H  B  (IT'PTR ) 

620  811  IF (NEXT. EO.LAST)  GO  TO  314 

630  NEXT* NX TIT M (NEXT) 

640  IF (NEXT .GT.O)  CO  TO  807 

650  812  WRIT£(6,813)  1900, SEPF(I900) ,SEPL(I900) .NEXT, SEXREF(NCXT) , 

660  &  NXTIT.M(NEXT) 

670  813  FORMAT ('OTROUBLE  IN  SEXREF',616) 

680  STOP 

690  814  IF(NXTITM(NEXT) .EQ.O)  GO  TO  815 

700  WRITE(6,S140) 

710  3140  F0R.‘1AT ( '  NEXT  =  O') 

720  GO  TO  812 

730C  COMPUTE  If  OF  UNITS  OF  THE  SE  REQUIRED 

740  315  IF(ITYPE.E0.1 .OR . ITYPE.EQ.5 .OR. JTYPE.EQ. 1 )  GO  TO  317 

750C  FOR  PECULIAR  SE 

760  IF  (JTYPE.GT.l)  GOTO  816 

770  IF  (USEHRS  (1900). CT.’.IRAVSE(I9O0))  GOTO  318 

780  RF,OMT(I900)*0 

790  GOTO  819 

800  816  TOTHRS-T>P»RS+USEHRS(I900) 

810  REOMT (1900)*  AIN'T (  2 .+(  (TOTHRS *OS  )  /OPHRS  (1900)  )  )  + 

820  &  AIN'T  (  (TOTHRS* ( 1  .-OS  )  )  /OPHRS  ( 1900) ) 

830  IF(TMRQT.CT.O)  REQ:iT(I900)=REQMT(I900)-TMRQT 

840  GOTO  819 

850  318  RF.QMT ( 1900)*  AINT (  l  +  (  (USEHRS  ( 1 900)-HRAVSF.(I900) )  /OPHRS  ( 1900) )  ) 

860  31°  S EUR ( 1900 )* (USEHRS ( I900)+T‘THRS ) / ( (REOMT ( I900)+TMR0T) *OPHRS (1900 ) ) 

870  IF(ITYPE.NE.4.AND .ITYPE.NE.8 .AND . JTYPE.NE.4)  CO  TO  840 

880  REQMT(I900)=1 

890  Tf!ROT=0 

900  SEUR(I900)-0. 

910  GO  TO  840 

920C  FOR  COMMON  SE 

930  317  REOMT (I 900 )*0. 

940  IF  (JTYPE.EO.l)  GOTO  821 

950  AVHRS*(OPHRS (1900 )-3SYHRS (1900) )*NSECI ( 1900 )-TMHR3 

960  820  IF(AVHRS.GT. USEHRS (1930))  GO  TO  830 

970  REQMT(I900)*REQMT(I900)+1 . 

980  AVHRS*AVHRS+OPHRS (1900) 

990  GO  TO  820 

1000  321  TOTHRS*T!rnRS-HjSF.HRS  ( 1900) +BSY HRS  (1900) 

10  in  REQMT(l900)-2.+AT'IT((T0THRS*nS)/o?HRS(I"00) )  + 

1029  6  AINT ( (TOTHRS* ( 1 .-OS) ) /OPT PS ( 1900 ) ) 

103Q  REQMT(I900)=REQMT(I900)-T;ROT-"SFCI(I900) 


s\'. 


i*a 


1040 

1050c 

1060 

1070 

1030C 

1090C 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180C 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1230 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1330 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1430 

1490 

1500 

1510 

1520 

1530 

1540C 


on to  330 


830  SEUR( 1900 ) = (USEHRS ( I900)+\’SECI  (1900 ) *BSYHRS  (T900)+TMMRS )  / 

&  (NSF.CI  (1900) *OPHRS(  1900 )  +  (RF.QMT (1900) +-TMT>nT)*0?UR8(l9Q9)) 

COMMUTE  SE  ACOUISITrON  &  OPERATION'S  COST 
840  SEACQaREOMT ( 1900) *CAf)R ( 1900) 

SEOP:,'3-USE‘!RS(I90O)/OPHRS(I9OO)*CODB(I9OO)*PIl'P 
SEFAC=FDK (1900) 

IF(TMROT.CT.O)  SEFAC=0 
SECOS  r ( 1900 ) =S EACO+S  EOPN5+S  EFAC 
IF(ITYPF..CT. 4. AND. ITYPE.LT. 8)  CO  TO  350 
IF  (ITYPE.EQ.9.AN0.JTYPC.E0.4)  COTO  350 
SECOST (I900)*SECOST (1900 ) /M 
SAVE  TOTAL  COST  FOR  THE  NET CORK  SE  COST  ARC 
850  SECOST  ( 1900)  =AIN'T (S  FCOST  ( 1900  )+0 . 5 ) 

IF(IND.EQ.O)  CO  TO  330 
I F(LIN'ES .LT . 55 )  CO  TO  355 
I P  AC  E=*  IPAG  E+l 

WRITE(6,80l)  (SYSNAM(ISUB)  , ISUB*l  ,3)  ,DATE,IPAGE 
WRITE(6,802)  CHARS  (5)  .CHARS  (6) 

WRITE (6 ,803)  PIUP 
WE  ITE ( 6 ,804 ) 

LINES  =*7 
355  choi:rs=»o. 

IFdTYPE.EO.  1  .0R.ITYPE.E0.5 .OR.JTYPE.~0.  I ) 

&  CPOURSmMSECT  ( I9O0)*ESY’TRS  (1900) 

REOMT(I900)-REOyT(n00)+T'fRQT 
USEHP.S  (1900  )*USE HRS  (  I900)+THHRS 
S  EACQ-REQMT  (1900)  *CA0  P,  (1900 ) 

SEOPNS* (USEMRS (1900) /OPHR5 (1900) )*CODB(I903) *PIUP 
S  EFAC*FD  3( 1900 ) 

TS  ECST  ( 1900)  =TSECST(Iori9)+S  ECOST  (1900 ) 

TSECST  (I 900 )=AI NT (TSECST (1900) +0.5) 

WRITE  (6, 860)  SEC0DE(T900) ,  (XSE(I990  ,  ISUR)  ,  ISUB*1  ,3)  ,N’SECI(I90n)  , 
&  REQMT  ( 1900 )  .CHOURS .USEHRS  (1900)  ,SEUR(I900)  ,SEACO  .SF.OPNS  , 

4  SF.FAC  .TSECST  (1900) 

860  FORMAT ( '  ' , 3X, 15 ,3X , 3A4 , 3X , 14 ,4X , 14 ,2X ,2F7 . 1 ,2PF6 . 1 ,0P4Fl 1 .0) 
LINES-LINFS+1 

IF(I900.LE.NT)SF.)  GO  TO  330 
IF(T900.CT.r:0SE)  COTO  375 
WRTTE(6,870) 

370  FORMAT ('  BASE') 

LINES*LI?!ES+1 
GOTO  380 

375  IF(I900.NF..:JDSE+‘!BSE)  GOTO  380 
1/',ITF.(6 ,877 ) 

377  FORMAT ( '  CIRF') 

LIMES-LINES+1 
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